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SYSTEMS IDENTIFICATION USING 

A MODIFIED NEWTON-RAPHSON METHOD - A FORTRAN PROGRAM 


By Lawrence W. Taylor, Jr. 

Langley Research Center 

and 

Kenneth W. Iliff 
Flight Research Center 

SUMMARY 

A digital computer program written in FORTRAN is offered which computes a 
maximum-likelihood estimate of the parameters of a linear, state space model. For the 
case considered, the maximum-likelihood estimate can be identical to one which minimizes 
simultaneously the weighted mean- square difference between the computed and measured 
response of a system, and the weighted square of the difference between the estimated and 
a priori parameter values. A modified Newton-Raphson or quasilinearization method is 
used to perform the minimization which typically requires several iterations. The modi- 
fication of the Newton-Raphson method was made in the interest of reduced computation 
and program simplicity. A starting technique is used which insures convergence for any 
initial values of the unknown parameters. A Cramer-Rao bound is used to indicate the 
variance of the estimated parameter values. Although the primary application of the pro- 
gram has been to determine aircraft stability derivatives from flight data, it is directly 
applicable to identification of any system which can be described by a linear, constant- 
coefficient model. The intent of this paper is to describe the program and its operation 
in sufficient detail to enable the user to apply the program to his particular problem with 
a minimum of difficulty. 


INTRODUCTION 

Problems in systems identification, in which models are constructed from response 
data, continue to receive considerable attention because of their importance and because 
of the difficulties involved. A few years ago linear regression methods were employed 
almost exclusively for modeling systems. The most often used methods were least 
squares (ref. 1) and Shinbrot's method (ref. 2). Unfortunately, these methods often 
required measurements of variables not generally available and, at best, would match the 



differential equations of motion to the response data. Too often the response calculated 
from the resulting models did not match the measured response well. It was recognized 
that it would be better to match the model response to the measured response per se, but 
that problem was difficult to solve because it was nonlinear. Only after experience had 
been gained in solving large numbers of nonlinear, simultaneous equations were these sys- 
tems identification problems solved by using a modified Newton-Raphson method (ref. 3) 
(or quasilinearization (ref. 4)). A second order was used because of its superior conver- 
gence properties. Several algorithms (refs. 5 to 9) have now been developed and numer- 
ous applications have been made (refs. 10 and 11, for example), particularly to the prob- 
lem of determining stability derivatives of airplanes. 

The modified Newton-Raphson (quasilinearization) method is applicable for model- 
ing any linear system with constant coefficients. The algorithm includes features which 
enable the use of a priori information such as wind-tunnel measurements. The initial 
conditions and bias terms can also be treated as unknown parameters. 

The purpose of this paper is to make available in a form for general use the pro- 
gram that has been successfully used for determination of aerodynamic derivatives in 
references 10 and 11. It is intended that the accompanying description, program, and 
example test case would enable a prospective user to run the test case first and then sub- 
stitute the data for his particular application with relative ease. 

SYMBOLS 

Program variables are given in a separate section, 
systems matrix (MX by MX) 
control matrix (MX by MU) 
parameter vector (MC by 1) 
a priori parameter vector (MC by 1) 
estimate of parameter vector (MC by 1) 
weighting matrix for measured response 
weighting matrix for a priori estimate vector 
expected value 


A 

B 

c 

c 0 

c 

Dl 

E{> 
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state transformation matrix (MZ by MX) 
control transformation matrix (MZ by MU) 
time increment 
identity matrix 

time and iteration indices, respectively 

cost function to be minimized 

number of time points 

probability 

roll rate 

yaw rate 

record length, sec; as superscript, transpose 
time, sec 
control vector 
state vector 

time derivative of state vector 
calculated response vector 
measured response vector 
sideslip angle 
aileron deflection 


rudder deflection 



<p 


bank angle 


V c ( ) first variation (gradient) with respect to parameter vector 

V c 2 ( ) second variation (gradient) with respect to parameter vector 

Dots over symbols denote derivatives with respect to time. 

The program variables used in the definitions of the symbols are as follows: 

MC number of unknown parameters 

MU number of control variables 

MX number of state variables 

MZ number of response variables 

SYSTEMS IDENTIFICATION 

The systems identification problem of determining the parameters of a linear, 
constant-coefficient model will be considered for several types of estimates. It will be 
shown that maximum-likelihood estimates can be identical to those which minimize the 
mean-square response error. The subject algorithm, therefore, can be used to obtain a 
variety of similar estimates. Attention is also given to the calculation of the gradient 
that is involved in the algorithm and to the Cramer-Rao bound which indicates the vari- 
ance of the estimates. 


Problem Statement 

The problem considered is that of determining the values of certain model param- 
eters which are best with regard to a particular criterion, if the input and noise measure- 
ments of the response of a linear, constant- coefficient system are given. The system to 
be identified is defined by the following equations: 


x = Ax + Bu 

(1) 

y = Fx + Gu + b 

(2) 

z = y + n 

(3) 


4 


where 


x state vector 

u control vector 

y calculated response variable 

b constant-bias vector 

n noise vector 

z measured response variable 

The unknown parameters will form a vector c. The matrices A, B, F, and G and the 
vectors b and x(0) are functions of c. 

Minimum-Response Error Estimate 

One criterion that is often used in systems identification is the mean-square differ- 
ence between the measured response and that given by the model. A cost function which 
is proportional to the mean- square error can be written as 

N 

J = I ( z i - yif D i( z i - yi) 

i=l 

where Dl is a weighting matrix and i is a time index. The summation approximates 
a time integral. The estimate of c is then 

c = ARG MIN(J) (5) 

c 

where ARG MIN means that vector c which minimizes the cost function J. 

Linearize the calculated response y with respect to the unknown parameter vec- 
tor c: 

Yi = y iQ + V c Yi(c - c 0 ) (6) 
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where 


y i 0 

nominal response calculated by using Cq 

v c yi 

gradient of 

y with respect to c 

c 0 

nominal c 

vector 


Substituting y i into the expression for J and solving for the value of c which 
minimizes J yields 



— — l 

N 

V ' 

-1 

N 

'T ' 

c = c 0 + 

2 v cyi DiVcYi 
i=l 

L J 


Z Voy, D x( z i - y io ) 


If this relationship is applied iteratively to update the calculated nominal response 
and its gradient with respect to the unknown parameter vector, the minimum-response 
error estimate c will result. The method has been called quasilinearization (refs. 6 
and 7) and modified Newton-Raphson (refs. 5 and 10). The latter seems more appropriate 
since the Newton-Raphson (ref. 3) method would give 


where 


c k + l = c k + 


— I — 1 i — 


v c 2 Jk 


v c J k 


(8) 


N 

V c J k = -2 y V c y ik T Di^ - y ik ) 


(9) 


N N 

V c 2 J k = 2 y V c yj J^Vcy, + 2 £ - y. 


i=l 


i=l 


(10) 


The second term of diminishes as the response error - y^ diminishes. The 

modified Newton-Raphson method is identical to quasilinearization if the term is neglected. 
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Maximum Conditional Likelihood Estimate 


Another criterion that is often used is to select c to maximize the likelihood of 
the measured response when c is given: 

c = ARG MAX {P(Z |c)} (11) 

c 

If the noise n has a Gaussian distribution, the probability distribution of the mea- 
sured response at any single time point can be written as 



1_ 

(MZ)/2 

(2tt) 





( 12 ) 


where 

Mj = E 

and MZ is the number of elements in the response vector. 

If the values of n are uncorrelated at different time samples, the conditional 
probability (given the value of c) of the entire set of measured response z can be writ- 
ten as 





N 

P(z |c) = 

N(MZ)/2 

v2t t) 

Mj 

exp 

N/2 

'i 1 ^i' y i) T “ 1 " 1 ( I i' y i) 


(13) 


The maximum conditional likelihood estimate of the unknown parameters will be the 
set of values of c which maximize P(z|c), if it is recognized that y is a function of c. 
If it is noted that the maximization of P(z|c) occurs for the value of c which minimizes 
the exponent, the maximum conditional likelihood estimate is the same as that which mini- 
mizes the mean-square response error provided the weighting matrix Dj equals Mj - ^. 

Maximum Unconditional Likelihood (Bayesian) Estimate 
The unconditional probability of z can be expressed as 

P(z) = P(z|c)P(c) (14) 
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The probability of c relates to the a priori information available for c before 
use is made of the measurements z. If P(c) is Gaussian, 


P(c) = 


(2 n) 


(MC)/2, ,1/2 


exp 




-f( c - c o) Tm 2 _1 ( c " c o) 


where 


M 2 = E {( c - c o)( c - c o) r 


(15) 


c 0 


= E{c) 


The unconditional probability of z is then 

N 


P(z) = - 


1 


[(MC)+N(MZ)]/2 | N/2. .1/2 

' 'm, m , 1 


exp 


(2tt) 


- \ l{ z i - n) TM i~\ z i - h) - |( c - c o) Tm 2 _1 ( c - c o) 

i=l 


(16) 


Again the expression is maximized by minimizing the exponent. Setting the gradient 
with respect to c equal to zero and solving yields 



r H 

-1 

f— — 1 


N 


N 

c k+l " c k + 

J, *cV I ' M l' lv c y i k + M 2' 1 




i_l t. 


i=l 


(17) 


An identical estimate would have resulted if a weighted sum of mean-square 
response error and the mean-square difference of c and its a priori value are mini- 
mized, provided the weighting matrices used equaled the appropriate inverse error 
covariance matrices. Consequently, the same algorithm can be used for both estimates 
(ref. 11): 


c 





(18) 
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Computation of the Response Variable and Its Gradient 

It is necessary first to solve the state equation (1) for the state variable x. The 
method used is that given by the recursive formula: 

x. +1 = (PHDxj + (APHI)(B) ^ + - Ul+1 (19) 


where 


PHI = e^ = I + Ah + + — — + . . . 

2: 3! 


APHI 


J 0 


A u2 a 2v,3 a 3 k 4 
e At dt = Ih + — + AJl_ + A JL. + 
2 : 3: 4: 


Through repeated application of the recursive formula and with the initial state and the 
control variable u known, it is possible to compute the state variable x. With x 
known, the response variable y can be computed, in turn, by using equation (2). 

The method used to compute the gradient of the response variable y with respect 
to the unknown parameter vector c is based on an equation obtained by differentiating 
the state equation (1) with respect to c 


9x _ ^ 9x 
9c 9c 




(20) 


Because of the similarity of the equation and the state equation (1), it is evident that 
it can be solved with the recursive formula 


(ft) =< PHI >(ft) + <APHI)M(Xi + x i+1 ) + ffiV^l ± l (21) 

' 'i+1 ' l+l 


Repeated application of the recursive formula will yield the partial derivative of the state 
variable with respect to c. The equation for the gradient of the response variable y 
with respect to c is obtained by differentiating equation (2) 


£y 

9c 



9F 

+ x 


9c 



9b 

4 

9C 


( 22 ) 
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It is often the case that many of the state -variable elements are measured directly. 
In that event, it is very useful to compute the gradient for the first iteration by using the 
measured-state variable instead of the computed- state variable. This technique usually 
removes the need for initial estimates of the unknown parameter vector c to such a 
degree that a null vector for c can be used initially. 

Variance of the Estimates 

A very important aspect of systems identification is the quality of the estimates. 
Since the estimates themselves can be considered to be random numbers, it is useful to 
consider their variance, or specifically their error covariance matrix. The Cramer-Rao 
bound (ref. 3) provides an estimate of the error covariance matrix and is identical to that 
developed in the following discussion. 

Note that when the maximum unconditional likelihood estimate of c converges, it 
is necessary for the following equation to be satisfied: 

N 

X v cy 1 TM i' 1 (* i - y s ) - M 2 _1 ( e - c o) = 0 <23 > 

i=l 


If c true is introduced, and y^ is linearized with respect to c, it follows that 



After solving for c - c^ rue , 

N 

1 v cy i TM f 1 (n i ) - M 2 "Vtrue - c o) 
i=l 

where 



Q 


N 

^ v cyi TM i" lv cyi + m 2 

i=l 


°i z i ^i,true 
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Squaring and taking the expected value yields 



C ' 


N 


N 

( £ - c true)( e - c true) 

}=Q' 1 

2 v c yi T Mf ^cYi + M^ 1 

Q' 1 = Q' 1 = 

1 V c y i T M 1 - 1 V c y i + M 2 - 1 



i=i 


i=l 


which is the inverse that is required in obtaining the unconditional maximum-likelihood 
estimate of c. If the a priori information is not used, M2 - * will be the null matrix so 
that 



The resulting estimate of the parameter error covariance matrix is useful in assessing 
the quality of the estimated parameter values. 
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PROGRAM DESCRIPTION 


The purpose of this digital computer program is to determine by means of a modi- 
fied Newton-Raphson method the unknown parameters of a model used to describe a sys- 
tem. In order to accomplish this purpose, it is necessary to input the data, control infor- 
mation, starting values, and so forth. Recursive formula (17) represents the method on 
which the algorithm is based. Because of the large number of values involved, matrices 
are used throughout the program. In order to conserve memory and to facilitate changing 
the sizes of the arrays, a vector format has been used for the matrices. The first and 
second elements of the vector are the number of rows and columns, respectively. The 
third and fourth elements of the vector are the number of elements and type, respectively. 
Starting with the fifth element of the vector the elements of the matrix are entered in 
order starting with the first row. 

For example, the matrix 

1.2 -3.7 0.5 

0.09 15.3 0.0_ 

would be represented by 

2., 3., 6., 1., 1.2, -3.7, 0.5, 0.09, 15.3, 0., 

Several subroutines are used to aid in the manipulation of the matrices the descrip- 
tions of which are given in appendix A. 

Program Variables 

The program variables are as follows: 

A stability matrix of dynamic equation (1) 

AP array of parameter information: (1) a priori value, (2) variance - -*-, and 

(3) scale factor 

ft At! 

APHI integral of transition matrix, \ e A dti 

J 0 

B control matrix of dynamic equation (1) 
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DC 


incremental change in parameter vector 


DUD,DUM dummy arrays 

DZN bias in transformation equation (2) 

DZ1 initial value of Y 

D1S weighting matrix for difference in response 

D1 row vector of diagonal elements of matrix D1S 

EE weighting matrix for difference in parameter values 

ENN number of time points 

F state transformation matrix in equation (2) 

FINK scaling parameter to avoid overflow and underflow during matrix inversions 

G control transformation matrix in equation (2) 

HH time increment 

I general index 

IG array identification (1) A, (2) B, (3) DZ1, (4) DZN, (5) F, (6) G 

II row index 

IJ vector element index corresponding to matrix element I,J 

IP array of indices which identify array and element for AP 

ITR number of iterations 

J general index 

JI vector element index corresponding to matrix element J,I 
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JJ 


column index 


JK parameter index 

JKM number of unknown parameters plus 1 

JKMM1 number of unknown parameters 

K general index 

LL iteration index 

LLM1 iteration index minus 1 

MU number of control elements 

MX number of state elements 

MXP1 = MX + 1 

MZ number of response elements 

NN number of time points 

NNM1 = NN - 1 

NP index of parameter information sets 

NPM number of parameter information sets 

PHI state transition matrix, 

SU scale factor for control 

SUM matrix of simultaneous equations 

SZ scale factor for response 

TT record length 

U control data 

U1 control vector 

XIA factor applied to diagonal of sum to aid inversion when sum is nearly singular 

X1,X2,XJ state vectors 
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XJI gradient of response variable Y with respect to parameter vector c 

Y calculated response 

Z measured response 


Example Problem 

An example problem is considered both to aid in explaining the operation of the pro- 
gram and to act as a test case. The problem is that of determining certain stability 
derivatives and biases with regard to the response of an airplane to a known input. The 
equations used to compute the state are the following (eq. (19)): 



r 

/3 






— | 

M 


— 


-j 

L P 

L r 

L 0 

0 


L *a 

L6 r 

LO 

N P 

N r 

N /3 

0 

f. 

+ 

N *a 

N 6 r 

N 0 

Y P 

-1 

Y /3 

Y 0 



Y 6a 

Y S r 

Y 0 

1 

0 

0 

0 

V/ 


0 

0 

0 



The response variables are given by (eq. (2)): 



1 

0 

0 

0 


r 

0 

0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

0 

1 

0 

/p\ 

0 

0 

0 

0 

0 

0 

1 

r 

4 - 

0 

0 

0 

L P 

L r 

L /3 

0 

P 


L S a 

L6 r 

LO 

N P 

N r 

N |3 

0 

\7 


N ^a 

N 6r 

N 0 

Y P 

-1 

Y i3 

Y <t> 


Y 5a 

Y a r 

Y o 



+ DZN 
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where L p , L r , L^, N p , N r , N^g, Y p , Y^, and Y<p are coefficients of the A and 
F matrices and L 6a> L 6r , L 0 , N 5a , N 6r , N 0 , Y Y 5r> and Y 0 are coefficients 

of the B and G matrices. The unknown parameter vector c is 



where Y^, Yg a , and Yg r are considered to be known. The initial values of the 

unknown parameters are chosen to be zero and hence the A, B, F, and G matrices will 
be null matrices (only for convenience) except for the elements which are known. The 
time histories of the control inputs and responses are shown in figure 1. Since the 
responses are noise free except for numerical errors, the parameter values reached by 
using the modified Newton-Raphson algorithm (ref. 5) would be expected to be very close 
to their true values. 
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Memory 


The total storage required for the program including its subroutines is 44 547g 
fixed-point words. This requirement will change, of course, as the dimensioned arrays 
are changed to meet the user's specific requirements. The following expressions are 
offered as an aid in changing the array dimensions. The dimensions can exceed those 
given, of course. 


A(MX • MX + 4) 

G(MZ • MU + 4) 

AP(NPM, 5) 

IP(NPM, 4) 

APHI(MX • MX + 4) 

PHI(MX • MX + 4) 

B(MX • MU + 4) 

SU(MU + 4) 

DC(JKM + 3) 

SZ(MZ + 4) 

DUD (MX • MX + 4) 

U(MU,NN) 

DUM(JKM • JKM + 4) 

U1(MU + 4) 

DZ1(MZ + 4) 

XJ(MZ + 4) 

DZN(MZ + 4) 

XJI(JKM • MZ + 4) 

D1(MZ + 4) 

X1(MZ + 4) 

D1S(MZ • MZ + 4) 

X2(MZ + 4) 

EE(MZ + 4) 

Y(MZ + 4) 

F(MZ • MX + 4) 

Z(MZ,NN) 


Input 

The entire input for the program is in the form of data cards. It is typical, how- 
ever, that the response data are recorded on magnetic tape. It will be necessary, there- 
fore, for the user to alter the program to fit his needs. The following description of the 
input data cards uses an example listed in appendix B. 

First card .- The first card contains two numbers: (1) NN, the number of data time 
points, and (2) ITR, the number of iterations. The first number must be right justified to 
column 10, the second to column 20. 

Second card .- The second card contains a single number corresponding to HH, the 
time increment of the periodically sampled response of the system. The number can be 
positioned anywhere within the first ten columns but must contain a decimal point. The 
value for the scaling parameter FINK is located in the second field and can be positioned 
anywhere within columns 11 to 20. FINK is used to avoid overflow and underflow during 
matrix inversion. 
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Matrix format (general).- The matrices and vectors are loaded by means of sub- 
routine load which calls for the following format. The first card contains (1) the matrix 
name in the first seven columns, (2) the number of rows right justified to column 10, and 
(3) the number of columns right justified to column 20. The next card contains the ele- 
ments of the first row of the matrix with one floating point number in each of seven fields 
of 10 columns. Additional cards are used as necessary to complete the first row. The 
same procedure is used for each succeeding row with each row starting in the first field 
of a card. 

Matrix A.- Matrix A is the systems matrix of the dynamic equation (1), the order 
of the coefficients can be any order that is consistent with the state variables. The order 
of the elements of the state vector for the example test case is p, r, j3, and <p. 

Matrix B.- The matrix B is the control matrix of the dynamic equations. The order 
of the coefficients of the B matrix must be consistent with that of the state and control 
variables. The order of the elements of the control vector for the example test case is 
6 a , S r , and 1. The constant 1 control variable is introduced to account for possible 
bias in the state. This bias is not to be confused with bias in the measured response. 

Matrix D1S.- The matrix D1S is the weighting matrix for differences in the response. 
For maximum -likelihood estimates, it should equal the inverse of the appropriate error 
covariance matrix. The matrix D1S can be a row vector of the diagonal elements or can 
be the complete square matrix. The program will recognize which you have chosen to 
use and will interpret D1S accordingly. 

Matrix F.- The matrix F is the lower part of the transformation matrix for the 
state that is not an identity matrix. If the particular problem does not involve responses 
which are not state variables, an F matrix is not needed. In this case, indicate zero rows 
and zero columns on the card for the F matrix. 

Matrix G.- The matrix G is the transformation matrix for the control. If the 
G matrix is not needed, indicate zero rows and zero columns on the card for the G matrix. 

Vector SZ.- The vector SZ contains the scale factors for the response variables. 

Vector SU.- The vector SU contains the scale factors for the control variables. 

Array IP,AP.- Array IP,AP is a set of both fixed and floating point numbers on a 
single card which contains the information which relates each unknown parameter to the 
element of each matrix or vector involved. 

(1) The first card of a set which pertains to a common parameter will have a 1 in 
column 10. Subsequent cards which pertain to the same parameter will have a 0 in col- 
umn 10 and must follow the first card. 
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(2) The number in column 20 of each card identifies each array in turn involving 
the parameter. The code follows: 

1 A matrix 


2 

3 

4 

5 

6 


B matrix 
DZ1 vector 
DZN vector 
F matrix 
G matrix 


(3) The number in column 30 of each card identifies the row of element number in 
the case of a vector. 

(4) The number in column 40 of each card identifies the column or is put to zero in 
the case of a vector. 

(5) The floating-point number in columns 41 to 50 is the a priori value of the param- 
eter. Only the first card for a particular parameter must have this information. 

(6) The floating-point number in columns 51 to 60 is the standard deviation of the a 
priori value of the parameter. Only the first card for a particular parameter must have 
this information. 

(7) The floating-point number in columns 61 to 70 is the scale or coefficient of the 
parameter for the particular element. 

The first card of the IP,AP array in the input listing (appendix B) is given in the 
label LP. The 1 in the first field indicates the card to be the first in a group pertaining 
to the parameter LP. The 1,1, and 1 in the second, third, and fourth fields indicate that 
the parameter is contained in the first matrix (A), first row, and first column. The large 
number 99999. in the next field says the variance of the a priori value of the parameter is 
infinite. The 1. in the next field indicates the scale or the ratio of the element to the 
parameter. If an a priori value was used, it would appear in the next field. 

The next card has the same label LP but has a 0 or blank in its field that indicates 
that it is concerned with the preceding parameter. The 5, 1, and 1 indicate that the 
parameter is contained in the fifth matrix (F), first row, and first column, and the scale 
factor is 1. again. The end of the array IP,AP is indicated by a card with a negative 
fixed-point number in the first field of 10 columns. 

Response and cont rol data .- The response data can be read by using any mode or 
format that is desired. For the example test case, the data are read by using a 6E12.4 
format for the variables 6 a , 6 r , 1, p, r, / 3 , <f), p, r, and / 3 . Two cards are 

required for each of the 200 time points. 
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Output 


The first line of output which is listed in appendix B indicates the number of data 
points, number of iterations, and the time increment between data points. Next, the A, 

B, Dl, F, G, SZ, and SU matrices are listed as part of the input data. The user should 
check the listing of the input for errors. The A and B matrices given in the examples 
have all the unknown parameters set to zero initially for convenience and to demonstrate 
how insensitive the algorithm is to initial values of the unknown parameters. 

The next block of output is IP,AP a table of the input cards used (1) to control which 
of the model parameters are to be adjusted, (2) to indicate which matrix elements are 
affected by changes in the parameter, (3) to indicate the relationship of the changes or 
scale, (4) to indicate the parameter's a priori value, and (5) to indicate its weighting. 

The next line gives the value of JKM, the number of rows and columns of a matrix 
of simultaneous equations. The number of unknown model parameters to be adjusted is 
one less. In the same line is given the value of NPM, the number of control cards of 
IP,AP just described. In the example each parameter affects two matrices (A and F or 
B and G) so that NPM equals twice the number of unknown parameters. 

The next block of output is a listing of the control and response data which is part 
of the input. The number of control and response variables is indicated and time is given 
in the first column. 

The next series of blocks of output indicate the results of the repeated application 
of the Newton-Raphson method in solving for the unknown parameters. 

First the cost or fit error is given. The cost must always be positive and should 
decrease and stabilize as the iterations proceed. Although it is necessary for the cost 
to behave in the manner just described, such behavior is not in itself sufficient evidence 
of a valid solution. The vector EE*D which indicates how the various components of the 
response error contribute to the total cost or fit error is given. 

The next vector given is the root mean square of the fit error in the same units as 
the corresponding response variable. Although these indices are useful, it is strongly 
recommended that the user look at the plotted time histories to assess the quality of the 
fit better and to detect erroneous data values. 

A vector of the partial derivatives of the cost with respect to each of the unknown 
parameters is also given. Theoretically, this vector should go to zero as a solution is 
reached. Practically, it is sufficient for the absolute value of its elements to continue to 
decrease as the iterations proceed. The value of the determinant of the matrix which is 
inverted is given in the next line. The vector DC gives the changes in the coefficients 
which should decrease in amplitude as a solution is approached. The last several blocks 
of output listed in appendix B are repeated for each iteration until the time history which 
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corresponds to the last coefficient values given is listed. Next is given the regular 
information for the last iteration followed by the A, B, DZ1, DZN, F, and G matrices. 

The last array in the output provides a summary of the results, including the final 
estimate and the approximate standard deviation if the inverse of the noise covariance 
matrix is used for D1S. The a priori values and a priori standard deviation are provided 
for comparison. 


Trouble Shooting 

Although there is no replacement for the investigator's understanding of his partic- 
ular problem, it is useful to point out some of the more common difficulties encountered 
in trying to obtain a solution. 

Case 1.- In case 1, the fit error converges as do most of the parameters but a few 
parameters diverge. Perhaps the most common difficulty in systems identification is 
caused by trying to determine a greater number of parameters than are justified by the 
amount and quality of the measured data. The only solution is to (1) increase the infor- 
mation processed by adding to the data base or adding a priori information, (2) ask for 
fewer parameters to be identified, or (3) improve the quality of the data by increasing the 
signal-noise ratio or by using more suitable inputs which reduce the correlation between 
variables. Determine that the same parameter does not occur twice in the control cards. 

Case 2.- In case 2 an error message, "Matrices Incompatible for MULT, DUM- 
DUM . . ."is listed repeatedly. The problem is just as the error message says and is 
caused by inconsistent input matrices. The number of states MX is determined by the 
number of rows in A. Matrix B must have the same number of rows; matrices A and F 
must have MX columns. The number of controls MU is determined by the number of 
columns in B. Matrix G must have the same number of columns. The number of 
responses MZ is determined by the number of columns in Dl. Matrices F and G must 
have MZ - MX rows. The matrices listed immediately after the error message are 
those for which multiplication was attempted. 

Case 3.- In case 3, the fit error neither converges nor diverges but just drifts. 
Something is wrong. Check to see that the data is in the proper order and that the 
parameter control cards are correct. Also, a model form which is in error could cause 
this type of behavior. 

Case 4.- In case 4, the fit error converges to a point and then fluctuates wildly. In 
certain marginal cases the numerical errors associated with inverting a matrix which is 
nearly singular can cause the fit error and parameter values to fluctuate. This problem 
can be alleviated by (1) reducing the number of unknown parameters, (2) adding a priori 
information, or (3) increasing slightly the factor multiplying the diagonal of the SUM 
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matrix. Care must be taken if the third option is chosen since a factor too large will 
make the algorithm appear to converge when it is simply sluggish. 

Case 5.- In case 5, the fit-error parameter values converge and the gradient con- 
tinues to diminish. These characteristics indicate that a bona fide solution has been 
reached. Another indication of the quality of the estimate of the parameter is its stan- 
dard deviation. If the standard deviation of the parameter is not noticeably smaller than 
its original value, the additional information contained in the measured data did little to 
improve the estimate of its value. 


Program NEWTON 

The program NEWTON performs system identification in state space format by 
using a modified Newton-Raphson method. Based on the recursive relationship, 
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The flow chart follows: 







The listing for program NEWTON is 


PROGRAM NEWTON ( INPUT , OUTPUT ,TAPE5= I 'iPUT , TAPE6=0UTPUT ) 

C SYSTEM IDENTIFICATION PROGR AM. . . .STAB IL ITY DERIVATIVES FROM FLIGHT DATA.. 

DIMENSION Z( 7(400 ) , U(5,400), DZ1 ( 12) , DZN< 12 ) , I P C 60 ,4 ) , API60.5) 

1 1 EE ( 12 ) t DUD (68 ) » APHK68), X1I16I, X2(16), XJ(16), Ul(12), 

2PHU68), 01(12), A (6 8 ) , B(36), XJK733), SUM ( 733) ,DC( 30 ) , DUM( 733 ) 

3, F (60) , G(60) , SZ125), SUI25), Y( 16 ) , D1S ( t 8 ) 

111 FORMAT ( 7X » 13,31 10.5F10.4) 

606 FORMAT (6E12. 4) 

700 FORMAT( 8F1C.4) 

777 FORMAT (7110) 

1001 FORMATdlH THERE ARE ,I4,14H DATA POINTS, ,I2,38H ITERATIONS AND T 
1HE TIME INCREMENT IS ,F5.4,8H SECONDS) 

1002 FORMAT! 80H PARAMETER GROUP ROW COLUMN APRIORI WEIG 

1HT SCALE ANSWER > 

1003 FORMAT! 6H JKM *,I3,7H, NPM =,I3) 

1004 FORMAT ( 36H THE FIRST COLUMN IS TIME. THE NEXT ,I2,41H COLUMNS ARE 
1C0NTR0L VARIABLES. THE NEXT ,I2,32H COLUMNS ARE RESPONSE VARIABLE 

2S . ) 

1005 FORMAT ( 31H THE VALUE OF THE FIT ERROR IS ,E12.5,26H FOR THE HERAT 
1I0N NUMBER ,13) 

1006 FORMAT ( 39H THE SCALING PARAMETER, FINK*S VALUE IS.F10.4) 

1010 F0RMAT(F8.2,10E12.4) 

CIIIIIHIIIIIIHIIIIIIIIIIUninil INPUT DATA (START) * 

3 READ 777, NN, ITR 
READ 700, FH, FINK 
PRINT 1001, NN, ITR, HH 
PRINT 1006, FINK 
ENN=NN 

CALL LOAD (4, A, B, D1S, F) 

CALL LOAD ( 3 , G, SZ, SU, SU) 

CALL SPIT ( A , 1HA ) 

CALL SP IT ( B , 1HB) 

CALL SP IT ( CIS, 3HD1S) 

CALL SPIT (F , 1HF ) 

CALL SP IT ( G , 1HG) 

CALL SP IT ( SZ , 2HSZ) 

CALL SPITISU. 2HS0) 

NNM1 = NN - 1 
MU » B(2) ♦ .01 
MX = A ( 2 ) ■* .01 
MZ = D1S( 2) + .01 
MXP1 = MX ♦ 1 
IB = 01S( 1) + .01 
IF( IB-l )78. 79,78 

78 CALL SE T( 01 , l.MZ) 

IJ * 4 - MZ 

DO 80 I = l.MZ 
IJ = IJ + MZ + 1 

80 Dl( 1+4) = DIS(IJ) 

GO TO 81 

79 CALL MAKE (Cl, 01S) 

81 CALL SP IT ( C 1 , 2PC1) 

AO * 1.0/ENN 

CALL ADO( AC , 01, 0., Dl, Dl) 

CALL AOD(AC, D1S • 0., D1S, D1S) 

PRINT 1002 
JKM = 1 

DO 5 NP = 1,60 

READ 111, ( IP(NP, I ) , I - 1»4)»(AP(NP»J), J = 1,3) 

PRINTlll,( IP(NP, I ) , I = 1,4),(AP(NP,J), J = 1,3) 

AP( NP ,4 ) = AP ( NP , 1 ) 

AP( NP,5) = 0. 

I F ( IP ( NP , 1) )6,5,5 
5 JKM = JKM + IP (NP , 1 ) 
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6 CONTINUE 

JKMM1 = JKM - 1 
NPM = NP - 1 
PRINT 1003, JKM, NPM 
TT = - HH 

PRINT 1004, MU, HZ 
00 48 I = 1 , NN 

READ 606, ( U IK , I ) , K = 1 ,MU) , I ZI K, I ) , K = 1,MZ) 

TT = TT + HH 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CONDITION INPUT DATA (START)...* 

DO 63 J = 1 ,MZ 

63 Z(J,I) = Z< J,I)*SZ(J+4) 

DO 64 J = 1 , MU 

64 U(J,I) = U( J ,1 )*St ( J+4) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CONDITION INPUT DATA (END) * 

48 PRINT1010,TT,(U(K, I) , K = 1 , MU) , ( Z ( K , I 1 , K * 1,MZ) 

Cl 1 1 1 1 1 II II II II 1 1 II 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 INPUT CATA (END) $ 

CLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL ITERATION LOOP (START) * 

00 12 LL= l » ITR 

CSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS SET MATRICES (START) * 

CALL SET ( U1 • MU, 1) 

CALL SET( X 1 , MX, 1) 

CALL SET( X J ,MX , 1 ) 

CALL SETIY, MZ, 1) 

CALL SET(EE,1,MZ) 

CALL SET(DZ1,1,MZ) 

CALL S£T(0ZN,1,MZ ) 

CALL SET1SUM, JKM, JKM) 

CALL SET (X J I , JKM, MZ) 

CSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS SET MATRICES (END) S 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC COMPUTE STATE TRANSITION MATRIX. 

CALL EAT( A , HH, PHI, APHI, DUM, DUD) 

JK = 0 

DO 10 NP = l.NPM 
IG = IP(NP,2J 
JJ = IP(NP,3) 



IF( 

IP 

(NP, 

n- 

1)71 

, 7 C, 

71 

70 

JK 

= 
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I I 
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( JK- 

1 ) * JKM 

♦ JK 
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SUM 

( I 
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1 . 

/AP ( 
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)**2 

71 

I F ( 

IP 

(NP, 

2) 

- 3) 

10,1 

6,16 

16 

IF ( 

IG 

-5 )8 

,10 

, 10 




8 J I = (JK-D+MZ + 4 ♦ I P ( NP, 3 ) 

XJI ( JI) = AP ( NP , 3 ) 

I F ( LL- 1)10, 55, 10 
55 IF(IP(NP, 21-3)10,17,20 
17 DZH4+JJ) = AP(NP»l)*AP(NP»3) 

GO TO 10 

20 DZN ( 4+ J J ) = AP(NP,1)*AP(NP,3) 

10 CONTINUE 

DO 51 I = 1 , MU 

51 Ul( 1+4) = <U(I,1) + U(I,2))*.5 
DO 24 J = 1, MX 

24 X1(J*4) = Z(J,1) - CZKJ+4) - DZN(J+4) 

C 111111111 111 1111111111 1111 ill 11111 COMPUTE INITIAL STATE ( START )... * 

CALL MULT ( PHI , XI, X2) 

CALL MULT ( 8 , Ul, CUM) 

CALL MULT (APHI, DUM, DUO) 

CALL ADD( 1.0, X2, l.0,DUD,X2) 

DO 47 I = 1 , MX 
47 Y ( I *4 ) = X2 ( I +4) 

TT = 0. 

Cl 1111 111 1111 111 mill 11 11 1111 111 11 COMPUTE INITIAL STATE (END) * 

CTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT TIME LOOP (START) * 
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DO 41 I = l.NNMl 
TT » TT + hH 
JK =0 

CGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG COMPUTE GRADIENT ( START 

DO 4 NP = l i NPM 
JK = JK + IP(NP'l) 

IG = I P ( NP » 2 ) 

IF( IG-4)82,4,82 

82 II = I P ( NP » 3 ) 

JJ = I P (NP >4) 

JI =* (JK-l)*MZ+4 
IF ( IP( NP, 11)25,83,25 

25 00 49 J = 1 , MX 

49 X J ( J+4 ) = XJKJI + J) 

CALL MULT ( PHI , XJ, DUM) 

DO 29 J = 1 , MX 
29 XJKJI + J) = OUM ( J *4 ) 

IF ( MZ-MX) 9,83,9 
9 CALL MULT ( F , XJ, DUM) 

DO 13 J = MXP1 , MZ 

13 X JI ( JI + J ) = DUMJ-MX+4) 

83 IJ = II + 4 - MX 
I F ( IG-1 )4 , 26 , 34 

26 IF(LL-1)33,33,32 

32 DO 42 J * 1, MX 
IJ = IJ + MX 

X J ( J+4 } = APHI( IJ ) *( XI ( JJ+4)+X2(JJ+4) )*.5*AP(NP,3) 

42 XJK J I + J ) > XJKJI+J) + XJ( J+4) 

GO TO 23 

33 DO 43 J = 1 , MX 
IJ = IJ + MX 

XJ< J+4) = APHKIJ )*(Z( JJ, I)*Z( JJ, 1*1) )*.5*AP(NP,3) 

43 XJKJI + J) = XJKJI + J) * X J ( J+4 ) 

34 IFUG-2) 23,35,23 

35 DO 45 J = 1 »MX 
IJ = IJ + MX 

X J ( J+4) = APHI(IJ)*(U( JJ,I)+U( JJ.I+l) )*.5*AP(NP,3) 

45 XJKJH-J) = XJKJI+J) + X J ( J+4 ) 

23 IF ( MZ-MX) 7 ,4,7 
7 I F ( I G— 3 )84,84,85 

84 CALL MULT (F,XJ»DUM) 

DO 14 J = MXPl.MZ 

14 XJKJI + J) = DUM ( J -MX+4 ) ♦ XJKJI + J) 

85 IF ( IG -5)4,15,18 

15 XJK JI + MX + II) = XJK JI+MX + II ) + X2(JJ+4)*AP(NP,3) 

GO TO 4 

18 I F ( IG-6 )4 , 2 8 ,4 

28 XJKJI+MX+II) = X J I( JI+MX+I I ) + U( JJ, I + 1)*AP(NP,3) 

4 CONTINUE 

DO 53 J = 1 , MU 

53 Ul( J+4) = L(J*I+l) 

IFIMZ-MX) 11,27,11 

11 CONTINUE 

CALL MULT ( F , X2 , OUM) 

CALL MULT ( G , Ul, CUD) 

DO 54 J = MXPl.MZ 

54 Y( J+4) = DUM( J+4-RX) ♦ 0UD(J+4-MX) 

27 JI = ( JKM— 1 ) *MZ + 4 
DO 52 J = 1 , MZ 

XJKJI + J) = ZU.I + l) - Y ( J+4) 

52 EEIJ+4) = E E ( J +4 ) ♦ D1 ( J+4) *X J I( J I+J ) **2 

CGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG COMPUTE GRADIENT (END) * 

IF(LL-ITR)68, 67,68 
67 I F ( 1-1)59,61,59 
61 PRINT 1004, MU, MZ 
59 00 69 J = 1 , MZ 
69 Y < J +4 ) = Y( J+4)/SZ( J+4) 

PRINT 1010, TT, (U( J, 1+1), J = i,MU) , (Y( J+4) , J = 1,MZ) 
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68 IF( I-NNMX J72,4X ,72 

72 CALL HAKE (XI, X2 ) 

DO 30 J = l, MU 

30 UXJ+4) = (U(J,I+X> ♦ U( J, 1+2) >*.5 

CSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS COMPUTE STATE (START) * 

CALL MULT ( PHI , XX, X2) 

CALL MULT ( 8 , UX, DUD) 

CALL MULT ( APHI , OLC, OUM ) 

CALL A00( X.0,X2,X.C,0UM,X2) 

DO 2 J = X , MX 
2 Y(J+4) = X2 ( J+4) 

CSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS COMPUTE STATE (END) * 

CEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE SET UP SYSTEM OF EQUATIONS $ 

4X CALL TR ICK ( XJ I , OXS, SUM) 

CTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT TIME LOOP (END) * 

LLM X = LL-X 

PRINT X005, SUM( JKM* JKM +4 ) , LLMX 
CALL SPIT(EE»4HEE*D) 

DO 60 I * I »MZ 

60 EE ( I +4 ) * SQRT(EE(I+4J/(DX(I+4)*ENN)>/SZ(I+4> 

CALL SPIT ( E E , 3FRMS) 

CODODDDDODDDODDDCDDDDDDDDDDDODODODD...MULT DIAS BY X.0005 TO CONDITION. $ 
CALL SET( DUM, I, JKM) 

DO X J = X , JKM 
JJ * JKM* J ♦ 4 
IJ * ( J- X ) * JKM ♦ J + 4 
DUM ( J+4 ) = SUM ( J J ) 

X SUM ( I J ) = SUM ( I J ) * X. 0005 
CALL SP IT ( CUM, 4HC-RA0) 

53X FORM AT ( 55 H MOULD YOU BELIEVE ) 

PRINT 53X 

CSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS SOLVE SYSTEM OF EQUATIONS t 

CALL ADO( F INK, SUM, 0., SUM, SUM) 

IF(LL-ITR)74, 73,74 

73 CALL M I NR ( SUM , JKM, JKM, XJI) 

CALL SET( ODD , JKMMX, X) 

JJ = 4 

DO 75 J = X, JKMMX 
JJ = JJ + JKM 
75 DUD( J+4 ) = SUM ( J J ) 

74 CALL SOLVE! SUM , CC) 

CALL SP IT ( CC , 2FDC) 

CUUUUUUUUUUUUUUUUUUUUUUOUUUUUUUUUUU UPDATE COEFFICIENTS (START) .. ..$ 

JK = 0 

DO 12 NP = X.NPM 
IG = IP(NP,2) 

II = IPINP.3) 

JJ = IP(NP,4) 

JK = JK + I P( NP , 1 ) 

IF( JK) X2, X2 , 50 
50 IF(IP(NP,X)-X)57,56,57 

56 AP ( NP ,4 ) = AP(NP,4) + DCJJK+4) 

57 IF(IG-X)X2* 19*21 

X9 IJ = ( I I-X )*MX + JJ + 4 

A ( I J ) = A ( I J ) + 0C< JK+4)*AP(NP,3> 

GO TO 12 

21 IF( I G- 21X2,22, 37 

22 IJ = (I I- X ) *MU ♦ J J + 4 

8 ( I J ) = B ( I J ) + 0C( JK+4)*AP(NP,3) 

37 I F < IG-3 >12,36,38 

36 021( 11+4) = QZXII+4) + DC( JK+4 ) * AP ( NP , 3 ) 


I) 


1 
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♦ 0C( JK+4)*AP(NP,3) 


38 XFC IG— 4)12,39,31 

39 OZN ( 1 1 +4 ) = D Z N ( 1 1 *4 ) 

31 I F ( IG-5)12,40.44 

40 I J = (I I- 1 ) *MX ♦ J J + 4 

F l I J J = F ( I J ) + OC ( JK+4) * AP ( NP, 3 ) 
44 IF(IG-6)12.46,12 
46 IJ = ( I I-1)*MU ♦ JJ + 4 

G( I J ) = G(IJ) ♦ DC( JK+4)*AP(NP,3) 

CUUUUUUUUUUUUUUULUUUUUUUUUUUUUUUUUU .... 


UPDATE COEFFICIENTS (END). 


12 CONTINUE 

CLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL... ..ITERATION LOOP (END) 
CALL SP IT ( A , 1HA ) 

CALL SPITIB, IH8 ) 

CALL SP IT ( DZ1 « 3H0Z1) 

CALL SPIT ( C ZN« 3HCZN) 

CALL SP IT ( F , 1HF ) 

CALL SP IT ( G * 1HG ) 

CALL MAKEIOUM, XJI) 

CALL INVRIXJT, SUM) 

CALL MULT ( SUM , DUM, XJI) 

CALL SP IT ( X J I « 4HS-1S) 

CALL MULTISUM, DUO, XJI) 

CALL SPIT (XJI, 5H OC* ) 


PRINT 1002 
JK = 0 

DO 58 NP = l.NPM 
IF( IP(NP,l)-l)58, 77,58 
77 JK = JK. «• 1 

I J = ( JK-1)*JKMM1 ♦ JK + 4 
APINP.5) = SQRT(FINK/SUM( I J) > 

58 PRINT 1 11 , ( I P( NP , J) , J = 1,4) , ( APtNP, J) , J = 1,5) 
STOP 
END 


$ 

$ 


CONCLUDING REMARKS 

A FORTRAN program has been developed which computes maximum-likelihood 
estimates or determines the parameters of a linear model which minimize the weighted 
mean-square difference between the model and the measured responses. The responses 
can be any linear transformation of the state and control variables. A priori information 
can be included and the initial state values and biases can also be determined. Several 
different parameter estimates are described and shown to be either identical or similar. 
An estimate of the error covariance matrix of the estimated parameters which is identical 
to the Cramer-Rao bound is derived. 

The program is described to enable the prospective user to apply the method to his 
particular problem after first running the check case provided. A few trouble-shooting 
tips are given as well. Copies of the program and test case will be provided by the author 
upon request. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 20, 1972. 
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APPENDIX A 


SUBROUTINES USED IN PROGRAM 

The subroutines described in this appendix are used in the main program NEWTON 
for a variety of matrix operations. In most cases the operation is simply the name of 
the subroutine. Since the matrices have vector format and all arrays of the subroutines 
appear in the call statement, all dimensioning of the arrays is done in the main program. 
The subroutines have a variety of uses besides the program discussed. The user is 
cautioned that little error checking is performed beyond checking the matrices' 
compatibility. 
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APPENDIX A - Continued 


Subroutine ADD(P,A,Q,B,C) 

Purpose: To add two matrices or vectors after each has been multiplied by a scalar 
Formula : C = PA + QB 
Flow chart: 



Description: Subroutine ADD can operate on two compatible matrices in the vector for- 
mat to (1) add, (2) subtract, (3) multiply by a scalar, or (4) perform some combination 
of these operations. 


Pro gram variables : A 

B 

C 

P 

Q 


first matrix or vector 

second matrix or vector 

weighted sum or difference of A and B 

first scalar 

second scalar 


Subroutine listing: SUBROUTINE ACDCP* A» Q« B» C) 

C VECTOR FORMAT 


DIMENSION A ( I ) f Bill, 
N = A ( 3 ) * 4.01 

DO 1 I 

= 5 * N 

1 cm = 

P* A C I) + Q*BC II 

cm = 

A ( I) 

C ( 2 ) = 

A( 2 I 

C( 3 ) = 

AC 3) 

C C4 I = 
RETURN 
END 

AC 4 I 
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APPENDIX A - Continued 


llllllllllllllllll■ll■ll!;llllllllll 


Subroutine EAT(A,T,PHI,APHI,A2,A3) 

Purpose: To form a state transition matrix and its integral for use in computing time 
histories 


Formula: 


6 

PHI = I 

n=0 


(AT) n 

nl 


e 


AT 


= y A n- ^T n _ r 

" tx n: J < 


APHI 


,AT 


dt 


Flow chart: 



Description : Subroutine EAT employs a finite power series to approximate the state 
transition matrix PHI and its integral APHI. The state is calculated by using the fol- 
lowing formula: 


X i+1 = (PHI)X i + (APHI) (B) 


u i+l + u i 


2 
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APPENDIX A — Continued 


The subroutine uses six terms of the power series but can be easily modified to include 
as many as is necessary. Since it is also assumed that the control input is constant over 
the time interval, a more accurate method of integration may be necessary in some cases. 


P rogram variables : A 

A2 

A3 

APHI 

BB 

G 

I 

II 
JJ 
PHI 
T 


stability matrix of x = Ax + Bu 
dummy matrix 
dummy matrix 


integral of state transition matrix, 



dt 


floating-point variable equal to term index 
scalar coefficient in power series 
general index 

number of rows of A and therefore of PHI and APHI 
number of columns of A and must be equal to II 
state transition matrix, e-^ T 

time increment over which transition matrix is desired. 
The value of T must be kept small compared with 
shortest time constant of system and cannot in gen- 
eral, be made the running time. 


Subroutine listing: subroutine eatca, t, phi, aphi, A 2 , A3> 

C VECTOR FORMAT 

DIMENSION All), PhUl), APHI(I), A2 Cl) , A3C1) 
123 FORM A T ( 55H MATRIX IS NOT SQUARE FOR EAT, STUPID 
II * A ( 1 ) * .01 
JJ * AC 2) + .01 
IFC 1 1 — J J ) 4, 3 » 4 

3 CALL SETCPHI, II, II) 

CALL SETCAPhI, II, II) 

DO 1 I = 1,11 

I J =» C 1-1 )* 1 1 ♦ I *4 

1 PHICIJ) = 1.0 
CALL MAKE ( A 2 , PHI ) 

G = 1.0 

DO 2 I = 1,6 
BB = I 
G = G*T/BB 

CALL AD0C1.0, APHI, G, A2, APHI) 

CALL MULT l A , A2, A3) 

CALL MAKE ( A2, A3) 

CALL ACDC 1 • , PHI, G, A2 , PHI) 

2 CONTINUE 
GO TO 5 

4 PRINT 123 

5 RETURN 
END 
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APPENDIX A - Continued 
Subroutine INVR(A,B) 


Purpose: To invert a matrix 


Formula: B = A - * 


Flow chart : f g tart 


Subtract multiples of Kth row to make 
off-diagonal elements of Kth column zero 

t 

Perform identical operations on 
a matrix that is initially I 

] 

Form product of diagonals for the 
determinant. Divide by diagonal 
elements making original matrix I 

C^s KIN? ). 5° 

| Yes 

Print the value of | A j 



Description: Subroutine INVR inverts the matrix A by using a Gauss-Jordan pivotal 
method in which the matrix A is diagonalized by subtracting multiples of succeeding 
rows to make off-diagonal elements zero. Identical manipulations are made on a 
matrix which is initially an identity matrix the size of A. The diagonal elements of the 
diagonalized matrix are then used as divisors for the corresponding rows of the other 
matrix and thus make the latter matrix A -1 . The determinant of A is also determined 
and printed. 

Program variables: A matrix (square) to be inverted 

B inverse of A 

DETA determinant of A 
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APPENDIX A — Continued 


I 


row index 



J 

column index 


K 

pivot index 


IJ,KPI,KK, ^ 
KIX,KL,KIL/ 

general indices 


AKK1,AKIK 

general variables 

Accuracy: 

[a - 1 A - l| < 10" 10 

1— J off- diagonal elements 


(N = 15) 


(Accuracy, of course, is problem dependent. The value given reflects the author's 
experience on the Control Data 6600 computer with single-precision arithmetic.) 


Subr outine listing : subroutine invria, ei 

C VECTOR FORMAT , OIASCNAL ELEMENTS MUST BE NON-ZERO. 
OIMENSION A(l). B(l) 

101 FORMAT ( El 5. 7 ) 

N - Mil ♦ .01 
CALL SETIB, N, NJ 
00 6 I =* l.N 
IJ = t I— 1 »*N *1+4 

6 BIIJI » 1.0 
OETA - 1.0 
00 3 K • l.N 
KP1 * K + 1 

KK * ( K— 1 1*N ♦ K ♦ 4 
AKK1 - 1. 0/ A( KK 1 
OETA = 0ETA4AIKK) 

DO S J > l.N 
IJ » ( K— 1 J *N + J •* 4 
5 B(IJ) * B«IJ)*AKK1 
DO 1 J = K,N 
IJ - ( K— 1 )*N ♦ J ♦ 4 

1 AUJ) « Ail J)*AKK1 
DO 2 KI = l.N 

IF( K - KI 17,2,7 

7 K IK » ( K I— 1 1 *N ♦ K ♦ 4 
AKIK > A(KIK) 

DO 4 L * l.N 
KL « ( K— 1 1* N ♦ L ♦ 4 
KIL >= (KI- 1 )*N ♦ L ♦ 4 
B(KIL) « B ( KIL I - AK I K*B ( KL I 
4 CONTINUE 
I F ( K— N 18,9,8 

8 DO 10 L 1 KPl.N 

KL *= ( K— 1 J*N ♦ L ♦ 4 
KIL » l KI — 1 I *N ♦ L + 4 
AiKILI » A(KIL) - AKIK*A(KL) 

10 CONTINUE 

9 CONTINUE 

2 CONTINUE 

3 CONTINUE 
PRINT 101, OETA 
RETURN 

END 
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APPENDIX A - Continued 
Subroutine LOAD(M,A,B,C,D) 


Purpose : To read input cards of M arrays (matrices or vectors) 


Flow chart: 




Description : Subroutine LOAD inputs up to four matrices through the repeated use of sub- 
routine LODE. All arrays are put into the vector format described previously. 


Program variables : A 


B 


C 


D 


M 


Subroutine lis ti ng : 

DIMENSION All). B(l), Clli. 
CALL LODE ( A ) 

IF ( M-I ) 8 * 8 * 2 
2 CALL LODE(B) 

I F ( M-2 ) 8* 8 * A 
4 CALL LODEIC) 

IFIM-318, 8,6 
6 CALL LODE ( C ) 

8 RETURN 
END 


first array 
second array 
third array 
fourth array 
number of arrays 

SUBROUTINE LQAC(M, A, B, C, D) 
C VECTOR FORMAT 
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APPENDIX A - Continued 
Subroutine LODE (A) 

Purpose: To read input cards of a single array (matrix or vector) 


Flow chart: 


(start) 


Determine n,JJ 


1 


Read cards 



1 

( Last row? 

\ No 



Description: Subroutine LODE reads the first card to determine the number of rows 
and columns of the array. The required number of cards are read for each row until 
all rows are read. The array is put into the vector format previously described. 


Prog ram variables : A 

I 

II 
J 
JJ 

L 


input array 
row index 
number of rows 
column index 
number of columns 

general index 
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APPENDIX A - Continued 


Subroutine listing : SUBROUTINE LOOECA) 

C VECTOR FORMAT 

DIMENSION All) 

700 FORMAT I7F10.3) 

444 FORMAT (7X*I3*7X*I3) 

READ 444* II* JJ 
A ( 1 ) = II 
A (2) = J J 
A ( 3 ) * II*JJ 
A (4 ) = 0. 

IF( At 3) )2 *3*2 

2 CONTINUE 

L = 4 - JJ 
DO 1 I = 1,11 
L = L ♦ JJ 

1 READ 700* ( A ( L+ J ) , J = 1*JJ) 

3 CONTINUE 
RETURN 
END 
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APPENDIX A - Continued 


Subroutine MAKE(A,B) 
Purpose : To make array A equal to array B 

Formula: A = B 


Flow chart: 


Calculate N 




Descri ption: Subroutine MAKE sets all the elements of array A equal to those of array B, 
including the elements which indicate the size of the array. The vector format is used 
as previously described. 


Program var iables : A 

B 

I 

N 


new array 
original array 
general index 

number of elements including those giving number of rows, 
columns, and so forth 


Subroutine listing: SUBROUTINE MAKE(A, E) 

C VECTOR FORMAT 

DIMENSION All), B (1) 

N = B(1)*B(2) ♦ 4.01 
DO 1 I = ltN 
1 A ( I ) = Bill 
RETURN 
END 
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APPENDIX A - Continued 


Subroutine MINR(A,I,J,B) 
Purpose: To form the minor of matrix A 


Formula: 


B = MINOR [a] 

I,J 


Flow chart: 



Description: The minor of matrix A is formed by eliminating the Ith row and Jth column. 


Program varia bles : A 

B 

I 

J 

K,L,H,JJ 


matrix for which minor is sought 
resulting minor matrix 
row to be eliminated 
column to be eliminated 
general indices 
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APPENDIX A — Continued 


Subroutine listing: 


SUBROUTINE MINR ( A » I, J, 
C VECTOR FORMAT 

DIMENSION A ( 1 ) t B(l) 

M = AUi + .01 
N = A( 2) + .01 
IF( I >11,10.11 

10 I = 99999 
MM1 ■= M 
GO TO 12 

11 MM1 = M - 1 

12 IF( J ) 13 * 1 A * I 3 
14 J = 99999 

NM1 = N 
GO TO 9 

13 NM1 = N - 1 

9 DO 1 K = 1,M 
DO 1 L = 1,N 
I F ( I-K ) 2 « 8 . 3 

2 II « K - 1 
GO TO 4 

3 I I = K 

4 I F ( J— L ) 5» 8 » 6 

5 JJ = L- 1 
GO TO 7 

6 JJ = L 

7 B ( I I*NM1 + JJ - NR1 + 4) 

8 CONTINUE 
l CONTINUE 

B(l) = MM 1 
B ( 2 I = NM1 
B ( 3) = MM 1*NM1 
B (4 ) = 0. 

RETURN 

END 


A(K*N * L - N * 4) 



APPENDIX A — Continued 


Subroutine MULT(A,B,C) 

Purpose: To compute the product of two arrays (matrices or vectors) 

Formula: C = AB 
KK 

C ij = Z a ik b kj 
k=l 



Description: Subroutine MULT generates an array C which is the product of arrays A 
and B. Since the elements of C are changed during the operation, array C cannot be 
the same as either array A or B. 


Program variables: A 


first array in product 


B second array in product 


C product array 
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APPENDIX A — Continued 


I row index 

II number of rows in A 

J column index 

JJ number of columns in B 

KK number of columns in A 

LL number of rows in B 

IJ,IK,KJ general indices 


Subroutine listing: 


SUBROUTINE MULT ( A , B, 
VECTOR FORMAT - NEITHER 



DIMENSION 

A 

(1 

) , 

a 

(1) 


II = 

All) 

■f 

• 

01 




JJ = 

B ( 2 ) 

♦ 

• 

01 




KK = 

A ( 2 ) 

+ 

• 

01 




LL = 

B ( 1) 

+ 

• 

01 




CALL 

SET ( 

C * 

I 

I , 

j 

J) 


I F ( KK 

. - LL) 

12 

,11 

, 

12 

12 

PRINT 

123 






123 

FORMAT ( 55 

H 

MATRI 

C 

ES 

11 

00 l 

I = 

It 

II 





00 l 

J = 

1, 

J J 





OQ 1 

K a 

1, 

KK 





I J = 

( I-l 

)* 

J J 


j 

+ ■ 


IK = 

( I-l 

)* 

KK 

+ 

K 



KJ = 

(K-l 

)* 

J J 

+ 

J 

+ 

1 

C( I J J 

= C 

( I 

J) 

♦ 

A 

( IK 


cm 

= A ( 

1) 






C ( 2 ) 

= B ( 

2) 






C ( 3 ) 

= II 

*J 

J 





C(4) 

= 0. 







C) 

A NOR 
C( 1> 


B CAN BO EQUIVALENT TO C 


INCOMPATIBLE FOR MULT, CUM-DJM. 


RETURN 

ENO 



APPENDIX A - Continued 


Subroutine SET(A,II,JJ) 

Purpose: To generate an array with all its elements set to zero 
Formula : A = [oj 
Flow chart: 


Set first elements of A 



Description: Subroutine SET generates a null array of desired size. The array must 
still be dimensioned although it originates by calling the subroutine. A vector format 
is used as was described earlier. 

Program variables : A null array generated 

II number of rows 


JJ number of columns 


N total number of elements including special elements 


SUBROUTINE SET ( A t II» JJ) 
C VECTOR FORMAT 

DIMENSION A ( 1 1 
A ( 1 ) = II 
A (2) = JJ 
A ( 3 ) = II*JJ 
A ( 4 ) = 0. 

N = II*JJ ♦ 4 
DO I I = 5.N 
1 A ( I ) = 0.0 
RETURN 
END 








APPENDIX A — Continued 


Subroutine SOLVE(A,C) 

Purpose: To solve a linear system of algebraic equations Y = 


Formula: C = 


where 



Flow chart: 



Subtract multiples of upper rows 
to triangularize A 


I 


Last row? 


No 


( Last column? No 
1 


Solve for next element of C 
starting with last 


1 


Multiply data by diagonal element 


First row? 


Print determinant value 


Set first elements of C 








APPENDIX A — Continued 


Description : Subroutine SOLVE employs Gauss' pivotal or elimination method to first 
triangularize the matrix B and then solve for the elements of the solution vector C in 
reverse order. The product of the diagonal terms of the triangularized matrix is the 
determinant, the value of which is printed to indicate a proper range of amplitudes of 
numbers involved. The matrix B consists of A except for the right-hand column which 
is Y. Both matrix A and vector C use the vector format described earlier. 


Program variables : A 

C 

DETA 

N 


juxtaposition of matrix of coefficients of simultaneous 
equations B and constant vector Y 

solution vector 

determinant 

number of unknowns, simultaneous equations 


Subroutine listing: 


SUBROUTINE SOLVEIA, C) 

C VECTOR FORMAT 

DIMENSION All), C< I) 

101 FORMAT! E15. 7) 

N = All) - .99 
NM1 = N - 1 

NP1 = N + 1 

NP3 = N - 3 

00 6 K = 1 ,NM1 

AKK 1 = 1 . / A ( K*NF 1 ♦ K - NP3 ) 

K1 = K + 1 
00 6 M = K1,N 

AKIK = A(KI*NP1 + K - NP3)*AKK1 
DO 6 L = Ki.NPl 
KIL = M*NPl + L - NP3 
6 A(KIL) = AIKIL) - AK I K*A I K*NP 1 + L - NP3 ) 

OETA = A! 5) 

C ( N + 4 ) = A!NP1*NP1 - NP3 ) /A ( N*NP 1 «- 3) 

00 16 I = 2 » N 

Nil = N - I *■ 1 

C ( NI 1 + 4) = AIM 1*NP1 + 4) 

N 1 2 = Nil + 1 
DO 8 J = N 1 2 , N 

8 CINIH-4) = CINIH-4) - C ( J+4 ) * A I N 1 1*NP1 + J - NP3) 
CINI1+4) = CINI1 + 4 )/A(NU*NPl + Nil - NP3 ) 

16 OETA = DETA*A(I*NP1 + I - NP3 ) 

PRINT 101, OETA 

cm = 1.0 

C ( 2 ) = All) - 1.0 
C < 3 ) = C ( 2 ) 

C ( 4 ) = 0.0 

RETURN 

END 
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APPENDIX A — Continued 


Subroutine SPIT(A,B) 


Purpose : To print array 
Flow chart: 


A 



Description : Subroutine SPIT prints a heading which includes the number of rows and 
columns of array A and then prints the array row after row in a straightforward man- 
ner. A vector format is used that was described earlier 


Program variables: 

A 

array to be printed 


n 

number of rows 


jj 

number of columns 

Subroutine listing: 




SUBROUTINE SPITIA, B> 

C VECTOR FORMAT 

DIMENSION A(l) 

505 FORMAT! 5E12. 4) 

321 FORMAT ( 8H MATRIX ,A5,4HHAS , 14, 10H ROWS AND ,I4,8H COLUMNS) 
II = A(l) ♦ .01 
JJ = A<2 ) + .01 
PRINT 321, B, II, JJ 
L = 4 - JJ 
IF C 11)3,3,2 

2 DO l I = 1,11 
L = L + JJ 

1 PRINT 505 , ( A(L+J ) , J = 1,JJ> 

3 RETURN 
END 
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APPENDIX A — Continued 
Subroutine TRICK(A,B,C) 

Purpose: To form a matrix of coefficients of a system of simultaneous linear algebraic 
equations 

Formula : C n+ ^ = C n + A^B’A 
where 
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APPENDIX A - Concluded 


Description : Subroutine TRICK is designed for use in evaluating a matrix which is 

involved typically when a system of linear simultaneous algebraic equations is formed. 
A vector format which was described earlier is used for all the arrays. 


Prog ram variables : A 

B 

C 

II 

JJ 

I,J,IJ 


matrix to be "squared" 
weighting matrix 
simulation matrix 
number of rows in A and C 
number of columns in A and B 
general indices 


Subroutine listing: 


SUBROUTINE TRICKIA, B, C) 

C VECTOR FORMAT 

DIMENSION AID, B(l) , C C 1 * 

IB = BID -* .01 
J B = B(21 + .01 
II = A< II ♦ .01 
JJ = At 2) ♦ .01 
IFt ( I B— 1 ) * ( JB-l) 13,4,3 
4 CONTINUE 

DO 1 J = 1, II 
DO 1 I = J,II 
IJ = (1-11*11 + J + 4 
DO 2 K = 1,JJ 

2 C ( I J I = C < I J 1 * At I*J J*K-JJ+41*At J*JJ+K-JJ+4)*B(K+4) 
1 C< J*II ♦ I - II * 41 = C < 1 J 1 

RETURN 

3 CONTINUE 

DO 11 J = 1,11 
DO 11 I = J,II 
I J = t 1-1 )* I I + J + 4 
DO 12 K - 1,JJ 
IK = ( 1-1 )* J J ♦ K + 4 
DO 12 L = 1,JJ 
KL = ( K— 1 1 * J J + L + 4 
JL = <J-1)*JJ + L ♦ 4 
12 CtlJ) = CtlJ) *■ A t IK )*BtKL)*A( JL 1 
11 C( J*II + I - I I ♦ 4) = CtlJ) 

RETURN 

END 
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APPENDIX B 


INPUT AND OUTPUT LISTINGS 


The input and output listings of this appendix are provided as an aid to the user. 
The user should first compute the example discussed in the text; then the input can be 
changed to suit the particular needs of the user. 


Input Listing 


200 10 
.05 .01 

A 4 4 


-1.0 .0789 

1 .0 

R 4 3 


D 1 

1 

7 





1 .0 

1 .0 

1 .0 

1 .0 

1 .0 

1 .0 

1 .0 

F 

3 

4 





0. 

n. 

0. 

0. 




0. 

0. 

0. 

0. 





-1.0 


.0789 




G 

3 

3 





O. 

o. 

0. 





0. 

n. 

0. 





0. 

n. 

0. 





BZ 

1 

7 





1 .0 

1 .0 

1 .0 

1 .0 

1 .0 

I .0 

1 .0 

su 

1 

3 





1 . o 

1.0 

1 .0 





LP 

1 

1 

1 

1 

99999. 

1 .0 

LP 


5 

1 

1 


1 .0 

LR 

1 

1 

1 

2 

99999. 

1 .0 

LR 


5 

1 

2 


1 .0 

LB 

1 

1 

1 

3 

99999. 

1 .0 

LB 


5 

1 

3 


1 .0 

LDA 

1 

2 

1 

1 

99999. 

1 .0 

LDA 


6 

1 

1 


1 .0 

LDR 

1 

2 

1 

2 

99999. 

1 .0 

LDR 


6 

1 

2 


1 .0 

L0 

1 

2 

1 

3 

99999. 

1 .0 

LO 


6 

1 

3 


I .0 

NP 

1 

1 

2 

1 

99999. 

1 .0 

NP 


5 

2 

1 


1 .0 

NR 

1 

1 

2 

2 

99999. 

1 .0 

NR 


5 

2 

2 


1 .0 

NB 

1 

1 

2 

3 

99999. 

1 .0 

NB 


=; 

2 

3 


1 .0 
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APPENDIX B — Continued 


NOA 

1 


2 

NOA 



6 

NDR 

1 


2 

NOR 



6 

NO 

1 


2 

NO 



6 

YP 

1 


1 

YP 



5 

YB 

1 


1 

YB 



5 

YO 

1 


2 

YO 



6 

END 

-1 



0 


0 


2 

2944E-06 

1 

6415E-03 

0 


0 


8 

6888E-06 

1 

3345E-03 

0 


0 


1 

841 6E-05 

1 

01 18E-03 

0 


0 


3 

0670E-05 

6 

7648E-04 

0 


0 


4 

461 4E-05 

3 

31 53E-04 

0 


0 


5 

9385E-05 

-1 

9951 E-05 

0 


0 


7 

41 05E-05 

-3 

7489E-04 

0 


0 


8 

7888E-05 

-7 

3023E-04 

0 


0 


9 

9846E-H5 

-1 

0829E-03 

0 


0 


1 

0910E-O4 

-1 

4301 E-03 

0 


0 


1 

1 478E-04 

-1 

7688E-03 

0 


0 


1 

1 604E-O4 

-2 

0963E-03 

0 


0 


1 

1 206E-04 

-2 

4100E-03 

0 


0 


1 

0206E-04 

-2 

7074E-03 

0 


0 


8 

5290E-05 

-2 

9861 E-03 

0 


0 


6 

1 064E-O5 

-3 

2442E-03 

0 


5 

OOOOE-O 1 

9 

9991 E-04 

2 

0536E+00 

0 


1 

OOOOE+OO 

7 

0519E-03 

4 

0 1 94E+00 

4 

9000E+00 

1 

OOOOE+OO 

4 

0948E-02 

4 

3730E+0 1 

9 

8000E+00 

1 

OOOOE+OO 

1 

8395E-01 

8 

2525E+0 1 

9 

5500E+O0 

1 

OOOOE+OO 

5 

1 324E-01 

7 

8187E+01 

9 

3000E+O0 

1 

OOOOE+OO 

1 

0379E+00 

7 

3338E+0 1 

2 

5000E+00 

1 

OOOOE+OO 

1 

7207E+O0 

1 

4618E+01 

-4 

3000E+ 00 

1 

OOOOE+OO 


2 1 

2 1 

2 2 

2 2 

2 3 

2 3 

3 I 

3 1 

3 3 

3 3 

3 3 

3 3 

1 • OOOOE+OO 8.9372E-05 

-8.0941E-04 4.7933E-04 

1 • OOOOE+OO 1 .6384E-04 

— 7*471 IE — 04 5.1709E-04 

1 • 0000E+00 2.2256E-04 

-6.8167E-04 5.5005E-04 

1 • OOOOE+OO 2*6481 E— 04 

-6.1363E-04 5.7803E-04 

1 « 0000E + 00 2.9005E-04 

-8.4357E-04 6 • 0087E — 04 

1 • 0000E+00 2.9786E-04 

-4.7208E-04 6.1847E-04 

1 • OOOOE+OO 2.8800E-04 

-3.9974E-04 6.3076E-04 

1 • OOOOE+OO 2.6036E-04 

— 3* 271 4E — 04 6.3773E-04 

l.OOOOE+OO 2.1502E-04 

— 2 • 54 84 E— 04 6.3940E-04 

l.OOOOE+OO 1.5216E-04 

— 1 • 8342E — 04 6 .358 1 E — 04 

l.OOOOE+OO 7.2150E-05 

— 1 . 1342E— 04 6.2708E-04 

l.OOOOE+OO -2.4530E-05 

— 4 . 5389E— 05 6.1334E-04 

l.OOOOE+OO -1.3725E-04 
2-0174E-05 5.9477E— 04 

l.OOOOE+OO -2.6526E-04 
8. 2782E— 05 5.7158E-04 

l.OOOOE+OO -4.0768E— 04 
1.4198E-04 5 . 440 1 E— 04 

l.OOOOE + OO — 5.6353E — 04 
1.9735E-04 5.12 34 E — 04 

l.OOOOE+OO 5 . 092 1 E— 02 

— 1 • 8679E + 00 5.1327E-02 

l.OOOOE+OO 2.0320E-01 

— 3.7056E+00 2.0201 E— 0 1 

l.OOOOE+OO 1 . 3997E+00 

— 9. 181 2E + 00 6.2049E — 01 

l.OOOOE+OO 4 . 5598E+00 

— 1 . 4529E + 0 1 1.4713E+00 

l.OOOOE+OO 8 . 5798E + 00 

— 1 • 3875E + 0 1 2.5101E + 00 

l.OOOOE+OO 1.2370E+01 

-1.3118E+01 3 . 4939E+00 

l.OOOOE+OO 1 .4568E+01 

— 4 . 8520E+00 4.11 88E+00 

l.OOOOE+OO 1 . 3844E+0 1 


99999. 1.0 

1 .0 

99999. 1.0 

1 .0 

99999. 1 .0 

1 .0 

99999. 1.0 

1 .0 

99999. 1.0 

1 .0 

99999. 1 .0 

1 .0 


-4 

1 952 E— 05 

2.2927E-05 

5 

OOOOE-02 


-8 

0879E-05 

4 .7857E-05 

1 

0000E-01 


-1 

1661 E-04 

7.4556E-05 

1 

5000E-01 


-1 

4900E-04 

1 . 0278E-04 

2 

OOOOE-O 1 


- 1 

7794E-04 

1 . 3227E-04 

2 

5000E-0 1 


-2 

0334E-04 

1 . 6278E— 04 

3 

00 00 E— 0 1 


-2 

251 3E-04 

1 . 94 03E —04 

3 

5000E-0 1 


-2 

4331 E-04 

2.2577E-04 

4 

OOOOE-Ol 


-2 

5785E-04 

2.5772E-04 

4 

50 00 E— 0 1 


-2 

6880E-04 

2.B962E-04 

5 

OOOOE-Ol 


-2 

7622E-04 

3.21 2 IE-04 

5 

5000E-01 


-2 

80 1 BE-04 

3.5224E-04 

6 

OOOOE-Ol 


-2 

B080E-04 

3 . B24 7E— 04 

6 

5000E-01 


-2 

7821 E-04 

4. 1 164E-04 

7 

OOOOE-Ol 


-2 

7258E-04 

4 .3955E-04 

7 

5000E-01 


-2 

640BE-04 

4 .6598E-04 

8 

OOOOE-Ol 


-4 

7040E-02 

1 .4466E-03 

8 

5000E-01 


-1 

8650E-01 

7.4774E— 03 

9 

OOOOE-O 1 


-5 

0905E-01 

2 • 6685E— 02 

9 

5000E-01 


-1 

1 024E+00 

7.7664E-02 

1 

OOOOE+OO 


-1 

8129E+00 

1 .7739E-01 

1 

0500E+00 


-2 

4881 E+00 

3.2771E-01 

1 

1 OOOE+OO 


-2 

9374E+00 

5.201 3E-01 

1 

1 500E+00 


-2 

9745E+00 

7.2736E— 01 


51 


APPENDIX B — Continued 


2 . 4401 e+oo 
-4.3333E+00 
3.0772E+00 
-4.3667E+00 
3.601 1E+O0 
-4.4000E+00 
4.0076E+00 
-4.4333E+00 
4.2933E+00 
-4.4667E+O0 
4 .4555E+O0 
-4.5000E+00 
4.4922E+00 
-2.2000E+00 
4.41 OBE+nO 

1 . OOOOE— O 1 
4 .2494E+nO 
1 .OOOOE— 0 1 
4.0453E+00 
1 .OOOOE-n 1 
3 .8051 E+nO 
1 . OOOOE— O 1 
3.5264E+O0 
1 • OOOOE— 0 1 
3.2069E+00 
1 .OOOOE— O 1 
2.8450E+00 
1 • OOOOE— 0 1 
2.4390E+O0 
1 . OOOOE — O 1 
1 .9880E+O0 
1 . OOOOE-O 1 
1 .491 5E+ O0 
1 • OOOOE — n 1 
9.4952E-nl 
1 . OOOOE— 0 1 
3.6255E-01 
1 . OOOOE-0 1 
-2.6828E-01 
1 . OOOOE— 0 1 
—9.41 28E— 0 1 
1 .OOOOE— O 1 

- 1 • 6538E + 00 
1 .OOOOE-0 1 

-2.401 3E+00 
1 .0000E-O 1 
-3.1 785E+O0 
1 • OOOOE — 0 1 

— 3 • 9795E+ 0 0 
1 . OOOOE— 0 1 

— 4 • 7980 E+00 
1 . OOOOE— O 1 
— 5.6275E+00 
1 .OOOOE— 0 1 
— 6 • 46 1 2E+00 
1 • OOOOE— 0 1 
— 7.2918E+O0 
1 . OOOOE— 0 1 
-8.1 123E+O0 
1 • OOOOE— O 1 
-8.91 52E+00 
1 . OOOOE— 0 1 


— 4.3445E+01 
1 .OOOOE+OO 
-4.5324E+01 
1 .OOOOE+OO 
— 4 . 6956E+0 1 
1 .OOOOE + OO 
—4 • 8328E+0 1 
1 .OOOOE+OO 
-4.9431 E+01 
1 .OOOOE+OO 
-5.0255E+0 1 
1 .OOOOE+OO 
— 5 • 0796E+0 1 
1 .OOOOE + OO 
-3 . 2028E+0 1 
1 .OOOOE+OO 
-1 . 3358E+0 1 
8.0833E-0 1 
-1 .4440E + 0 1 
6.1667E-01 
-1 . 5428E + 0 1 
4.2500E-0 1 
-1 .6301 E+0 1 
2.3333E-0 1 
-1 . 7035E+0 1 
4. 1 667E-02 
-1 .76 1 2E+0 1 
-1 .5000E-01 
-1 .80 1 2E+0 1 
-3.4167E— 01 
-1 . 8220E+0 1 
-5.3333E-0 1 
-1 .82238 + 0 1 
-7.28008-01 
-1 . 8008E+0 1 
-9. 1667E-01 
-1 . 7566E+0 1 
-1.1 083E+00 
-1 .68898+0 1 
-1 .30008+00 
-1 . 5973E+0 1 
-1 . 3000E+00 
-1 .40278+0 1 
-1 .3000E+00 
-1 . 1 873E+0 1 
-1 .30008+00 
-9.53448+00 
-1 • 3000E+00 
— 7.0347E+00 
-1 . 3000E+00 
— 4.3986E+00 
-1 . 3000E+00 
-1 .651 3E+00 
-1 .3000E+00 
1 . 1814E+00 
-1 .30008+00 
4 . 0736E+0 0 
-1 • 300 OE + O 0 
6.99928+00 
-1 . 3000E+0 0 
9.9323E+00 
-1 .3000E+00 


3 . 3565E+ 00 
1 .OOOOE+OO 
3 • 7779E+00 
1 .OOOOE+OO 
4. 1 528E+00 
1 .OOOOE+OO 
4. 47868+00 
1 .OOOOE+OO 
4.75338+00 
1 .00008+00 
4 • 975 1 E+00 
1 .00008+00 
5. 1 427E+00 
1 .OOOOE+OO 
2.61 41 E+00 
1 . 00008 + 00 
8. 1 551 E-02 
1 . OOOOE+OO 
8.86558-01 
1 .OOOOE+OO 

1 . 6687E + 00 
1 .00008+00 

2 • 4236E + CO 
1 .00008+00 
3. 1 470E+00 
1 .00008+00 
3.8350E+00 
1 .OOOOE+OO 
4 . 4843E+00 
1 .OOOOE+OO 

09 1 6E + 00 
1 .00008+00 
5. 65428+00 
1 .00008+00 
6. 1699E+00 
1 .OOOOE+OO 
6 * 6366E+00 
1 .OOOOE+OO 
7.05308+00 
1 .00008+00 
7.41 80E+00 
1 .OOOOE + OO 
7.0 1 49E+00 
1 • OOOOE+OO 
6. 57 1 2E + 00 
1 .00008 + 00 
6 * 09 1 28+00 
1 .OOOOE+OO 
=:.5792E+00 
1 .OOOOE+OO 
5.0399E+00 
1 . OOOOE+OO 
4 » 4 780E + 00 
1 .OOOOE+OO 
3 • 8983E + 00 
1 . OOOOE+OO 
3 • 3058E + 00 
1 .OOOOE+OO 
2* 7053E+00 
1 • OOOOE+OO 
2.101 9E+00 
1 .OOOOE+OO 


4 .0874E+00 
1 * 1 624E+01 
3.71 34E+00 
9.31 61 E+00 
3 . 3093E+00 
6 • 9329E+00 
2 . 8788E + 0C 
4 . 4878E+ 00 
2 . 4256E+ 0 0 
1 .99458+00 
1 .95358+00 

— 5.3294 E— 0 1 
1 . 4 665E+ 00 

-2.6036E+00 
1 . 0744E+00 
-3.73798+00 
8.8490E-01 
-4 .4332E+00 
7.73708-01 
-5.18038+00 
6. 1 776E-01 
— 5.9740E+00 
4.1926E-01 
— 6.8080E+0C 
1 .8068E-01 
— 7.6749E+0C 
— 9.5236E— 02 

— 8.5662E+00 

— 4 . 0548E— 0 1 
— 9.4728E+00 
-7.4683E-01 
-1 .03858+01 
-1 .1 159E + 00 
-1 . 1 291 E+0 1 
-1 . 5090E+00 
-1 . 2 I 82E + 0 1 
-1 .92268+00 
-1 .30448 + 01 
-2.3527E+00 

— 1 .38678 + 0 1 
— 2 • 7955E + 00 
-1 .461 8E+ 0 1 
— 3.2276E+00 
-1 . 5266E+ 0 1 
-3 . 6262E+00 
-1 .5802E+01 
-3 • 9887E+0C 
-1 .621 7E+01 
—4 .31 28E+00 
-1 . 6503E+0 1 

— 4 .5965 E + 00 
-1 • 6655E+0 1 
— 4.8382E+00 
-1 .6667E + 01 
-5.0366E+00 
-1 . 6 535E + 0 1 
-5. I 908E+00 
-1 . 6259E+ 0 1 
— 5.3002E+00 
-1 . 5835E+0 1 
— 5.3646E+00 
-1 . 5266E+0 1 


i . 2000E+00 
—2 . 7959E+00 
1 . 2500 E+00 
-2 • 5974 E+00 
1 . 3000E+00 
-2.381 5E+00 
1 . 3500E+00 
-2. 1504E+00 
1 . 4000E+00 
-1 .9070E+0C 
1 . 4500E+00 
-1 . 6538E+00 
1 . 5000E+00 
-1 . 4598E + 00 
1 . 5500E+00 
-1 . 3924 E + 00 
1 . 6000E+00 
-1 .3681 E+00 
1 . 6500 E+00 
-1 . 3042E + 00 
1 • 7000E+00 
-1.201 7E+00 
1 • 7500E+00 
-1 . 0624 E+00 
1 . 8000E+00 
—8 . 8766E— 0 1 
1 . 8500E + 00 
—6 . 795 1 E— 0 1 
1 . 9000 E+ 00 
—4 . 3995E— 0 1 
1 . 9500E+00 
-1.71 12E-01 
2. OOOOE+OO 
1 .2468E— 01 
2.0500E+00 
4 .4504E— 01 
2 . 1 OOOE+OC 
7 . 8748E— 0 1 
2 . 1 500E+00 
1 . 1 495E+00 
2 . 2000E+00 
1.51 05E+00 
2 . 2500E+00 
1 . 8503E+00 
2 . 3000E+00 
2.1 670E+00 
2 . 3500E+0C 
2 . 4589E+00 
2 • 4000E+00 
2 • 7245E+00 
2 • 4500E+00 
2 . 9625E+00 
2 • 5000E+00 
3. 1 720E+00 
2 . 5500E+00 
3.3521 E+00 
2.6000E+00 
3 . 5024 E+00 
2 .6500E+00 
3 . 6226E+00 
2 . 7000E+00 
3.71 26E+00 


9 . 2251 E— 0 1 
1 . 0982E+00 
1 • 2530E+00 
1 . 3857E +00 
1 .4952E + 00 
1 • 5808E+00 
1 .6437E + 00 
1 .6920E+00 
1 . 7337E + 00 
1 . 7686E + 00 
1 .79468+00 
1 .80988+00 
1 .81 20E+00 
1 • 7996E + 00 
1 . 7709E+00 
1 .7244E+00 
1 . 6588E+00 
1 . 5730E+00 
1 • 4662E + 00 
1 . 3375E + 00 
1 . 1 868E+00 
1 .01 53E+00 
8.2474E-01 
6. 1 704E-01 
3.941 3E — 0 1 
1 . 5809E — 0 1 
■8.8965E-02 
• 3.44B4E— 01 
6. 0730E-01 
8.741 IE-01 
1 . 1 430E+00 


52 



APPENDIX B - Continued 


-9.6934E+00 
1 • 0000E— 0 1 

- 1 .0439E+01 
1 .OOOOE-Ol 

-1.1 146E+01 
1 .OOOOE-Ol 
-1.1 B07E+01 
1 .OOOOE-Ol 
-1 . 24 1 4E+n 1 
1 .OOOOE-Ol 
-1 .2962E+01 
1 .OOOOE-Ol 

- 1 . 3444E+0 1 
1 • OOOOE— O 1 

-1 . 3855E+0 1 
1 .OOOOE-Ol 
-1 .41 SSE + O 1 
1 .OOOOE-Ol 

- 1 .4440E+01 
1 .OOOOE-Ol 

-1 .4606E+01 
1 .OOOOE-Ol 

- 1 • 4682E+0 1 
1 .OOOOE-Ol 

- 1 • 4665E+0 1 
1 .OOOOE-Ol 

-1 .4552E+01 
1 . OOOOE— O 1 
-1 .4341E+01 
1 .OOOOE-Ol 

- 1 .4031 E + 0 1 
1 • OOOOE-0 1 

- 1 .3621 E+Ol 
1 .OOOOE-Ol 

-1.3111 E+0 1 
1 .OOOOE-Ol 

- 1 .2501 E+Ol 
1 .OOOOE-Ol 

-1 • 1 792E+0 1 
1 « OOOOE — O 1 
-1 • 0986E+O 1 
1 . OOOOE — o 1 

- 1 . 0086E+O 1 
1 .OOOOE— 0 1 

— 9.0935E+00 
1 . OOOOE— O 1 
-8.01 30E+00 
1 .OOOOE-nl 
-6.8481 E+00 
1 • OOOOE— O 1 
-5.6033E+00 
1 .OOOOE-Ol 
— 4.2833E+nO 
1 .OOOOE-nl 
— 2.8936E+00 
1 .OOOOE-nl 

- 1 « 4397E+00 
1 .OOOOE-Ol 
7.2358E— 02 
1 .OOOOE-Ol 


1 . 2847E+0 1 
-1 .3000E+00 
1 .5719E+01 
-1 • 3000E+00 
1 . 8524E+0 1 
-1 . 3000E+00 
2.1 237E+0 1 
-1 . 3000E+0 0 
2 . 3837E+0 1 
-1 . 3000E+00 
2 . 6303E+0 1 
-1 .3000E+00 
2.8615E+01 
-1 • 3000E+0 0 
3 . 0754E+0 1 
-1 . 3000E+00 
3 . 2705E+0 1 
-1 . 3000E+0 O 
3 .4451 E+0 1 
-1 . 3000E+00 
3 . 598 1 E+0 1 
-1 . 3000E+00 
3 . 7282E+0 1 
-1 .3000E + 00 
3.8346E+01 
-1 . 3000E+0 0 
3.91 65E+0 1 
-1 . 3000E + 00 
3 . 9735E+0 1 
-1 • 3000E+0 0 
4.0052E+01 
-1 . 3000E+00 
4.01 15E+01 
-1 • 3000E+0 0 
3 . 9925E + 0 1 
-1 . 3000E+00 
3.9485E+01 
-1 . 3000E + 00 
3.8801 E+Ol 
-1 . 3000E + 00 
3.7878E+0 1 
-1 .3000E+00 
3.6725E+0 1 
-1 .3000E+00 
3 . 5353E+0 1 
-1 . 3000E+00 
3 . 3773E+0 1 
-1 .300 OE + O 0 
3 . 2000E+0 1 
-1 . 3000E+00 
3 . 0047E+0 1 
-1 . 3000E+00 
2.7930E+0 1 
-1 . 3000E+00 
2 . 5667E+0 1 
-1 . 3000E+00 
2 . 3276E + 0 1 
-1 • 300 OE + O 0 
2 . 0776E+0 1 
-1 . 3000E+00 


1 . 5004E + 00 
1 • OOOOE+OO 
9.0566E-01 
1 .OOOOE+OO 
3* 2229E — 0 1 
1 .OOOOE+OO 
2.4515E-01 
1 .OOOOE+OO 
7 • 9231 E— 0 1 
1 .OOOOE+OO 
1 .3151 E+00 
1 .OOOOE+OO 
1 . 8096E+00 
1 .OOOOE+OO 
? . 2723E+00 
1 .OOOOE+OO 
2 . 6998E+00 
1 .OOOOE+OO 
3 • 0892E+00 
1 .OOOOE+OO 
3.4379E+00 
1 .OOOOE+OO 
3. 7435E+00 
1 .OOOOE + OO 
-4 • 0042E+ 00 
1 .OOOOE+OO 
— 4.2183E+00 
1 .OOOOE + OO 
— 4 • 3849E+00 
1 .OOOOE+OO 
— 4 . 5 030 E+00 
1 .OOOOE+OO 
— 4.5724E+00 
1 .OOOOE+OO 
-4.5931 E + 00 
1 .OOOOE+OO 
— 4.5654E+00 
1 .OOOOE+OO 
— 4.4902E+00 
1 . OOOOE+OO 
— 4 • 3685E+00 
1 . OOOOE+OO 
— 4.2019E+00 
1 . OOOOE+OO 
— 3.9922E+00 
1 . OOOOE+OO 
— 3.741 4E+00 
1 .OOOOE + OO 
-3.4521 E+00 
1 .OOOOE+OO 
-3. 1 270E+00 
1 .OOOOE+OO 
-2.7689E+00 
1 .OOOOE + OO 
— ? • 38 1 OE+00 
1 .OOOOE+OO 
- l . 9668E + 00 
I . OOOOE+OO 
- 1 . 5297E + 00 
1 .OOOOE+OO 


— 5.3840E+00 
-1 .4551 E+Ol 

— 5.3 590E+00 
-1 . 3695E+0 1 
— 5.2904E+00 
-1 .2701E+01 
-5. 1 794E+00 
-1.1 573E+01 
-5.0272E+00 
-1 .0319E+01 
— 4.8359E+00 
— 8 . 94 54 E+ 00 
— 4.6072E+00 
— 7.4605E+00 
— 4.3437E+00 
—5 • 8732E+ 00 
-4 .0477E + 00 
-4.1 934E+00 
— 3.7220E+00 
—2 .431 7E+00 

— 3.3 696 E+00 

— 5.991 3E— 0 1 

— 2 . 99 3 7 E+00 

1 . 2926E+ 00 

— 2 . 59 74 E+ 00 

3.2314E+00 

-2.1 841 E+00 
5.2049E+00 
-1 « 7574E + 0 0 
7.2006E+00 
-1 . 3206E+00 

9.2 059E+ 00 
-8.7740E-01 

1 . 1 208E + 0 1 
—4.3131 E— 0 1 
1 .31 94E+ 0 1 
1.41 38E-02 
1 .51 52E+0 1 
4 .5547E-0 1 
1 • 7070E+ 0 1 
8.8929E-01 
1 . 8936E + 0 1 
1.3! 23E+00 
2.0739E+01 
1 .721 3E+00 

2 . 2 468E + 0 1 
2.11 34E+00 
2.411 3E+01 
2 • 4856E+00 
2 . 5665E + 0 1 
2.8353E+00 
2.71 15E+01 
3. 1 601 E+00 
2 • 8456E+0 1 
3.4577E+00 
2 . 9680E+0 1 
3.7262E+00 
3.0781 E+Ol 
3 . 9638E + 00 
3. 1 756E+0 1 


2 . 7500 E+00 
3 . 7727E+00 
2 • 8000E+00 
3 . 8034E+00 

2 • 8500E+00 
3 . 8052E+00 
2 • 9000E+00 
3 . 7792E+00 
2.9500E+00 
3 . 7264 E+00 
3. OOOOE+OO 
3 . 6482E+00 
3 . 0500E+00 
3. 5460E+00 
3.1 0 OOE+OO 
3.421 5E+00 
3. 1 500E+00 

3 • 2766E+00 
3 . 2000E + 00 
3.11 33E+00 
3 . 2500E+00 
2 . 9336E+00 
3 • 3000E+00 
2 . 7397E+00 
3 . 3500E+00 
2 . 5339E+00 
3.4000E+00 
2.31 86E+00 
3.4500E+00 
2 . 0962E + 00 
3.5000E+00 
1 . B692E+00 
3.5500E+00 
1 • 6398E+00 
3.6000E+00 
1.41 06E+00 
3 . 65 00 E+00 
1.1841 E+00 
3.7000E+00 
9.6240E-01 
3.7500E+00 
7.4796E-01 
3 .8000E+00 
5 . 4293E— 0 1 
3 • 8500E+00 
3.4943E-01 
3 . 9000E+00 
1 .6943E-01 
3 .9500E+00 
4.8128E-03 
4. OOOOE+OO 

-1 .4271 E— 0 1 
4 . 0500 E + 00 
-2.7158E-01 
4 . 1 OOOE+OO 
— 3 . 8038E— 0 1 
4.1 500E+00 
—4 .67B8E— 01 
4 . 2000E+00 
-5.3302E-01 


-1.41 18E+00 
-1 .6782E+00 
-1 .940 1 E+00 
-2.1954E+00 
— 2.4422E+00 

— 2 . 6784E +00 
-2.9023E+00 
-3.11 22E+00 
-3. 3066E+00 
—3 • 4 84 OE+OO 

— 3.64 32 E+00 
— 3.7830E+00 
-3.9026E+00 
—4 .001 2E+00 
-4 . 0782E+00 
-4 . 1 332E+00 
—4 . 1 659E+00 
-4 . 1 763E + 00 
-4 . 1 645E+00 
-4. 1 309E+00 
-4 . 0758E+00 
— 3.9999E+00 
— 3.9039E+00 
— 3 . 78B9E+00 
-3.6557E+00 
-3.5058E+00 
— 3.3402E+00 
-3. 1 605E+00 
— 2.9681E+00 
-2.7646E+00 


APPENDIX B — Continued 


1 .6363E+00 
1 .OOOOE— 01 
3.2458E+00 
1 .OOOOE-ni 
4.8940E+00 
1 .OOOOE-01 
6.5742E+00 
1 .OOOOE-01 
8.2797E+00 
1 • OOOOE— 0 1 
1 • 0003E+0 1 
1 .OOOOE-01 
1 . 1 739E+01 
1 .OOOOE-01 
1 . 3479E+0 1 

1 .OOOOE-01 
1 .5217E+01 
1 .OOOOE— 0 1 
1 • 6946E+0 1 
1 .OOOOE-01 
1 .866 1 E+0 1 
1 • OOOOE— 0 1 
2 • 0356E+0 1 
1 • OOOOE— 0 1 
2.2024E+O1 
1 .OOOOE-01 

2.3660E+0 1 
1 . OOOOE— O 1 
2.5260E+01 
1 • OOOOE-o 1 

2.681 8E+0 1 
-2.3500E+00 
2.8322E+0 1 
— 4.8000E+00 
2 . 9727E+0 1 
-7.2500E+00 
3 . 098 1 E+0 1 
-9.7000E+O0 
3.2034E+01 
-9.5000E+O0 
3 • 2849E+0 1 
-9.3000E+00 
3. 3426E+0 1 
— 7.9000E+00 
3 • 37S3E+0 1 
-6.5000E+00 
3 . 3952E+0 1 
-5. 1 OOOE+OO 
3 • 3967E+0 1 
-3.7000E+O0 
3 . 3866E+0 1 
-2.3000E+00 
3. 3685E+0 1 
— 9 • OOOOE— 0 1 
3 . 3459E+0 1 
5 » OOOOE— O 1 
3.3225E+01 
1 • 9000E+00 
3.3019E+01 
3.3000E+00 
3.2874E+0 1 
4.7000E+00 
3 . 2825E+0 1 


1 .81 85E+01 
-1 • 3000E+00 
1 . 5525E+0 1 
-1 • 3000E+00 
1 .281 4E+0 1 
-1 • 3000E+00 
1 • 0074E+0 1 
-1 .3000E+00 
7 . 3254E+00 
-1 .3000E+00 
4 . 5875E + 0 0 
-1 • 3000E+00 
1 . 8806E+0 0 
-1 .3000E+00 
-7.7582E-01 
-1 . 3000E+00 
— 3 • 3630E+0 0 
-1 « 3000E+0 0 
-5.8629E+00 
-1 . 3000E+00 
— 8 • 25 84 E +00 
-1 .3000E+00 
-1 • 0533E+0 1 
-1 . 3000E+00 
-1 . 2673E+0 1 
-1 . 3000E+00 
-1 . 4663E+0 1 
-1 • 3000E+0 0 
-1 .6491 E+0 1 
-1 • 3000E + 0 0 
-1 .81468+01 
-1 • 3000E+0 0 
— 3.9592E+0 1 
-1 • 3000E+00 
-6.0442E+0 1 
-1 .3000E+00 
-8 . 0572E+0 1 
-1 . 3000E+00 
-9 . 9863E+0 1 
-9. OOOOE-O 1 
-9 . 4 956 E+0 1 
-5. OOOOE-O 1 
-8.951 0E+01 
-5. OOOOE-O 1 
-7 . 5464E+0 1 
-5 .OOOOE— 0 1 
-6.1 161E+01 
-5 • OOOOE— 0 1 
-4 . 6695E+0 1 
-5. OOOOE-O 1 
-3.2160E+01 
-5. OOOOE-O 1 
-1 • 7652E+0 1 
-5 . OOOOE— 0 1 
-3 . 2624 E+0 0 
-5. OOOOE-O 1 
1 .0915E + 01 
-5. OOOOE-O 1 
2.4792E+0 1 
-5. OOOOE-O 1 
3.8281 E+0 1 
-5. OOOOE-O 1 
5.1 300E+0 1 


-1 . 0733E+ 00 
1 . OOOOE+OO 
—6 .0151 E— 0 1 
1 .OOOOE+OO 
-1 . 1806E-01 
1 .OOOOE+OO 
3.731 5E— 0 1 
1 .OOOOE+OO 
8 .6823E-0 1 
1 .OOOOE+OO 
1 . 3633E+00 
1 .OOOOE+OO 
1 .8546E+00 
1 .OOOOE+OO 
2. 3383E+00 
1 .OOOOE+OO 
2.81 08E+00 
1 .OOOOE+OO 
7. 2685E+00 
1 . OOOOE+OO 
3 . 7083E+00 
1 . OOOOE+OO 
4. 1 268E+00 
1 .OOOOE+OO 
4.521 1E+00 
1 .OOOOE+OO 
4.8885E+00 
1 .OOOOE+OO 
5 • 2265E+ 00 
1 .OOOOE+OO 
8. 5328E + 00 
1 .OOOOE+OO 
8 • 5789E+ 00 
1 .OOOOE+OO 
1 . 1 536E+01 
1 .OOOOE+OO 
1 • 4379E+ 0 1 
1 .OOOOE+OO 
1 • 7086E+0 1 
1 .OOOOE+OO 
1 .51 40E+0 1 
1 .OOOOE+OO 
1 . 3098E + 0 1 
1 .OOOOE+OO 
1.11 10E+01 
1 .OOOOE+OO 
9. 0563E + 00 
1 .OOOOE+OO 
6 • 9558E+00 
1 .OOOOE+OO 
4 . 8261 E+00 
1 .OOOOE+OO 
? . 6853E+ 00 
1 .OOOOE+OO 
5*5161 E— O 1 
1 .OOOOE+OO 
-1 . 5574E+00 
1 .OOOOE+OO 
■3 • 6242E+00 
1 .OOOOE+OO 
■5. 6322E+00 
1 .OOOOE+OO 
-7 . 5649E+ 00 


4.1 692E+00 
3 . 2599E+0 1 
4 • 34 1 OE+OO 
3 • 3307E+0 1 
4.4784E+00 
3.3880E+01 
4.5807E+00 
3.431 5E+ 0 1 
4.6476E+00 
3.4613E+01 
4 .6 790E+00 

3 . 4 774E+0 1 
4.6 752E+00 
3.4801 E+ 0 1 
4 . 6365E + 00 
3 • 4698E + 0 1 
4 . 5637E+00 

3 . 4 4 67E + 0 1 

4 . 4 579E+ 00 
3.41 13E+01 
4 • 3202E+ 00 
3 . 3643E+ 0 1 
4.1 52 1 E+00 
3 . 3062 E+0 1 
3 . 9553E+00 
3.2378E+01 
3.7317E+00 
3.1 598E+ 0 1 

3 • 4834E+00 
3 . 0732E+0 1 
3.21 26E+00 
2.9286E+01 
2.81 05E+00 
2.6783E+01 
2.1 709E+00 
2.3255E+01 
1 . 30 1 OE+OO 
1 . 8740E+0 1 
2.1011 E— 0 1 

1 . 3868E+0 1 
•9.3039E-01 
9.2538E+00 
■1 .951 3E+00 
5.1 280E+00 
■2.8332E+00 
1 .71 13E+00 
■3.5574E+00 
9 . 8575E— 0 1 
■4.11 99E + 00 
2 . 9574E+00 
■4.51 79E+00 

4 .2025E+00 
4.7503E+00 
4 . 7249E+ 00 
4.81 68E+00 
4 .5326E+00 
4.71 87E+00 
3.6387E+00 
4 .4587E+00 
■2 . 0603E+00 
4 . 0403E+00 
1 .8 1 22E-0 1 
3 . 46B6E+00 


4.2500E+00 
—5.74 95E— 0 1 
4 . 3000E+00 
— 5 . 9298 E — 0 1 
4 . 35 00 E +00 
—5 . 8663E-0 1 
4 • 40 OOE + OO 
— 5 • 5560 E—0 1 
4 .4500E+00 
-4.9981 E—0 1 
4 .5000E+00 
-4.1933E-01 
4 . 5500E+00 
-3. 1447E-01 
4 .6000E+00 
-1 . 8569E— 0 1 
4 . 6500E+00 
— 3 . 3644 E — 0 2 
4.7000E+00 
1 . 4086E— 0 1 
4 . 7500E+00 
3 . 3683E-0 1 
4 .8000E+00 
5.531 3E-01 
4 • 8500E+00 
7.8849E-01 
4 . 9000E+00 
1 . 04 1 5E+00 
4 . 9500E+00 
1 .31 06E+00 
5. OOOOE+OO 
1 .6637E+00 
5 • 0500 E+00 
2. 1 670 E+00 

5.1 OOOE+OO 
2.81 53E+00 
5. 1500E+00 
3 . 6025E+00 
5 • 2000E+00 
4 . 4086E+00 
5 . 2500E+00 
5.11 49E+00 
5.3OOOE+O0 
5.7205E+00 
5 • 3500E+00 
6 • 2248E+00 
5.4000E+00 
6 . 6253E+00 
5 • 4500E+00 
6 • 9200E+00 
5 . 5000 E+00 
7.1 078E+00 
5. 5500E+00 
7. 1 887E+00 
5 • 6000E+00 

7.1 634 E+00 
5 . 6500E+00 
7 . 0337E+00 
5 • 7000E+00 
6 • 8020E+00 
5.7500E+00 
6.471 7E+00 
5 . 8000E+00 


-2.551 7E+00 
-2.331 1 E+00 
-2. 1 045E+00 
-1 . 8736E+00 
-1 . 64 03E +00 
-1 .4063E + 00 
-1 . 1734E+00 
-9.4323E-01 
-7. 1755E-01 
—4 . 9797E— 0 1 
-2.8604E-01 
— 8 . 3234E— 02 
1 .0905E-01 
2.8953E-01 
4.5702E— 01 
6 . 0836E-0 1 
7 . 3364E — 0 1 
8.21 15E-01 
8.5959E-01 
8.4121 E—0 1 
7.6876E-01 
6 • 4963E— 0 1 
4 . 8833E— 0 1 
2.9585E-01 
7 • 9348E— 02 
-1 . 5292E— 0 1 
-3.9265E-01 
-6.3158E-01 
-8.6155E-01 
-1 . 0745E+00 
-1 .2628E+00 
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6.1 OOOE+OO 
3.2904E+01 
7.5000E+00 
3.31 42E+0 1 
8.9000E+00 
3.3570E+01 
1 . 0300E+0 1 
3.4214E+01 
1 . 0400E+0 1 
3 . 5095E+0 1 
1 . 0500E+0 1 
3.6218E+01 
1 . 0600E+0 1 
3 . 7576E+0 1 
1 . 0700E+0 1 
3.91 65E+0 1 
9.5000E+00 
4 . 0972E+0 1 
8.3000E+O0 
4.2967E+01 
7.1 OOOE+OO 
4.51 15E+01 
5 . 9000E+O0 
4.7379E+01 
4 • 7000E+00 
4 • 9722E+0 1 
3.5000E+n0 
5.21 08E+01 
2.3000E+00 
5 • 4499E+0 1 
1 . 1 OOOE+OO 
5.6860E+O 1 

— 1 • OOOOE— 0 1 
5.9155E+01 

-1 • 3000E+ 00 
6.1 349E+0 1 
-2.5000E+O0 
6.341 OE+Ol 
— 3.7000E+O0 
6 . 5304E+ 0 1 

— 3 . 7 OOOE+OO 
6 • 7009E+0 1 

— 3.7000E+00 
6 . 8526E+0 1 
-3.7000E+O0 
6.9868E+0 1 
-3.7000E+n0 
7. 1 044E+01 
-3.7000E+00 
7.2063E+01 
-3.7000E+00 
7.2924E+0 1 
-3.7000E+00 
7 . 3623E+0 1 
-3.7000E+O0 
7.41 58E+0 1 
-3.7000E+00 
7.4523E+0 1 
-3.7000E+00 
7.471 6E+0 1 

— 3 » 7000E+O0 
7 • 4734E+0 1 

—3. 7 OOOE+nO 


-5.0000E-0 1 
6 . 3770E+0 1 
-5.0000E-01 
7.5618E+01 
-5 .OOOOE - 0 1 
8.6775E+0 1 
-5 • OOOOE- 0 1 
9.71 79E+0 1 
-5 .OOOOE- 0 1 
9.61 76E+0 1 
-5 .OOOOE— 0 1 
9.4529E+0 1 
-5 • OOOOE— 0 1 
9.2258E+01 
-5.0000E-0 1 
8.9388E+0 1 
—5 • OOOOE— 0 1 
7 . 5349E+0 1 
—5 .OOOOE— 0 1 
6 . 0987E+0 1 
—5 .OOOOE— 0 1 
4.6402E+01 
—5 . OOOOE— 0 1 
3.1 692E+0 1 
-5.0000E-0 1 
1 .6959E+01 
—5. OOOOE— 01 
2 • 3003E+00 
—5 • OOOOE— 0 1 
-1 .21 85E+0 1 
-5.0000E-01 
-2.6403E+0 1 
—5 .OOOOE— 0 1 
-4.0259E+0 1 
—5 . OOOOE— 0 1 
— 5 • 3666E+0 1 
—5 .OOOOE— 0 1 
— 6 • 6539E+0 1 
—5 . OOOOE— 0 1 
—7 . 8798E+0 1 
9 • OOOOE— 0 1 
—7 . 4825E+0 1 
2 . 3000E+0 0 
—7 . 0550E+0 1 
3 • 7000E+0 0 
-6.61 19E+01 
5.1 OOOE+OO 
—6 • 1 677E+0 1 
5 • 0900E+00 
—6.31 64E+0 1 
5 . 0800E+00 
— 6.4655E+0 l 
5.0700E+00 
-6.61 25E+0 1 
5 . 0600E+00 
— 6 . 7549E+0 1 
5 . 0500E+0 0 
—6 • 8902E+0 1 
5 . 0400E+00 
-7.0161E+01 
5 . 0300E+0 0 
-7. 1 306E+01 
5.0200E+00 


1 • OOOOE+OO 
— 9.4072E+00 
1 .OOOOE+OO 
-1.11 44E+0 1 
1 .OOOOE+OO 
-1 .2763E+01 
1 .OOOOE+OO 
-1 • 4250E+0 1 
1 .OOOOE+OO 
-1 .41 23E+01 
1 .OOOOE+OO 
-1 . 3872E+0 1 
1 .OOOOE+OO 
-1 . 350 1 E+0 1 
1 .OOOOE+OO 
-1 . 30 1 3E+0 1 
1 .OOOOE + OO 
-1 . 0942E+0 1 
1 .OOOOE+OO 
-8.7943E+00 
1 .OOOOE+OO 
— 6 • 5886E+00 
1 .OOOOE+OO 
-4.3440E+00 
1 . OOOOE+OO 
— 2 • 0794E+ 00 
1 .OOOOE+OO 
1 . 8606E— 0 1 
1 .OOOOE + OO 
2.4335E+00 
1 .OOOOE+OO 
4 . 6445E+00 
1 .OOOOE+OO 
6.S012E+00 
1 .OOOOE+OO 
8 . 8860E+00 
1 .OOOOE+OO 
1 . 0883E+0 1 
1 . OOOOE + OO 
1 . 2775E + 0 1 
1 .OOOOE+OO 
7. 9600E+00 
1 .OOOOE+OO 
3« 1 239E+00 
1 .OOOOE+OO 
-1 . 7055E+00 
1 .OOOOE+OO 
-6.501 1E+00 
1 .OOOOE+OO 
— 5 • 9687E+00 
1 .OOOOE+OO 
—5 . 4368E+00 
1 .OOOOE+OO 
—4 .91 04E+00 
1 .OOOOE+OO 
— 4 • 3942E+00 
1 .OOOOE+OO 
— 3.8929E+00 
1 .OOOOE+OO 
-3*41 09E+00 
1 .OOOOE+OO 
— 2 • 9524 E+ 00 
1 .OOOOE+OO 


3 • 0602E+ 00 
2 • 7496E+00 
6 • 5475E+00 
1 • 8904E+00 
1 .061 OE + O 1 
•8.9924E-01 
1 .52 1 2E+0 1 
2. 1490E-01 
2 . 0049E+0 1 
1 .3832E+00 
2.481 9E+ 0 1 
2.5371 E+00 
2 .949 1 E+ 0 1 
3 • 6680E+ 00 
3.4035E+01 
4 • 7673E+00 
3.81 55E+ 0 1 
5 . 7679E+ 00 
4. 1 565E+01 
6.6045E+00 
4 • 4 250E + 0 1 

7 . 2 722E + 00 
4 .6203E+01 
7.7677E+00 
4 • 74 1 9E+ 0 1 
8 . 0891 E+00 
4 . 7900E + 0 1 

8.2 358E+ 00 
4 . 7652E+0 1 
8.2086E+00 
4.6687E+01 
8 . 0097E+00 
4.5018E+01 
7.6425E+00 
4 . 2668E+0 1 
7.111 7E+00 
3.9661 E+01 
6.4233E+00 
3 . 6 025E+ 0 1 
5.5844E+00 
3.2 1 83E+01 
4 .7999E+00 
2 . 8547E+0 1 
4 • 2 695E + 0 0 
2.51 30E+0 1 
3.991 4E+ 00 
2.1 935E+0 1 
3 . 9620E+00 
1 • 88 1 5E+0 1 
4 .0321 E+00 
1 .5619E+01 
4 . 0562E+00 

1 . 2 34 9E+ 0 1 
4.0348E+00 
9 • 0073E+00 
3.9691 E+00 
5 . 5957E+00 
3.8601 E+00 
2.1 187E+00 
3.7096E+00 

-1 .41 85E+00 
3.51 94E+00 
—5 • 0096E+00 


6 . 0470E+00 • 

5 • 8500E+00 

5 . 5328E+00 ■ 
5 • 9000E+00 
4 . 9346E+00 • 
5.9500E+00 
4 . 25B8E+00 ■ 
6. OOOOE+OO 
3 . 5489E+00 

6 • 0500 E+00 
2 .8485E+00 
6.1 OOOE+OO 
2.1 637E+00 
6. 1 500E+00 
l . 5004 E+00 
6 . 2000E+00 
9.01 15E-01 
6 . 2500E+00 
4.0744E— 01 
6 • 3000E+00 
2 . 265 1 E— 02 
6 . 3500E+00 

—2 . 5O0OE— 0 1 
6 . 4000E+00 
-4.11 44E-0 1 
6. 4500E+00 
— 4.5876E— 01 
6 . 5000E+00 
-3.931 7E-01 
6 . 5500 E + 00 
-2. 1604E-01 
6 . 6000 E+00 
7 . 0346E-02 
6.6500E+00 
4.6284E-0 1 
6 . 7000E+00 
9.5745E-0 1 
6 • 7500 E+00 
1 . 5493E+00 
6 . 8000E+00 
2 . 0679E + 00 
6 • 8500E+00 
2.3451 E+00 
6 . 9000E+00 
2.3805E+00 
6 . 9500E+00 
2. 1 752E+00 
7. OOOOE+OO 
1 • 8635E + 00 
7. 0500E + 00 
1 . 5784E+00 
7. 1 OOOE+OO 
1 .31 97E+00 
7. 1 500E+00 
1 • 0872E+00 
7 . 2000 E+00 
8.8004E-01 
7.2500E+00 
6.9754E-01 
7 • 3000E+00 
5. 3856E-0 1 
7 . 3500E+00 
4.01 84 E— 0 1 


1 .4187E+00 
1 .5351 E+00 
1 • 6053E+00 
1 .6227E+00 
1 . 5827E +00 
1 .4846E+00 
- 1 . 3294E +00 
-1.1 183E+00 
— 8.5440E— 01 
—5 . 445 1 E— 0 1 
-1 .9700E-01 
1 .796 1 E— 0 1 
5.7664E-01 
9.8538E-01 
1 . 397 1 E+00 
1 .8032E+00 
2. 1 950E+00 
2.564 5E + 00 
2.9034E+00 
3 • 204 OE + 00 
3.4629E+00 
3.6889E+00 
3 • 8946E + 00 
4 . 0928E+00 
4 . 292 8E + 00 
4 .4952E+00 
4 .6977E+00 
4 . 8979E + 00 
5.0938E+00 
5 . 2832E + 00 
5.464 1E+00 
5.6346E+00 
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7.4573E+01 
-3.7000E+00 
7 . 4232E+0 1 
-3.7000E+00 
7.3708E+01 
-3.7000E+00 
7.2999E+01 
-3.7000E+O0 
7.21 05E+01 
-3.7000E+00 
7.1023E+01 
-3.7000E+00 
6.9755E+0 1 
-3.7000E+00 
6.8301 E+O 1 
-3.7000E+00 
6.6661 E+Ol 
-3 » 7000E+ OO 
6 • 4837E+0 1 

— 3.7000E + nO 
6.2831 E+Ol 

— 3.7000E+O0 
6 • 0645E+O 1 
-3.7000E+n0 
5 . 8283E + 0 1 
— 3.7000E+P0 
5.5746E+0 1 

— 3.7000E+00 
5 . 3026E+0 1 

— 3.7000E+00 
5.011 3E + 01 

-3.7000E+O0 
4 .6999E+0 1 
-3.4429E+00 
4.3681E+01 
-3. 1 857E+00 
4.01 74E+0 1 

— 2 • 9286 E+OO 
3 . 6497E+0 1 

-2.6714E+n0 
3.2669E+0 1 
-2.4143E+O0 
2 .871 4E+0 1 
-2.1 571 E+OO 
2.4655E+0 1 

— 1 .90O0E + n0 
2.0516E+01 

— 1 . 6429E + 00 
1 .6322E + 0 1 

— 1 . 3857E + 00 
1 .2099E+O 1 

-1.1 286E+nO 
7.8726E+O0 
-8.7143E-01 
3.6687E+00 
-6. 1 429E-01 
-4.8689E— nl 
—3.571 4E— 0 1 

— 4 • 5686E + 00 

— 1 • OOOOE — O 1 
— 8.5514E+00 


-7.2315E+01 
5.01 00E+00 
—7 .31 70E+0 1 
5.0000E+00 
— 7 • 3854E+0 1 
4.9900E+00 
-7.4353E+01 
4 . 9800E+0 O 
—7 . 4652E+0 1 
4.9700E+00 
-7.4740E+0 1 
4 .9600E + 00 
—7 .46 1 0E+0 1 
4.9500E+O0 
-7.4252E+01 
4.9400E+00 
— 7.3663E+01 
4 . 9300E + 0 0 
-7.2839E+0 1 
4 .9200E + 00 
-7. 1 780E+0 1 
4.91 00E+00 
—7 . 0486E+0 1 
4.9000E+00 
-6.8961 E+0 1 
3.4000E+00 
—7 . 334 1 E+0 1 
1 .9000E+00 
-7.7231E+01 
4 . OOOOE— 0 1 
—8 . 0479E+0 1 
-1.1 000E+00 
—8 . 294 1 E+0 1 
-1 .0357E+00 
— 7 • 5956E+0 1 
-9.71 43E-0 1 
—6 .8263E+0 1 
-9.0714E-01 
— 5 . 9939E+0 1 
— 8.4286E— 0 1 
-5. 1 062E+0 1 
-7.7857E-0 1 
-4. 1 71 6E+0 1 
-7. 1 429E-0 1 
-3.1 985E+0 1 
-6.5000E-0 1 
-2. 1 957E+0 1 
-5.8571 E-0 1 
-1 . 1 71 9E+0 1 
-5.21 43E-0 1 
-1 . 3592E + 0 0 
-4.571 4E-0 1 
9 . 0 339E+0 0 
-3 .9286E-0 1 
1 .9374E+0 1 
—3 . 2857E— 0 1 
2.9575E+0 1 
-2.6429E-0 1 
3.9555E+01 
-2.0000E-0 1 
4.9235E+01 


— 2 • 52 1 3E+00 
1 . 0000E+00 
—2 .121 2E+00 
1 . OOOOE+OO 
-1 . 7554E+00 
1 .OOOOE + OO 
-1 . 4269E + 00 
1 .OOOOE+OO 
-1 . 1383E + 00 
1 .OOOOE+OO 
—8*91 97E— 0 1 
1 .OOOOE+OO 

— 6 • 8964E — 0 1 
1 .OOOOE+OO 

— 8 . 3285E— 0 1 
1 • OOOOE+OO 
— 4.2267E— 01 
1 .OOOOE+OO 

— 7 . 5980E — 0 1 
1 .OOOOE + OO 

-3.4452E-01 
1 .OOOOE+OO 
-3.7673E-0 1 
1 .OOOOE+OO 
— 4 . 5595E— 0 1 
1 .OOOOE+OO 
4 . 9857E+00 
1 .OOOOE+OO 
1 . 0292E + 0 1 
1 .OOOOE + OO 
1 . 5434E+0 1 
1 .OOOOE+OO 
2 • 0384E+ 0 1 
1 .OOOOE + OO 
1 . 8982E + 0 1 
1 .OOOOE+OO 
1 . 744 1 E+0 1 
1 .OOOOE + OO 
1 . 5775E+0 1 
1 .OOOOE+OO 
1 . 3999E+0 1 
1 • OOOOE+OO 
1 .21 28E+01 
1 .OOOOE+OO 
].01 79E+01 
1 . OOOOE+OO 
8. 1 667E+00 
1 . OOOOE+OO 
6. 1093E+00 
1 .OOOOE+OO 
4 . 0228E+00 
1 . OOOOE+OO 
1 .924 1 E+OO 
1 . OOOOE+OO 
-1 .7034E-01 
1 . OOOOE+OO 
—2 • 244 1 E+OO 
1 • OOOOE+OO 
— 4.2812E+00 
1 . OOOOE+OO 

— 6 • 2662E + 00 


3.291 7E+00 

— 8 • 64 74 E+ 00 
3. 0290E+00 

-1 .2324E+0 1 
2 . 7338E + 00 
-1 • 6030E+ 0 1 
2.4 092E+00 
-1 . 9756E+ 0 1 
2.0581 E+OO 

— 2 . 3491 E+ 0 1 
1 . 6838E+ 00 

-2.7226E+01 
1 . 2897E + 00 

— 3 • 0949E+ 0 1 
8.7925E-01 

-3.4647E+01 
4 .8601 E-01 
-3.831 1 E+Ol 
2.3589E-02 
-4. 1 927E+01 
-4.1438E-01 
—4 • 5485E + 0 1 
— 8.5428E— 01 
— 4.8972E+01 
-1 . 2925E+ 0 0 
-5.2531 E+Ol 
-1 .8771 E + OO 
— 5.6298E+0 1 
— 2.75C7E+00 

— 6 . 0244E + 0 1 

— 3 . 8997 E+ 00 
-6 . 4332E+0 1 
-5.3088E+00 
“6.8308E+0 1 
— 6 . 7898E + 00 
-7. 1916E+01 
-8. 1 593E+00 
— 7.5124E+01 
-9. 4081 E+OO 
-7.7901 E+0 1 
- 1 . 0528E+0 1 
—8 . 0222E+0 1 
-1 .151 3E+0 1 

— 8 . 2 066E+ 0 1 
-1 .2356E + 0 1 
-8.341 6E+0 1 
-1 • 3053E+0 1 

— 8 . 4 259E + 0 1 
-1 .360 1 E + 0 I 

— 8 • 4 586E+ 0 1 
-1 • 3997E + 0 1 

— 8 . 4 394 E + 0 1 
- 1 .424 1 E+0 1 
-8.3683E+01 
-1 . 4 332E+ 0 1 
— 8.2459E+01 
- 1 . 4272E+0 1 
—8 . 0730E+0 1 
-1 . 4 064 E+0 1 
—7. 8509E+0 1 
-1 .371 1 E+Ol 


7 . 4000E+00 
2.8591E-01 
7 . 4500E+00 
1 .8915 E— 0 1 
7 . 5000E+00 
1 . 0975E— 0 1 
7 . 5500E+00 
4 .5794E-02 
7 . 6000 E+OO 
— 4.7846E— 03 
7 • 6500E+00 
-4.41 38E-02 
7.7000E+00 
— 7.4509E— 02 
7 . 7500E+00 
-9.8201 E-02 
7 . 8000E+00 
-1.1 756E-0 1 
7. 8500E+00 
-1 . 3497E— 0 1 
7 . 9000E+00 
-1 . 5281 E-0 1 
7 .95 00 E+OO 
-1 . 7343E— 0 1 
8. OOOOE+OO 
— 5.9748E— 02 
8 • 0500E+00 
3 . 2276E-0 1 
8. 1 OOOE+OO 
9.6659E-0 1 
8. 1500E+00 
1 • 862BE+00 
8.2000E+00 
2 • 8476E+00 
B.2500E+00 
3 . 7587E+00 
8 . 3000E+00 
4 . 5896E+00 
8 . 3500E+00 
5 . 3344 E + OO 
8 . 4000E+00 
5 . 9879E+00 
8.4500E+00 
6 . 5459E+00 
8 . 5000E+00 
7 . 0047E+00 
8 . 5500E+00 
7.361 8E+00 
8 . 6000E+00 
7.61 52E+00 
8 . 6500E+00 
7.7639E+00 
8 . 7000E+00 
7 . 8077E+00 

8 . 7500E+0O 
7.7472E+00 
8.8000E+00 
7 • 5839E+00 
8 . 8500E+00 
7.31 99E + 0 0 
8 . 9000 E + OO 


5 . 7927E+00 
5.9369E+00 
6. 0656E+00 
6.1 774E+00 
6.271 OE+00 
6.3454E+00 
6 . 3997E+00 
6.4331E+00 
6 . 4452E+00 
6.4354E+00 
6.4037E+00 
6 . 3500E+00 
6.271 7E+00 
6. 1 568E+00 
5.991 4E+00 
5,761 9E +00 
5.4590E+00 
5 . 0848E+00 
4 .6452E+00 
4 . 1 462E+00 
3 . 5947E+00 
2 . 9974E+00 
2.3616E+00 
1 .694 7E +00 
1 . 004 1 E+OO 
2 . 9759E —0 1 
-4 . 1 732E-01 
-1.1 330E+00 
- 1 . 84 20E+00 
-2.5369E+00 


56 
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-9.5455E-02 
-1 .2412E+01 
-9.0909E-02 

— 1 .61 33E+01 
— 8.6364E— 02 
-1 . 9696E+0 1 
-8. 1 818E-02 

— 2 . 3087E+O 1 
-7.7273E-02 
-2.6291 E+Ol 
— 7.2727E— 02 
— 2.9296E+0 1 
-6.8182E-02 

— 3.2089E + 0 1 
— 6.3636E— 02 
— 3 • 4662E+0 1 
— 5.9091E— 02 

— 3 • 7006E+O 1 
-5.4545E— 02 
-3.91 13E+01 

— 5.0000E— 02 
-4. 0980E+0 1 
-4.5455E-02 
-4.2602E+0 1 
— 4 • 0909E— 02 

— 4 • 3978E+0 1 
— 3.6364E— 02 
-4.51 09E+0 1 
—3.181 8E— 02 
-4.5994E+01 
— 2.7273E— 02 

— 4 • 6639E+0 1 

— 2 • 2727E— 02 
-4.7047E+01 
-1 .8182E-02 
—4 . 7226E+0 1 

— 1 . 3636E — 02 
—4.71 83E+0 1 
— 9.0909E— 03 

— 4 . 6928E + 0 1 
— 4.5455E— 03 

— 4 . 6472E + 0 1 
-5.6932E-1 3 
-4.5825E+01 


1 .9091 E — 0 1 
5.6250E+0I 
1 .81 82E-01 
6 • 2867E+0 1 
1 • 7273E— 0 1 
6.9036E+01 
-1 • 6364E— 0 1 
7.471 3E+0 1 
-1 .5455E— 01 
7.9855E+0 1 
-1 .4545E-0 1 
8 • 4426E+0 1 
-1 . 36 36E— 0 1 
8 • 8396E+0 1 
-1 . 2727E— 0 1 
9. 1738E+01 
-1 . 1818E-01 
9 • 4432E+0 1 
-1 . 0909E— 0 1 
9.6461E+01 
-1 • OOOOE— 0 1 
9.7817E+01 
— 9.0909E— 02 
9 . 8493E+0 1 
-8. 1818E-02 
9.8492E+0 1 
-7.2727E— 02 
9.7819E+0 1 
-6.3636E-02 
9.6485E+01 
— 5.4545E— 02 
9 . 4508E+0 1 
—4 . 54 55 E— 02 
9. 1907E + 01 
— 3.6364E— 02 
8 • 8709E+0 1 
— 2.7273E— 02 
8.4944E+0 1 
-1 .81 82E-02 
8.0646E+0 1 
— 9.0909E— 03 
7.5854E+0 1 
-1.1 386E-1 2 
7 . 0607E+0 1 


1 • OOOOE+OO 

— 7.6922E + 00. 
1 .OOOOE+OO 

— <9. 0462E+00 
1 . OOOOE+OO 

-1 .0318E+01 
1 .OOOOE+OO 
-1 . 1 500E+0 1 
1 .OOOOE+OO 
-1 . 2582E+0 1 
1 .OOOOE+OO 
-1 . 3558E+0 1 
1 .OOOOE+OO 
-1 . 4421 E+0 1 
1 .OOOOE+OO 
-1 .5166E+01 
1 . OOOOE+OO 
-1 .5789E+0 1 
1 • OOOOE+OO 
-1 . 6285E+0 1 
1 . OOOOE+OO 
- l . 6654E+0 1 
1 . OOOOE+OO 
- l . 6893E+0 1 
1 .OOOOE+OO 
-1 . 7002E+0 1 
1 . OOOOE+OO 
-1 .6981 E+01 
1 .OOOOE+OO 
-1 • 6834E+0 1 
1 . OOOOE + OO 
-1 .656 1 E + 0 1 
1 .OOOOE + OO 
-1 .6168E+01 
1 .OOOOE+OO 
-1 . 5658E + 0 1 
1 .OOOOE + OO 
-1 . 5036E+0 1 
1 . 00005 + 00 
-1 .431 05 + 01 
1 .00005 + 00 
- 1 . 34855+01 
1 • OOOOE+OO 
-1 . 2570E+0 1 


-7.5870E+01 
-1 . 3235E+0 1 
—7 . 2890E+0 1 
-1 . 2 659E+0 1 
-6.9591 E+01 
-1 .1987E+01 
-6.5995E+01 
-1 . 1 228E+0 1 
—6 .21 28E+0 1 
-1 .03S7E+01 
— 5.8019E+01 
— 9.4723E+00 
— 5 . 3696E+0 1 
— 8.4924E+00 
-4.9190E+01 

— 7.4 556E+00 
— 4.4533E+01 

— 6 . 3706E + 00 
— 3 . 9758E+0 1 
— 5.2467E+00 

— 3.4898E + 0 1 
-4 . 0930E+00 

— 2 . 9987E + 0 1 
-2.91 90E+00 
-2.5060E+01 
-1 . 7340E+00 
-2.01 49E + 0 1 

— 5.47558 — 0 1 
-1 .5289E + 0 1 

6.3121E — 01 
- 1 .051 1 E+01 
1 .79315+00 

— 5 • 8485E+ 00 
2 . 9294 E+ 00 

-1 • 3307E+00 
4 .0315E+00 
3.01 29E+00 
5.091 3E+00 
7. 1 5485+00 
6.1011 5+00 
1 . 1 0695+01 
7.05365+00 
1 . 4733E + 0 1 
7 .9422E + 00 


6 • 9707E+00 
8.9500E+00 
6.551 9E+00 
9. OOOOE + OO 
6 . 0675E+00 
9.0500E+00 
5.5216E+00 
9. 1000E+0P 
4.9191 E+00 
9. 1 500E+00 
4 . 2652 E+00 
9.2000E+00 
3 . 5653E+00 
9 . 2500E+00 
2.8251 E+00 
9 . 3000E+00 
2.0507E+00 
9 • 3500E+00 
1 . 2483E+00 
9 . 4000E+00 
4 .2431 E-0 1 
9 . 4500E+00 
-4 . 1 489E-0 1 
9 . 5000E+00 
-1 . 2628E+00 
9 . 5500E+00 
-2.11 29E+00 
9.6000E+00 
— 2 . 9588 E+00 
9 • 6500E+00 
—3 • 7942E+00 
9 . 7000E+00 
-4.61 29E+0O 
9 . 75 00 E+00 
-5. 4090E+00 
9 . 8000E+00 
-6 • 1 768E+00 
9 . 8500E + 00 
-6.91 09E+00 
9 . 9000E+00 
-7.6062E+00 
9 . 9500E+00 
-8.2579E+00 
1 . OOOOE + O 1 


-3.21 10E+00 
-3 .8588E+00 
—4 .4753E+00 
—5 . 056 1 E + 00 
—5 • 5967E + 00 
-6 . 0935E+00 
—6 . 54 29E + 00 
—6.941 8E+00 
-7.2876E+00 
— 7 . 5782E+00 
-7.81 18E+00 
—7 .9871 E+00 
-8 • 1 035E + 00 
-8 . 1 605E+00 
-8 . 1 584 E + 00 
-8. 0977E+00 
— 7.9795E + 00 
-7.B053E+00 
-7.5770E+00 
-7.2970E+00 
-6.9679E+00 
— 6.5927E+00 


U1 

co 


Output Listing 


THERE ARE 200 DATA POINTS, 10 ITERATIONS AND THE TIME INCREMENT IS .0500 SECONDS 
THE SCALING PARAMETER, FINKAS VALUE IS .0100 


MATRIX A 

- 0 . 

- 0 . 

- 0 . 

1 .00 OOE+OO 
MATRIX B 
- 0 . 

- 0 . 

- 0 . 


HAS A ROWS AND 

- 0 . - 0 . 

- 0 . - 0 . 

-1 . OOOOE+OO -0. 

-0. -C. 

HAS 4 ROWS AND 

- 0 . - 0 . 

-0. -C. 

- 0 . - 0 . 


4 COLUMNS 

- 0 . 

- 0 . 

7.8900E-02 

- 0 . 

3 COLUMNS 


MATRIX D1S 
l.OOOOE+OO 
l.OOOOE+OC 
MATRIX F 
0. 

0. 

- 0 . 

MATRIX G 

0. 

0 . 

0. 

MATRIX SZ 
L.OOOOE+OO 
l.OOOOE+OO 
MATRIX SU 
l.OOOOE+OO 
MATRIX 01 
l.OOOOE+OO 
l.OOOOE+OO 


HAS 1 ROWS 

l.OOOOE+OO 1 
l.OOOOE+OO 
HAS 3 ROWS 

0. 0 
0. 0 
-l.OOOOE+OO -0 
HAS 3 ROWS 

0. 0 
0. 0 
0. C 

HAS 1 ROWS 

l.OOOOE+OO 1 
l.OOOOE+OO 
HAS 1 ROWS 

l.OOOOE+OO 1 
HAS 1 ROWS 

l.OOOOt+OO 1 
l.OOOOE+OO 


ANO 7 COLUMNS 
•OOOOE+OO l.OOOOE+OO 


ANO 


ANO 


4 COLUMNS 

0. 

0 . 

7.8900E-02 
3 COLUMNS 


ANO 7 COLUMNS 

.OOOOE+OO l.OOOOE+OO 

AND 3 COLUMNS 

.OOOOE+OO 
ANO 7 COLUMNS 

.OOOOE+OO l.OOOOE+OO 


l.OOOOE+OO 


l.OOOOE+OC 


l.OOOOE+OO 
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PARAMETEK 

GROUP 

ROW 

COLUMN 

apr i or I DEVIATION 

SCALE 

1 

1 

1 

1 

- 0 . 000099999.0000 

1.0000 

-0 

5 

1 

1 

- 0.0000 

- 0.0000 

1.0000 

1 

1 

i 

2 

- 0 . 000099999.0000 

1.0000 

-0 

5 

1 

2 

- 0.0000 

- 0.0000 

1.0000 

1 

1 

1 

3 

- 0 . 000099999 . 000 C 

1.0000 

-0 

5 

1 

3 

- 0.0000 

- 0.0000 

1.0000 

1 

2 

1 

1 

- 0 . 000099999.0000 

1.0000 

-0 

6 

1 

1 

- 0.0000 

- 0.0000 

1.0000 

1 

2 

1 

2 

- 0 . 000099999.0000 

1.0000 

-0 

6 

1 

2 

- 0.0000 

- 0.0000 

1 .0000 

1 

2 

1 

3 

- 0 . 000099999.0000 

1.0000 

-0 

6 

1 

3 

- 0.0000 

- 0.0000 

1.0000 

1 

1 

2 

1 

- 0 . 000099999.0000 

1.0000 

-0 

5 

2 

1 

- 0.0000 

- 0.0000 

1.0000 

1 

1 

2 

2 

- 0 . 000099999.0000 

1.0000 

-0 

5 

2 

2 

- 0.0000 

- 0.0000 

1.0000 

1 

1 

2 

3 

- 0 . 000099999.0000 

1.0000 

-0 

5 

2 

3 

- 0.0000 

-o.oooc 

1.0000 

1 

2 

2 

1 

- 0 . 000099999.0000 

1.0000 

-0 

6 

2 

l 

- 0.0000 

- 0.0000 

1.0000 

1 

2 

2 

2 

- 0 . 000099999.0000 

1 .0000 

-0 

6 

2 

2 

- 0.0000 

-o.oooc 

1.0000 

1 

2 

2 

3 

- 0 . 000099999.0000 

1.0000 

-0 

6 

2 

3 

- 0.0000 

- 0.0000 

1.0000 

1 

1 

3 

1 

- 0 . 000099999.0000 

1.0000 

-0 

5 

3 

1 

- 0.0000 

-o.oooa 

1.0000 

1 

1 

3 

3 

- 0 . 000099999.0000 

1.0000 

-0 

5 

3 

3 

- 0.0000 

- 0.0000 

1.0000 

1 

2 

3 

3 

- 0 . 000099999.0000 

1.0000 

-0 

6 

3 

3 

- 0.0000 

- 0.0000 

1.0000 

-I 

-0 

-0 

-0 

- 0.0000 

-o.oooc 

- 0.0000 


CJ1 

CO 
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03 

O 


JKM = 16, NPM = 30 


THE fIRST COLUMN IS TIME. THE NEXT 3 COLUMNS ARE CONTROL VARIABLES. 

0.00 0. 0. 1. OOOOE+OO 8.9372E-05 -4. 1952E-05 

.05 0. 0. 1. OOOOE+OO 1.6384E-04 -8.0879E-0S 

.10 0. 0. 1. OOOOE+OO 2.2256E-04 -1.1661E-04 

.15 0. 0. 1. OOOOE+OO 2. 6481E-04 -1.4900E-04 

.20 0. 0. 1. OOOCE+OO 2.9J05E-04 -1.7794E-04 

• 25 0. 0. 1. OOOOE+OO 2. 97866-04 -2.0334E-04 

.30 0. 0. 1 • OOOOE+OO 2.8800E-04 -2.2513E-04 

■ 35 0. 0. 1. GOOOE+OO 2.6036E-04 -2.4331E-04 

.40 0. 0. 1. OOOOE+OO 2. 1502E-04 -2.5785E-04 

•45 0. 0. 1. OOOOE+OO 1.5216E-04 -2.6880E-04 

.50 0. 0. 1. OOOOE+OO 7.2150E-05 -2.7622E-04 

•55 0. 0. 1. C000£+00 -2.4530E-05 -2.8018E-04 

•60 0. 0. 1. OOOOE+OO -1.3725E-04 -2.8080E-04 

•65 0. 0. 1. OOOOE+OO -2.6526E-04 -2.7821E-04 

•70 0. 0. 1. OOOOE+OO -4.0768E-04 -2.7258E-04 

•75 0. U. 1. OOOOE+OO -5.6353E-04 -2.6408E-04 

-80 0. 5.0000E-01 1. OOOOE+OO 5.0921E-02 -4.7040E-02 

.85 0. 1. OOOOE+OO 1. OOOOE+OO 2.0320E-01 -1.8650E-01 

.90 4.9000E+00 1. OOOOE+OO 1. OOOOE+OO 1.3997E+00 -5.0905E-01 

.95 9. 8000E+00 1. OOOOE+OO 1. OOOOE+OO 4.5598E+00 -1.1024E+0Q 

1.00 9 . 5500E+00 1. OOOOE+OO 1. OOOOE+OO 8.5798E+00 -1.8129E+00 

1.05 S.3000E+00 l. OOOOE+OO 1. OOOOE+OO 1.2370E+01 -2.4881E+00 

1.10 2.5000E+00 1. OOOOE+OO 1. OOOOE+OO 1.4568E+01 -2.9374E+00 

1-15 -4.3000E+00 1. OOOOE+OO 1. OOOOE+OO 1.3844E+01 -2.9745E+00 

1.20 -4.3333E+00 1. OOOOE+OO 1. OOOOE+OO 1.1624E+01 -2.7959E+00 

1.25 -4.3667E+00 1. OOOOE+OO 1. OOOOE+OO 9.3161E+00 -2.5974E+00 

1.30 -4.4000E+00 1. OOOOE+OO 1. OOOOE+OO 6.9329E+00 -2.3815E+0Q 

1.35 -4.4333E+00 1. OOOOE+OO 1. OOOOE+OO 4..V878E +00 -2. 1504E + 00 

1.40 -4.4667E+00 1. OOOOE+OO 1. OOOOE+OO 1.9945E+00 -1.9070E+00 

1.45 -4.5000E+00 1. OOOOE+OO 1. OOOOE+OO -5.3294E-01 -1.6538E+00 

1.50 -2.2000E+00 1. OOOOE+OO l.OOOOt+OO -2.6036E+00 -1.4598E+00 

1.55 1.0000E-01 1. OOOOE+OO 1. OOOOE+OO -3.7379E+00 -1.3924S+00 

1.60 1.0000E-01 8.0833E-01 1. OOOOE+OO -4.4332E+00 -1.3681E+03 

1.65 l.OOOOE-Ol 6 • 1667E-0 1 1. OOOOE+OO -5.1803E+00 -1.3042E+00 

1.70 l • OOOOE-O 1 4.2500E-01 1. OOOOE+OO -5.9740E+00 -1.2017E+0Q 

1.75 1.0000E-01 2.3333E-01 1. OOOOE+OO -6.8080E+00 -1.0624E+00 

1.80 1.0000E-01 4. 1667E-02 1. OOOOE+OO -7.6749E+00 -8.8766E-01 

1.85 1 • OOOOE-O 1 -1.5000E-01 1. OOOOE+OO -8.5662E+00 -6.7951E-01 

1.90 l.OOOOE-Ol -3 .4 167E-0 1 1. OOOOE+OO -9.4728E+00 -4.3995E-01 

1.95 l.OOOOE-Ol -5.3333E-01 1. OOOOE+OO -1.0385E+01 -1.7112E-01 

2.00 l.OOOOE-Ol -7.2500E-01 1. OOOOE+OO -1.1291E+01 1.2468E-01 

2.05 l.OOOOE-Ol -9 • 1667E-01 1. OOOOE+OO -1.2182E+01 4.4504E-O1 

2.10 l.OOOOE-Ol -1.10e3E+00 1. OOOOE+OO - 1 . 3044E + 01 7.8748E-01 

2.15 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.3867E+01 1.1495E+00 

2.20 l.OOOOE-Ol -1.3000E+00 l. OOOOE+OO -1.4618E+01 1.5105E+00 

2.25 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.5266E+01 1.8503E+00 

2.30 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.5802E+01 2.1670E+00 

2.35 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.6217E+01 2.4589E+00 

2.40 l.OOOOE-Ol -1.3000E+00 l. OOOOE+OO -1.6503E+01 2.7245E+00 

2.45 1. OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.6655E+01 2.9625E+00 

2.50 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.6667E+01 3.1720E+00 

2.55 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.6535E+01 3.3521E+00 

2.60 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.6259E+01 3.5024E+00 

2.65 l.OOOOE-Ol -1.30G0E+Q0 1. OOOOE+OO -1.5835E+01 3.6226E+00 

2.70 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.5266E+01 3.7126E+00 

2.75 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.4551E+01 3.7727E+00 

2.80 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.3695E+01 3.8034E+00 

2.85 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.2701E+01 3.8052E+00 

2.90 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.1573E+01 3.7792E+00 

2.95 l.OOOOE-Ol -1.3000E+00 1. OOOOE+OO -1.0319E+01 3.7264E+00 


THE NEXT 7 COLUMNS ARE RESPONSE VARIABLES. 

2.2927E-05 2.2944E-06 1.6415E-03 -8.0941E-04 4.7933E-04 

4.7357E-05 8.6888E-06 1.3345E-03 -7.4711E-04 5.1709E-04 

7.4556E-05 1 . 841 6E-05 1.0118E-03 -6.8167E-04 5.5005E-04 

1.0278E-04 3 • 0670E-05 6.7648E-04 -6.1363E-04 5.7B03E-04 

1.3227E-04 4.4614E-05 3.3153E-04 -5.4357E-04 6.0087E-04 

1.6278E-04 5.9385E-05 -1.9951E-05 -4.7208E-04 6.1847E-04 

1.9403E-04 7.4105E-05 -3.7489E-04 -3.9974E-04 6.3076E-04 

2.2577E-04 8.7888E-05 -7.3023E-0<> -3.2714E-04 6.3773E-04 

2.5772E-04 9.9846E-05 -1.0829E-03 -2.5484E-04 6.3940E-04 

2.3962E-04 1.0910E-04 -1.4301E-03 -1.8342E-04 6.3581E-04 

3.2 12 IE-04 1. 1478E-04 -1.7688E-03 -1.1342E-04 6.2708E-04 

3.5224E-04 1.1604E-04 -2.0963E-03 -4.5389E-05 6.1334E-04 

3.8247E-04 1.1206E-04 -2.4100E-03 2.0174E-05 5.9477E-04 

4.1 164E-04 1 . 0206E-04 -2.7074E-03 8.2782E-05 5.7158E-04 

4.3955E-04 8.5290E-05 -2.9861E-03 1.4198E-04 5.4401E-04 

4.6598E-04 6.1064E-05 -3.2442E-03 1.9735E-04 5.1234E-04 

1.4466E-03 9. 9991E-04 2.0536E+00 -1.8679E+00 5.1327E-02 

7.4774E-03 7.0519E-03 4.0194E+00 -3.7056E+00 2.0201E-01 

2.6685E-02 4.0948E-02 4.3730E+01 -9.1812E+00 6.2049E-01 

7.7664E-02 1.8395E-01 8.2525E+01 -1.4529E+01 1.4713E+00 

1.77 39E-0 1 5 . 1324E-0 1 7.8187E + 01 -1.3875E + 01 2.5101E+00 

3.2771E-01 1 • 0379E+00 7.3338E+01 -1.3118E+01 3.4939E+00 

5.2013E-01 1 • 7207E +00 1.4618E + 01 -4.8520E + 00 4.1188E+00 

7.2736E-01 2.4401E+00 -4.3445E+01 3.3565E+00 4.0874E+00 

9.2 25 IE-01 3. 0772E+00 -4.5324E + 01 3.7779E + 00 3.7134E + 00 

1.0982E+00 3 • 6011E+00 -4.6956E+01 4.1528E+00 3.3093E+00 

1.2530E+00 4. 0076E+00 -4.8328E+01 4.4786E+00 2.B788E+00 

1.3857E+00 4.2933E+00 -4.9431E+01 4.7533E+00 2. 4256E+00 

1.+952E+00 4.4555E+00 -5.0255E+01 4.9751E+00 1.9535E+00 

1.5808E+00 4.4922E+00 -5.0796E+01 5.1427E+00 1.4665E+00 

1.6437E+00 4.4108E+00 -3.2028E+01 2.6141E+00 1.0744E+00 

1.6920E+00 4. 2494E+00 -1.3358E+01 8.1551E-02 8.8490E-01 

1.7337E+00 4. 0453E+00 -1.4440E+01 8.8655E-01 7.7370E-01 

1.7686E+00 3 . 805 1E+00 -1.5428E+01 1.6687E+00 6.1776E-01 

1.7946E+00 3. 5Z64E+00 -1 .6301E+01 2.4236E+00 4.1926E-01 

1.8098E+00 3.2069E+00 -1.7035E+01 3.1470E+00 1.8068E-01 

1.8120E+00 2. 8450E+00 -1.7612E+01 3.8350E+00 -9.5236E-02 

1.7996E+00 2.4390E+00 -1.8012E+01 4.4843E+00 -4.0548E-01 

1.7709E+00 1 . 9B80E+00 -1.8220E+01 5.0916E+00 -7.4683E-01 

1.7244E+00 1.4915E+00 -1.8223E+01 5.6542E+00 -1.1159E+00 

1.6588E+00 9.4952E-01 -1.8008E+01 6.1699E+00 -1.5090E+00 

1.5730E+00 3 . 6255E-0 1 -1.7566E+01 6.6366E+00 -1.9226E+00 

1.4662E+00 -2.6828E-01 -1.6889E+01 7.0530E+00 -2.3527E+00 

1 .33 75E + 00 -9.4128E-01 -1.5973E + 01 7.4180E + 00 -2.7955E+00 

1.1868E+00 -1.6538E+00 -1.4027E+01 7.0149E+00 -3.2276E+00 

1 .0 153E+00 -2.4013E+00 -1.1873E+01 6.5712E+00 -3.6262E+00 

8.2474E-01 -3 • 1785E+00 -9.5344E+00 6.0912E+00 -3.9887E+00 

6 . 1704E-01 -3.9795E+00 -7.0347E+00 5.5792E+00 -4.3128E+00 

3.9413E-01 -4.7980E+00 -4.3986E+00 5.0399E+00 -4.5965E+00 

1.5809E-01 -5. 6275E+00 -1.6513E+00 4.4780E+00 -4.8382E+00 

-8.8965E-02 -6.4612E+00 1.1814E+00 3.8983E+00 -5.0366E+00 

-3.4484E-01 -7. 2918E+00 4.0736E+00 3.3058E+00 -5.1908E+00 

-6.0730E-01 -8.1123E+00 6.9992E+00 2.7053E+00 -5.3002E+00 

-8.741 1E-0 1 -8. 9152E+00 9.9323E+00 2.1019E+00 -5.3646E+00 

-1.1430E+00 -9 • 6934E+00 1.2847E+01 1.5004E+00 -5.3840E+00 

- 1.4118E+00 -1.0439E+01 1.5719E+01 9.0566E-01 -5.3590E+00 

-1.6782E+00 — 1 . 11 46 E +01 1.8524E+01 3.2229E-01 -5.2904E+00 

- 1 . 9401E + 00 -1 . 1807E+0 1 2.1237E + 01 -2.4515E-01 -5.1794E+00 

-2.1954E+00 - 1 . 2414E+0 1 2.3837E+01 -7.9231E-01 -5.0272E+00 

-2.4422E+00 rl.2962E+0l 2.5303E+01 -1.3151E+00 -4.8359E+00 
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3.00 1. 00006-01 -1.30006*00 

3.05 1.00006-01 -1.30006+00 

3.10 1.00006-01 -1.30006+00 

3.15 1.00006-01 -1.30006+00 

3.20 l. 00006-01 -1.30006+00 

3.25 1.00006-01 -1.30006+00 

3.30 1.00006-01 -1.30006+00 

3.35 1.00006-01 -1.30006+00 

3.40 l. 00006-01 -1.30006+00 

3.45 1.00006-01 -1.30006+00 

3.50 1.00006-01 -1.30006+00 

3.55 1.00006-01 -1.30006+00 

3.60 1.00006-01 -1.30006+00 

3.65 1.00006-01 -1.30006+00 

3.70 1.00006-01 -1.30006+00 

3.75 1.00006-01 -1.30006+00 

3.80 1.00006-01 -1.30006+00 

3.85 1.00006-01 -1.30006+00 

3.90 1.00006-01 -1.30006+00 

3.55 1.00006-01 -1.30006+00 

4.00 1.00006-01 -1.30006+00 

4.05 1. 00006-01 -1.30006+00 

4.10 1.00006-01 -1.30006+00 

4.15 1.00006-01 -1.30006+00 

4.20 1.00006-01 -1.30006+00 

4.25 1.0000E-01 -1.30006+00 

4.30 1.00006-01 -1.30006+00 

4.35 1.0000E-01 -1.30006+00 

4.40 1.00006-01 -1.30006+00 

4.45 l. 00006-01 -1.3000E+00 

4.50 l. OOOOE-O 1 -1.30006+00 

4.55 1.0000E-01 -1.30006+00 

4.60 1.00006-01 -1. 30006+00 

4.65 1.00006-01 -1.30006+00 

4.70 1.0000E-01 -1.30006+00 

4.75 1.00006-01 -1.30006+00 

4.80 l. 00006-01 -1.30006+00 

4.85 1.00006-01 -1.30006+00 

4.90 1.00006-01 -1.30006+00 

4.95 1.00006-01 -1.30006+00 

5.00 -2.35006+00 -1.3000E+00 

5.05 -4.80006+00 -1.30C0E+00 

5.10 -7.25006+00 -1.3000E+00 

5.15 -9.7000E+00 -1.30006+00 

5.20 -9.5000E+00 -9.0000E-01 

5.25 -9.30006+00 -5. 00006-01 

5.30 -7.90006+00 -5.00006-01 

5.35 -6.5000E+00 -5.00006-01 

5.40 -5.1000E+00 -5.00006-01 

5.45 -3.70006+00 -5.00006-01 

5.50 -2.30006+00 -5.0000E-01 

5.55 -9.00006-01 -5.00006-01 

5.60 5.00006-01 -5.00006-01 

5.65 1.9000E+00 -5.0000E-01 

5.70 3.30006+00 -5.00006-01 

5.75 4.7000E+00 -5.0000E-01 

5.80 6.10006+00 -5.00006-01 

5.85 7.5000E+00 -5.00006-01 

5.90 8.90006+00 -5.00006-01 

5.95 1.03006+01 -5.00006-01 

6.00 1.04006+01 -5.00006-01 

6.05 1.05006+01 -5.00006-01 


1.00006+00 -8.9454E+00 3.64826+00 -2.6784E+00 
1.00006+00 -7.46056+00 3.5460E+00 -2.9023E+00 
1.00006+00 -5.87326+00 3.42156+00 -3.1122E+00 
1.00006+00 -4.19346+00 3.27666+00 -3.30666+00 
1.00006+00 -2.43176+00 3.11336+00 -3.48406+00 
1.00006+00 -5.9913E-01 2.93366+00 -3.5432E+00 
1.00006+00 1.2926E+00 2.73976+00 -3.7830E+00 
1. 00006+00 3.23146+00 2.5339E+00 -3.9026E+00 
l.OOOOE+OO 5.20496+00 2.3186E+00 -4.0012E+00 
1. OOOOE+OO 7.2006E+00 2.0962E+00 -4.0782E+00 
1.00006+00 9.20596+00 1.86926+00 -4.1332E+00 
1.00006+00 1.12086+01 1.63986+00 -4.16596+00 
1.00006+00 1 .31946+01 1.41066+00 -4.17636+00 
1.00006+00 1.51526+01 1.18416+00 -4.16456+00 
1.00006+00 1.70706+01 9.62406-01 -4.1309E+00 
1.00006+00 1.89366+01 7.47966-01 -4.0758E+00 
1.00006+00 2.07396+01 5.42936-01 -3.9999E+00 
l.OOOOE+OO 2.24686+01 3.49436-01 -3.90396+00 
1.00006+00 2.41136+01 1.69436-01 -3.7889E+00 
l.OOOOE+OO 2.56656+01 4.81286-03 -3.6557E+00 
1. 00006+00 2.71156+01 -1.42716-01 -3.5058E+00 
l.OOOOE+OO 2.84566+01 -2.71586-01 -3.3402E+00 
1.00006+00 2.96806+01 -3.8Q38E-01 -3.1605E+00 
l.OOOOE+OO 3.07816+01 -4.67B8E-01 -2.56816+00 
1. GOQOE+OO 3.17566+01 -5.33026-01 -2.7646E+00 
1.00006+00 3.25996+01 -5.7495E-01 -2.55176+00 
1.00006+00 3 . 3307E+01 -5.92986-01 -2.3311E+00 
1.00006*00 3.3880E+01 -5.86636-01 -2.1045E+00 
l.OOOOE+OO 3.43156+01 -5.5560E-01 -1.8736E+00 
l.OOOOE+OO 3.46136+01 -4.99816-01 -1.6403E+00 
1.00006+00 3.47746+C1 -4.19336-01 -1.4063E+00 
1.00006+00 3.48016+01 -3.14476-01 -1.17346+00 
l.OOOOE+OO 3.46986+01 -1.85696-01 -9.4323E-01 
l.OOOOE+OO 3. 4467E+0 l -3.36446-02 -7.1755E-01 
l.OOOOE+OO 3.41136+01 1.40866-01 -4.9797E-01 
1.00006+00 3.36436+01 3.3683E-01 -2.8604E-01 
l.OOOOE+OO 3.30626+01 5.53136-01 -8.32346-02 
l.OOOOE+OO 3.23786+01 7.8849E-01 1.09056-01 
1.00006+00 3.15986+01 1.04156+00 2.8953E-01 
l.OOOOE+OO 3.07326+01 1.3106E+00 4.57026-01 
1.00006+00 2.92866+01 1.66376+00 6.C836E-01 
l.OOOOE+OO 2.67836+01 2.16706+00 7.33646-01 
l.OOOOE+OO 2.32556+01 2.81536+00 8.2115E-01 
l.OOOOE+OO 1.8740E+01 3.60256+00 8.59596-01 
l.OOOOE+OO 1.38686+01 4.4086E+00 8.4121E-01 
l.OOOOE+OO 9.25386+00 5.11496+00 7.68766-01 
1.00006+00 5.12806+00 5.72056+00 6.4863E-01 
l.OOOOE+OO 1.71136+00 6.22486+00 4.8833E-01 
l.OOOOE+OO -9.85756-01 6.6253E+00 2.95856-01 
1.00006+00 -2.95746 +00 6.92006+00 7.93486-02 
1.00006+00 -4.2025E+00 7.1078E+00 -1.52926-01 
1.00006+00 -4.72496+00 7.18876+00 -3.92656-01 
l.OOOOE+OO -4.53266+00 7.1634E+00 -6.3158E-01 
l.OOOOE+OO -3.6387E+00 7.0337E+00 -8.6155E-01 
l.OOOOE+OO -2.0603E+00 6.80206+00 -1.07456+00 
1.00006+00 1.81226-01 6.47176+00 -1.26286+00 
1.00006+00 3.06026+00 6.04706+00 -1.41B7E+00 
l.OOOOE+OO 6.54756+00 5.53286+00 -1.53516+00 
1.00006+00 1.06106+01 4.93466+00 -1.6053E+00 
l.OOOOE+OO 1.52126+01 4.25886+00 -1.6227E+00 
l.OOOOE+OO 2.00496+01 3.5489E+00 -1.5B27E+00 
1.00006+00 2.48196+01 2.84856+00 -1.48466+00 


1.34446+01 2.86156+01 -1.8096E+00 -4.5072E+00 

-1 . 3855E+01 3 .0754E+01 -2.2723E+00 -4.34376+00 

-1.41886+01 3.2705E+01 -2.6998E+00 -4.0477E+00 

-1.4440E+01 3 .445 1E+01 -3.0892E+00 -3.7220E+00 

- 1. 4606E+01 3.59816+01 -3.4379E+00 -3.3696E+00 

-1.46826+01 3.72826+01 -3.7435E+00 -2.9937E+00 

-1.4665E+01 3.8346E+01 -4.0042E+00 -2.59746+00 

-1.45526+01 3.9165E+01 -4.2183E+00 -2.18416+00 

-1.4341E+01 3.9735E+01 -4.3849E+00 -1.75746+00 

-1.40316+01 4.0052E+01 -4.5030E+00 -1.3206E+00 

-1.36216+01 4.01 15E+01 -4.57246+00 -8.7740E-01 

-1.3111E+01 3.99256+01 -4.5931E+00 -4.3131E-01 

-1.25016+01 3.9485E+01 -4.5654E+00 1.4138E-02 

-1. 1792E+01 3.8801E+01 -4.49026+00 4.5547E-01 

-1.09866+01 3.7878E+01 -4.3685E+00 8.8929E-01 

-1 . 0086E+01 3.6725E+01 -4.2019E+00 1.3123E+00 

-9.0935E+00 3.5353E+01 -3.99226+00 1.7213E+00 

-8.0130E+00 3.3773E+01 -3.74146+00 2.1134E+00 

-6 . 84816+00 3.2000E+01 -3.4521E+00 2.4856E+00 

-5.6033E+00 3.0047E+01 -3.12706+00 2.83536+00 

-4.28336+00 2.7930E+01 -2.76896+00 3.16016+00 

-2. 8936E+00 2.5667E+01 -2.3810E+00 3.4577E+00 

-1.43976+00 2.3276E+01 -1.96686+00 3.72626+00 

7. 2358E-02 2.0776E+01 -1.52976+00 3.96386+00 

1.63636+00 1.8185E+01 -1.0733E+00 4.16926+00 

3. 2458E+00 1.5525E+01 -6.0151E-01 4.34106+00 

4. 8940E+00 1.28 146+01 -1.1806E-01 4.4784E+00 

6. 5742E+00 1.0074E+01 3.7315E-01 4.5807E+00 

8. 2797E+00 7.3254E+00 8.6823E-01 4.6476E+00 

1.00036+01 4.58756*00 1.3633E+00 4.6790E+00 

1 . 1739E+01 1 • 8806E+00 1.8546E+00 4.6752E+00 

I . 3479E+0 l -7.7582E-01 2.33836+00 4.6365E+00 

1 . 5217E+01 -3.3630E+00 2.81086+00 4.5S37E+00 

1.6946E+01 -5.8629E+00 3.2685E+00 4.4579E+00 

1. 8661E+01 - 8. 25 84E + 00 3.7083E+00 4.3202E+00 

2.03566+01 -1 .05 33E + 01 4.1268E+00 4.1521E+00 

2.2024E+01 -1.2673E+01 4.52116+00 3.95536+00 

2.3660E+01 -1.46 63E + 01 4.8885E+00 3.73176 + 00 

2.5260E+01 - 1 . 6491E+01 5.2265E+00 3.4834E+00 

2.6818E+01 -1.8146E+01 5.532BE+00 3.2126E+00 

2 . 8322E+01 -3.9592E+01 8.57896+00 2.8105E+00 

2 • 9727E+01 -6.0442E+01 1.1536E+01 2.1709E+00 

3 . 0981 E +01 -8.0572E+01 1.4379E + 01 1.3010E+00 

3.2034E+01 -9.9863E+01 1.7086E+01 2.10116-01 

3. 2849E+01 -9.4956E+01 1.51406+01 -9.30396-01 

3.34266+01 -8.95106+01 1.3098E+01 -1.9513E+00 

3.37836+01 -7.54646+01 1.1110E+01 -2.8332E+00 

3.3952E+01 -6.11616*01 9. 05636+00 -3.5574E+00 

3. 3967E+01 -4.6695E+01 6.9558E+00 -4.1199E+00 

3.3866E+01 -3.21 60E+01 4.8261E+00 -4.5179E+00 

3. 3685E+01 -1.7652E+01 2.6853E+00 -4.7503E+00 

3.3459E+01 -3.2624E+00 5.5161E-01 -4.8168E+00 

3.3225E+01 1.0915E+01 -1.5574E+00 -4.7187E+00 

3. 3019E+0 1 2.4792E+01 -3.6242E+00 -4.4587E+00 

3 • 2S74E+01 3.8281E+01 -5.6322E+00 -4.04036+00 

3.28256+01 5.1300E+01 -7.56496+00 -3.46866+00 

3. 2904E+01 6.3770E+01 -9.4072E+00 -2.7496E+00 

3. 3142E+01 7.5618E+01 -1.11446+01 -1.8904E+00 

3.3570E+01 6.6775E+01 -1.2763E+01 -8.9924E-01 

3 . 4214E+0 1 9.7179E+01 -1.4250E+01 2.1490E-01 

3. 5095E+01 9.6176E+01 -1.4123E+01 1.38326+00 

3.6218E+01 9.4529E+01 -1.3872E+01 2.5371E+00 
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6.10 1.0600E+01 -5.0000E-01 

6.15 1.0700E+01 -5.0000E-01 

6.20 9. 5000E+00 -5.0000E-01 

6.25 8.3000E+00 -5.0000E-01 

6.30 7. 1000E+00 -5.0000E-01 

6.35 5 . 9000E *-00 -5.0000E-01 

6.40 4.7000E+00 -5.0000E-01 

6.45 3. 5000E+00 -5.0000E-01 

6.50 2. 3000E+00 -5.0000E-01 

6.55 1.1000E+00 -5 . OOOOE-O 1 

6.60 -l.OOOOE-Ol -5.0000E-01 

6.65 -1.3000E+00 -5.0000E-01 

6.70 -2.5000E+00 -5.0000E-01 

6.75 -3.7000E+00 -5.0000E-01 

6.80 -3.7000E+00 9.0000E-01 

6.85 -3.7000E+00 2.3000E+00 

6.90 -3. 7000E+00 3.7000E+00 

6.95 -3. 7000E+00 5.1000E+00 

7.00 -3.7000E+00 5.0900E+00 

7.05 -3.7000E+00 5.0800E+00 

7.10 -3.7000E+00 5.0700E+00 

7.15 -3. 7000E+00 5.0600E+00 

7.20 -3.7000E+00 5.0500E+00 

7.25 -3.7000E+00 5.0400E+00 

7.30 -3.7000E+00 5.0300E+00 

7.35 -3.7000E+00 5. 0200E+00 

7.40 -3.7000E+00 5.0100E+00 

7.45 — 3.7000E+00 5. 00006+00 

7.50 -3.7000E + 00 4.9900E+00 

7.55 -3.700QE+00 4.9800E+00 

7.60 -3.7000E+00 4.9700E+00 

7.65 -3.7000E+00 4.9600E+00 

7.70 -3.7000E+00 4.9500E + 00 

7.75 -3.7000E+00 4.9400E+00 

7.80 — 3 . 70Q0E+Q0 4.9300E+00 

7.85 -3.7000E+00 4.9200E+00 

7.90 -3. 7000E+00 4.9100E+00 

7.95 -3.7000E+00 4.9000E + 00 

8.00 -3.7000E+00 3.4000E+00 

8.05 -3.7Q00E+00 1.9000E+00 

8.10 -3.7Q00E+00 4.0000E-01 

8.15 -3.7000E+00 -1.1000E+00 

8.20 -3.4429E+00 -1.0357E+00 

8.25 -3. 18576+00 -9.7143E-01 

8.30 -2.9286E+00 -9.0714E-01 

8.35 -2.6714E+00 -8.4286E-01 

8.40 -2.4143E+00 -7.7657E-01 

8.45 -2.1571E+00 -7.1429E-01 

8.50 -1 .90006+00 -6.5000E-01 

8.55 -1.6429E+00 -5.8571E-01 

8.60 -1.3857E+00 -5.2143E-01 

8.65 - l. 1286E+00 -4.5714E-01 

8.70 -8.7143E-01 -3.9286E-01 

8.75 -6.1429E-01 -3.2857E-01 

8.80 -3.5714E-01 -2.6429E-01 

8.85 — l . OOOOE-O 1 — 2.000QE-01 

8.90 -9.5455E-02 -1.9091E-01 

8.95 -9.0909E-02 -1.8182E-01 

9.00 -8.6364E-02 -1.7273E-01 

9.05 -8.1818E-02 -1.6364E-01 
9. 10 -7.7273E-02 -1.5455E-01 

9.15 -7.2727E-02 -1.4545E-01 

9.20 -6.8182E-02 -1.3636E-01 


1. OOOOE+OQ 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
1 . OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 
l.OOOOE+OO 


2.9491E+01 
3.4035E+01 
3 . 8155E+0 1 
4. 1565E+0 1 
4.4250E+01 
4.6203E+01 
4 . 7419E+01 
4.7900E+01 
4. 7652E+01 
4.6687E+01 
4.5018E+01 
4.2668E+01 
3.9661E+01 
3.6025E+01 
3.2163E+C1 
2 . 8547E+0 1 
2.5130E+01 
2. 1935E+01 
1.8815E+01 
1.5619E+01 
1.2349E+01 
9 • 0073E+00 
5.5957E+00 
2 . 1187E+00 
-1.4185E+00 
-5.0096E+00 
-8.6474E+00 
-1.2324E+01 
-1.6030E+01 
-1.9756E+01 
-2.349 IE + 01 
-2.7226E+01 
-3.0949E+01 
-3.4647E+01 
-3.8311E+01 
-4. 1927E+0 1 
-4. 5485 E +01 
-4. 8972E+01 
-5.2531E+01 
-5.6298E+01 
-6.0244E+01 
-6.4332E+01 
-6.8308E+01 
-7.1916E+01 
-7.5124E+01 
-7.7901E+01 
-8.0222E+01 
-8.2066E+01 
-8. 3416E+01 
-8.4259E+01 
-8.4586E+01 
-8.4394E+01 
-8.3683E+01 
—8. 2459E+01 
-8. 0730E+01 
-7.8509E+01 
-7.5870E+01 
-7.2890E+01 
-6.9591E+01 
-6.5995E+01 
-6.2128E+01 
-5.8019E+01 
-5.36966+01 


2. 1637E+00 
1.5004E+00 
9.01 15E-0 1 
4. 0744E-0 1 
2.2651E-02 
-2.5080E-01 
-4. 1144E-01 
-4.5876E-01 
-3.9317E-01 
-2.1604E-01 
7. 0346E-02 
4.6284E-01 
9 . 5745E-0 1 
1 -5493E+00 
2.0679E+00 
2.3451E+00 
2.3805E+00 
2.1752E+00 
1 . 8635E+00 
1.5784E+00 
1.3 197E+00 
1.0872E+00 
8 .8004E-01 
6 . 9754E-0 1 
5.3856E-01 
4. 0 184E-01 
2 • 859 IE-01 
1.89 15E-0 1 
1.0975E-01 
4.5794E-02 
-4.7846E-03 
-4.4138E-02 
-7.4509E-02 
-9. 82016-02 
-1 . 1756E-0 1 
-1 . 3497E-0 1 
-1.5281E-01 
-1 .7343E-01 
-5. 8748E-02 
3 . 2276E-0 1 
9 • 6659E-0 1 
1 .8628E+00 
2. 8476E+00 
3.75876+00 
4.5 896E+00 
5.3344E+00 
5 . 9 879E + 00 
6.5459E+00 
7.0047E+00 
7 . 3618E+00 
7.6152E+00 
7.7639E+00 
7.8077E+00 
7. 7472E + 00 
7 « 5839E+00 
7. 3199E+Q0 
6 . 9707E+00 
6.5519E+00 
6.Q675E+00 
5.5216E+00 
4.9191E+00 
4.2652E+00 
3.5653E+00 


- 1.3294E+00 
-1.1183E+00 
-8.5440E-01 
-5.445 IE-01 
-1.9700E-01 

1.7961E-01 
5.7664E-01 
9 .85386-01 
1.3971E + 00 
1.8032E+00 
2.1950E+00 
2.5645E+00 
2.9034E+00 
3.2040E+00 
3.4629E+00 
3.6889E+00 
3.8946E+00 
4.0928E+00 
4.2928E+00 
4.4952E+00 
4.6977E+00 
4.8979E+00 
5 .0938E+00 
5.2832E+00 
5.4641E+00 
5.6346E+00 
5.7927E+00 
5.9369E+00 
6.0656E+00 
6. 1774E+00 
5.2710E+00 
6.3454E+00 
6.3997E+00 
6.4331E+00 
6.4452E+00 
6.4354E+00 
6.4037E+00 
6.3500E+00 
6.2717E+00 
6.1568E+00 
5. 9914E+00 
5.7619E+00 
5.4590E+00 
5 .0848E+-00 
4.6452E+00 
4.1462E+00 
3.5947E+00 
2.9974E+00 
2.36 16E+00 
1.6947E+00 
1.0041 E+00 
2.9759E-01 
-4.1732E-01 

- 1. L330E + 00 
-1 .84206+00 
-2. 53696+00 
-3.2110E+00 
-3.8588E+00 
-4.47536+00 
-5.0561E+00 
-5.5967E+00 
-6.0935E+00 
-6.5429E+00 


3. 7576E+01 
3 . 9165E+0 1 
4. 0972E+01 
4.2967E+01 
4. 51 15E+0 l 
4.7379E+01 
4.9722E+01 
5.2108E+01 
5.4499E+01 
5.6860E+01 
5. 9155E+01 
6. 1349E+01 
6. 3410E+01 
8.5304E+01 
6 . 7009E+0 l 
6. 8526E+01 
6.9B68E+0 1 
7. 1044E+01 
7.2063E+01 
7 • 2924E+01 
7. 3623E+01 
7. 4158E+01 
7.4523E+01 
7. 4716E+0 1 
7.4734E+01 
7 • 4573E+01 
7.4232E+01 
7. 3708E+01 
7 • 2999E+01 
7.2105E+01 
7.1023E+01 
6 • 9755E+01 
6.8301E+01 
6.6661E+01 
6.4837E+01 
6.2831E+01 
6. 0645E+01 
5 ■ 8283E+0 1 
5. 5746E+01 
5 . 3026E+0 1 
5. 0113E+01 
4. 6999E+01 
4.3681E+01 
4.0174E+01 
3 . 6497E+0 1 
3.2669E+01 
2.87 14E+01 
2. 4655E+01 
2.0516E+01 
1 . 6322E+01 
1 • 2099E+0 1 
7 . 8726E+00 
3.6687E+00 
-4.8689E-01 
-4.5686E+00 
-8.5514E+00 
-1.2412E+01 
-1.6133E+01 
-1. 9696E+01 
-2 . 3087E+01 
-2.6291E+01 
-2 • 9296E+01 
—3. 2089E+01 


9.22586+01 
8.9388E+01 
7. 53496+01 
6.0987E+01 
4.6402E+01 
3.1692E+01 
1.6959E+01 
2.3003E+00 
-1.21B5E+01 
-2.6403E+01 
-4.0259E+01 
-5.3666E + 01 
-6.6539E+01 
-7.8798E + 01 
-7.4825E+01 
-7.0550E+01 
-6.6119E+01 
-6. 1677E+01 
-6.3164E+01 
-6.4655E+01 
-6.6125E+01 
-6.7549E+01 
-6.8902E+01 
-7.0161E+01 
-7.1306E+01 
-7.23 15E+01 
-7.3170E + 01 
-7.3854E+01 
-7.4353E + 01 
-7.4652E+01 
-7.4740E+01 
-7 .46 10E+01 
-7.4252E+01 
-7.3663E+01 
-7.2839E+01 
-7.1780E+01 
-7.0486E + 01 
-6.8961E+01 
-7.3341E+01 
-7. 723 IE + 01 
-8.0479E+01 
-8.2941E+01 
-7 .59 56E + 01 
-6.8263E+01 
-5 .99 39E + 01 
-5.1062E+01 
-4.17 16E+01 
-3.19 85E+01 
-2.1957E+01 
-1.1719E+01 
-1.3592E+00 
9.0339E+00 
1.9374E+01 
2.9575E+01 
3.9555E+01 
4.9235E+01 
5.6250E+01 
6.2867E+01 
6.90366+01 
7 .47 13E + 01 
7.9855E+01 
8.44266+01 
8.8396E+01 


-1.3501E+01 
-1.3013E+01 
-1 . 0942E+01 
-8.7943E+00 
-6.5886E+00 
-4.3440E+00 
-2 • 0794E+00 
1.8606E-01 
2.4335E+00 
4.6445E+00 
6.8012E+00 
8. 8860E+00 
1.0883E+01 
1.2775E + 01 
7.9600E+00 
3. 1239E+00 
-1.7055E+00 
-5.5011E+00 
-5.9687E+00 
-5.4368E+00 
-4.9104E+00 
-4. 3942E+00 
-3.8929E+00 
-3.4109E+00 
-2.9524E+00 
-2.5213E+00 
-2.1212E+00 
-1.7554E+00 
-1.4269E+00 
-1.1383E+00 
-8.91 97E-01 
-6.8964E-01 
-5.3285E-01 
-4.2267E-01 
-3.5980E-01 
-3.4452E-01 - 
-3.7673E-01 - 
-4.5595E-01 - 
4.9857E+00 - 
1.02 92E + 01 - 
1. 5434E+01 - 
2. 0384E+01 - 
1.8982E+01 - 
1.7441E+01 - 
1.5775E+01 - 
1 • 3999E+01 - 
1.2128E+01 - 
1.0179E+01 - 
8.1667E+00 - 
6.1 093E+00 - 
4.0228E+00 - 
1.9241E+00 - 
-1.7034E-01 - 
-2.2441E+00 - 
-4.28 12E+00 - 
-6.2662E+00 - 
-7.6922E+00 - 
-9.0462E+00 - 
-1.0318E+01 - 
-1.1500E+01 - 
-1.2582E+01 - 
-1 . 3558E + 01 - 
-1.4421E+01 - 


3.6680E+00 
4.7673E+00 
5.7679E+00 
6 . 6045E+00 
7.2722E+00 
7.7677E+00 
8.0891E+00 
8.2358E+00 
8 .2086E+00 
8 .0097E+00 
7.6425E+00 
7.1117E+00 
6.4233E+00 
5 .5 844E+00 
4.7999E+00 
4.2695E+00 
3.9914E+00 
3.9620E+00 
4 .03216+00 
4.0562E+00 
4.0348E+00 
3 .9691E+00 
3. 8601E+00 
3.7096E+00 
3.5194E+00 
3.2917E+00 
3.0290E+00 
2.7338E+00 
2.4092E+00 
2.0581E+00 
1.6838E+00 
1.2897E+00 
8.7925E-01 
4.5601E-01 
2.3589E-02 
4.1438E-01 
8.5428E-01 
1.2925E+00 
1.8771E+00 
2.7507E+00 
3.8997E+00 
5.3088E+00 
6.7898E+00 
8. 1593E+00 
9.4081E+00 
1.0528E+01 
1.1513E+01 
1.2356E+01 
1.3053E+01 
1.3601E+01 
1.3997E+01 
1 .42416+01 
1.4332E+01 
1.4272E+01 
1.4064E+01 
1 .37116+01 
1.3235E+01 
1.2659E+01 
1.1987E+01 
1.1228E+01 
1.0387E+01 
9.4723E+00 
8.4924E+00 
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9.25 -6.3636E-02 -1.2727E-01 l.OOOOE+OO -4.9190E+01 2.8251E+00 -6.9418E+00 -3.4662E+01 

9.30 -5.9091E-02 -1.1818E-01 1.0000E+00 -4.4533E+01 2.0507E+00 -7.2876E+00 -3.7006E+01 

9.35 -5.45456-02 -1.0909E-01 1.0000E+00 -3.9758E+01 1.2483E+00 -7.5782E+00 -3.9113E+01 

9.40 -5.0000E-02 -1.0000E-01 l.OOOOE+OO -3.4898E+01 4.2431E-01 -7.8118E+00 -4.0980E+01 

9.45 -4.5455E-02 -9.0909E-02 l.OOOOE+OO -2.9987E+01 -4.1489E-01 -7.9871E+00 -4.2602E+0L 

9.50 -4.0909E-02 -8.1818E-02 l.OOOOE+OO -2.5060E+01 -1.2628E+00 -8.1035E+00 -4.3978E+01 

9.55 -3.6364E-02 -7.2727E-02 l.OOOOE+OO -2.0149E+01 -2.1129E+00 -8.1605E+00 -4.5109E+01 

9.60 -3.1818E-02 -6.3636E-02 l.OOOOE+OO -1.5289E+01 -2.9588E+00 -8.1584E+00 -4.5994E+01 

9.65 -2.7273E-02 -5.4545E-02 l.OOOOE+OO -1.0511E+01 -3.7942E+00 -8.0977E+00 -4.6639E+01 

9.70 -2.2727E-02 -4.5455E-02 l.OOOOE+OO -5.8485E+00 -4.6129E+00 -7.9795E+00 -4.7047E+01 

9.75 -1.8182E-02 -3.6364E-02 l.OOOOE+OO -1.3307E+00 -5.4090E+00 -7.8053E+00 -4.7226E+01 

9.80 -1.3636E-02 -2.7273E-02 l.OOOOE+OO 3.0129E+00 -6.1768E+00 -7.577 OE+OO -4.7183E+01 

9.85 -9.0909E-03 -1.8182E-02 l.OOOOE+OO 7.1548E+00 -6.9109E+00 -7.2970E+00 -4.6928E+01 

9.90 -4.5455E-03 -9.0909E-03 l.OOOOE+OO 1.1069E+01 -7.6062E+00 -6.9679E+00 -4.6472E+01 

9.95 -5.6932E-13 -1.1386E-12 l.OOOOE+OO 1.4733E+01 -8.2579E+00 -6.5927E+00 -4.5825E+01 

THE VALUE CF THE FIT ERROR IS 5.19084E+03 FOR THE ITERATION NUMBER 0 
MATRIX E£*D hAS 1 ROWS AND 7 COLUMNS 

1.0968E+03 1.0882E+01 1.2861E+01 1.2049E+03 2.7745E+03 

6.3254E+01 2.7582E+01 

MATRIX RMS HAS 1 ROWS AND 7 COLUMNS 

3.3118E+01 3.2989E+00 3.5863E+00 3.4712E+01 5.2674E+01 

7.9533E+00 5.2519E+00 

MATRIX GRAD HAS 1 ROWS AND 16 COLUMNS 
9.75376+02 1.0513E+01 -2.24916+02 2.1901E+02 -2.1106E+C2 

1.9261E+02 -*1.09776+01 3.1955E + C1 2.0065E+01 -2.7480E + 01 

4.0500E+00 2 .417 1E+01 9.39466+01 2.6661E-01 -2.5070E+00 

5.19C8E+03 

WOULD VCU BELIEVE 

1.5117555E-05 

MATRIX OC FAS 1 ROWS AND 15 COLUMNS 

-4.0576E-01 6. 27526-01 -5.5459E+C0 9.6505E+00 -5.3042E+0C 

— 1.2882E+00 6.1596E-02 -1.3353E-C1 1.6400E+00 -1.5716E+00 

-2.0338E+00 -8. 24336-02 3.3916E-02 8.2978E-03 -5.4547E-01 

THE VALUE OF THE FIT ERROR IS 5.35909E+02 FOR THE ITERATION NUMBER 1 
MATRIX EE*0 HAS 1 ROWS AND 7 COLUMNS 

8.1382E+01 8.06956+00 2.4167E+00 2.5188E+01 3.9272E+02 

1.9191E+01 6.9371E+00 

MATRIX RMS HAS 1 ROWS AND 7 COLUMNS 

9.02 12E+00 2 .84076+00 1.5546E+00 5.0188E+00 1.9817E+01 

4.3807E+00 2. 6338E+00 

MATRIX GRAD HAS 1 ROWS ANO 16 COLUMNS 
7.4488E+01 -3.6222E+00 -6.2248E-C2 -1.3093E+01 2.1587E+01 

1.56616+01 -7 .43646+02 3.9285E+01 5.5023E+01 2.9144E+01 

-7.7703E+00 6.2957E+01 7.9773E+01 2.3470E+01 1.1470E+01 

5.35916+02 

WOULD YOU BELIEVE 

1.59859926-04 

MATRIX DC HAS 1 ROWS ANO 15 COLUMNS 

5.77706-01 1.4666E+00 -3.8994E+00 -7.3725E-01 1.2303E+01 

-1.60136+00 -2.24276-01 -7.7426E-C1 -2.0399E-01 4.53316-02 

— 1.8603E+00 1 . 0177E+00 6.C8796-C2 -6.8057E-01 8.6466E-01 

THE VALUE OF THE FIT ERROR IS 2.67663E+02 FOR TH6 ITERATION NUMBER 2 
MATRIX EE*D HAS l ROWS ANO 7 COLUMNS 

3.6351E+01 3.9183E+00 3.7904E+C0 1.4265E+02 7.2200E+01 

5.10116+00 3.65216+00 

MATRIX RMS HAS 1 ROWS AND 7 COLUMNS 

6.02926+00 1.9795E+00 l. 94696+00 1.19446+01 8.4971E+00 

2.25866+00 1.91116+00 


9.17386+01 -1.5166E+01 -7.45566+00 
9.4432E+01 -1.5789E+01 -6.3706E+00 
9.6461E+01 -1.62856+01 -5.2467E+00 
9.7817E+01 -1.6654E+01 -4.09306+00 
9.8493E+01 -1.5893E+01 -2.9190E+00 
9. 8492E+01 -1.7002E+01 -1.7340E+00 
9.78 19E+01 -1.6981E+01 -5.4755E-01 
9.6485E+01 -1.6834E+01 6.31216-01 
9.4508E+01 -1.6561E+01 1.7931E+00 
9.1907E+01 -1.6168E+01 2.9294E+00 
8.8709E+01 -1.5658E+01 4.0315E+00 
8.4944E+01 -1.50366+01 5.0913E+00 
8.0646E+01 -1.4310E+01 6.1011E+00 
7.5854E+01 -1.3485E+01 7.0536E+00 
7.0607E+01 -1.2570E+01 7.9422E+00 


03 

00 
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MATRIX C-RAC HAS 1 ROWS AND 16 COLUMNS 

-7.9885E+Q1 1.63976+01 7. 75586+01 -3.68346+01 2.52436+01 

2.82306+01 1 . 8387E+03 3.3001E+C1 2.4093E+02 6.85726+01 

3.60C8E+01 1.42356+02 -3.3650E+02 -3.6445E+02 -7.45816+01 

2.67666+02 

WOULD YOU BELIEVE 

7.66613946-07 

MATRIX DC HAS 1 ROWS AND 15 COLUMNS 

-4.85736-01 -2.40426+00 -2.1477E+C0 -7.51566-01 -2.24746+00 
4.07166+00 1.7355E-01 9.04766-01 8.5963E-01 4.16266-01 

2.4093E-C1 -8. 69476-01 -1.8133E-02 3.18986-01 1.8371E-01 

THE VALUE OF THE FIT ERROR IS 1.11355E+02 FOR THE ITERATION NUMBER 3 
MATRIX EE*D HAS 1 ROWS ANO 7 COLUMNS 
1.51706+01 1.97056+00 1.2194E-01 7.87616+01 1.43276+01 

7.64156-01 2.40936-01 

MATRIX RMS HAS 1 ROWS AND 7 COLUMNS 

3.8949E+00 1.40376+00 3.4920E-C1 8.87476+00 3.7851E+00 

8.7416E-01 4.90846-01 

MATRIX GRAD HAS 1 ROWS AND 16 COLUMNS 
-3.1804E+C2 —8 .8816E+01 -5.8652E+C0 1.45046+01 -3.4808E+0C 

-4.82156+01 -1.64846+03 -5.5075E+C2 3.24476+01 -6.27636+01 

4.38896+01 -2.3004E+02 4.47476+02 2. 08536+01 4.2124E+01 

1.11366+02 

WOULD YOU BELIEVE 

3.11648006-04 

MATRIX DC HAS 1 ROWS ANO 15 COLUMNS 

3.78916-02 3.13806-01 1.0566E-C1 4.58326-02 -5.53116-01 

-8.04686-01 -1.32626-02 -1.83346-01 -1.32726-01 -2.92556-02 
-7.7024E-02 -3.4351E-02 5.01116-03 1.1390E-02 -4.3531E-01 

THE VALUE CF THE FIT ERROR IS 2.72967E+00 FOR THE ITERATION NUMBER 4 
MATRIX 6E*0 HAS 1 ROWS ANO 7 COLUMNS 

2.10436-01 7.64566-02 9.7859E-C4 2.3095E+00 1.02436-01 

2.7835E-Q2 2.04226-03 

MATRIX RMS HAS 1 ROWS ANO 7 COLUMNS 

4.5872E-01 2.76516-01 3.12826-02 1.51976+00 3.2005E-01 

1.66846-01 4.51916-02 

MATRIX GRAD HAS 1 ROWS ANO 16 COLUMNS 
-2.89906+01 -1.02376+01 4.7531E-C1 1.46536+00 7.73126-01 

-7.4411E+00 -1.60446+02 -5.74216+01 2.32246+00 -3.1222E-01 

7.22396+00 -3.71766+01 1.25C56+C0 2.00226-01 1.13486+01 

2.72S7E+00 

WOULD YOU BELIEVE 

8.6634005E-05 

MATRIX DC HAS 1 ROWS ANO 15 COLUMNS 

2.18016-03 5.82266-02 9.48426-02 1.48506-02 -2.5872E-02 

-2.7B43E— 01 3.3919E-03 -2.30066-02 5.66756-03 -2. 74526-04 

-3.04106-02 -4.4265E-02 1.2779E-C3 2.47706-02 -8.4584E-02 

THE VALUE OF THE FIT ERROR IS 4.410526-03 FOR THE ITERATION NUMBER 5 
MATRIX EE*D HAS 1 ROWS AND 7 COLUMNS 

6.7047E-04 2.2670E-04 9.21826-05 4.56416-04 2.68496-03 

2. 06986-04 7.28706-05 

MATRIX RMS HAS 1 ROWS AND 7 COLUMNS 

2.5893E-0 2 1.5057E-02 S.6012E-C3 2.13646-02 5.1816E-02 

1.43876-02 8.53646-03 

MATRIX GRAO HAS 1 ROWS ANO 16 COLUMNS 
2.29896-01 -9 . 2070E-02 -2.2273E-C1 1.9753E-01 -4.3745E-02 

-8.41 19E-02 -9.07506+00 -2.31976-01 1.36946-01 -5.80276-01 

-1.53266-01 -3 . 7310E-01 6.2925E+C0 1.29766+00 1.3671E-01 

4.4105E-03 

WOULD YOU BELIEVE 

8.08708046-05 
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MATRIX OC HAS 1 ROWS AND 15 COLUMNS 

-1.1106E-04 1.3073E-03 -6.5665E-C3 9.1840E-04 5.4555E-03 

-7.5666E-02 2.5109E-04 8.4918E-04 1.2352E-03 -1.0992E-03 

-1.0369E-03 8.6470E-03 6.2564E-05 1.8393E-03 -L.0767E-02 

THE VALUE OF THE FIT ERROR IS 1 .027606-04 FOR THE ITERATION NUMBER 6 
MATRIX EE*D HAS 1 ROWS ANO 7 COLUMNS 
1.1862E-05 2.4446E-06 6.4610E-06 2.3854E-05 5.1881E-05 

4.2667E-06 1.9910E-06 

MATRIX RMS HAS 1 ROWS AND 7 COLUMNS 

3.4441E-03 1.5635E-03 2.5418E-03 4.8841E-03 7.2029E-03 

2.0656E-03 1.4U0E-03 

MATRIX GRAD HAS i ROWS ANO 16 COLUMNS 
1.3120E-01 -1.3166E-02 -1.S781E-C2 6.5262E-03 -4.2880E-04 

-5.4349E-04 -7.8444E-01 7.8726E-C2 -3.6398E-02 1.1753E-03 

-4.3023E-02 7.7035E-03 -3.2960E-C1 1.8323E-01 -3.2993E-04 

1.0276E-04 

WOULD 70U BELIEVE 

8.3062154E-05 

MATRIX OC HAS 1 ROWS ANO 15 COLUMNS 

-8.6889E-05 -1.4371E-04 -1.0037E-03 -9.8075E-05 -1.8201E-04 
-1.7278E-02 3.4672E-05 3.0734E-C4 1.7953E-04 2.7525E-05 

1.89216-04 2.6055E-03 -9.0888E-06 2.8012E-04 -1.5325E-03 

THE VALUE OF THE FIT ERROR IS 2.29906E-05 FOR THE ITERATION NUMBER 7 
MATRIX £E*D HAS 1 ROWS AND 7 COLUMNS 

1.7275E-06 8.4306E-08 1.7895E-06 1.7192E-05 1.8254E-06 

2.2031E-07 1.5132E-07 

MATRIX RMS HAS 1 ROWS AND 7 COLUMNS 

1.3143E-03 2.9036E-04 1.3377E-03 4. 1464E-03 1.3511E-03 

4.6937E-04 3.8900E-04 

MATRIX GRAO HAS 1 ROWS ANO 16 COLUMNS 
7.8532E-03 -1.0569E-03 -5.5453E-04 -4.0673E-04 1.5895E-04 

-2.5042E-04 -3.7550E-02 5.2844E-03 -4.7062E-03 1.9848E-03 

-2.8049E-03 1.2471E-03 -6.3983E-02 1.C544E-02 -1.3571E-05 

2.2991E-05 

WOULD YOU BELIEVE 

8.3433431E-05 

MATRIX OC HAS l ROWS ANO 15 COLUMNS 

-3.6997E-05 -8.3164E-06 -2.436SE-C4 -2.9114E-06 -1.1388E-04 
-3.3925E-03 1.0085E-05 6.0567E-C5 6.0934E-05 1.7531E-05 

4. 8660E-05 5.1169E-04 -4.9164E-06 5.8552E-05 -2.8182E-04 

THE VALUE OF THE FIT ERROR IS 2.CS472E-C5 FOR THE ITERATION NUMBER 8 
MATRIX EE*D HAS 1 ROWS ANO 7 COLUMNS 

1.4983E-06 3.4114E-08 1.5057E-06 1.7077E-05 6.1100E-07 

5.9819E-08 1.6 148E-07 

MATRIX RMS HAS 1 ROWS AND 7 COLUMNS 

1.2241E-03 1.8470E-04 1.2271E-C3 4.1324E-03 7.8166E-04 

2.4456E- 04 4.0184E-04 

MATRIX GRAO HAS 1 ROWS ANO 16 COLUMNS 
5.9765E-C4 -6.2485E-05 9.5939E-06 -8.3012E-05 2.8771E-05 

-5.0517E-05 -1.5591E-04 5.0592E-04 -4.4072E-04 2.8803E-04 

-1.8818E-04 2.4321E-04 -7.6072E-C3 4.7133E-04 9.5277E-06 

2.0947E-05 

WOULD YOU BELIEVE 

8.3486853E-C5 

MATRIX OC HAS 1 ROWS ANO 15 COLUMNS 

-8.9911E-06 -5.2281E-06 -7.C370E-C5 1.0552E-06 -3.1564E-05 

-6.5232E-0A 2.3883E-06 1.3300E-05 1.8493E-05 3.4129E-06 

1.1277E-05 9.5692E-05 -1.4811E-C8 1.0794E-05 -5.7084E-05 

THE FIRST COLUMN IS TIME. THE NEXT 3 COLUMNS ARE CONTROL VARIABLES. THE NEXT 7 COLUMNS ARE RESPONSE VARIABLES. 
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3.5258E-04 
3.5741E-04 
3.4530E-04 
3. 1610E-04 
2.6985E-04 
2.0669E-04 
1.2693E-04 
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-8.0488E-05 
-2.0688E-04 
-3.4734E-04 
5.1043E-02 
2.0323E-01 
1.3980E+00 
4 .55632+00 
8.5764E+00 
1.2367E+01 
1.4567E+01 
1.3846E+01 
L.1626E+01 
9. 3 179E + 00 
6.9347E+00 
4.4896E+00 
1 . 996 3E +00 
-5.3113E-01 
-2.6027E+00 
-3.73785+00 
-4.4330E+00 
-5.18015+00 
-5.9738E+00 
-6. 8078E+00 
-7.6746E+00 
-8. 5660E+00 
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-1. 5266E+01 
- 1 .455 1E + 01 
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-1.27005+01 
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-2.2723E+00 
-2.6998E+00 
-3. 0892E+00 
-3.4379E+00 
-3.7435E+00 
-4.0041E+00 
-4.2183E+00 


5.75636-04 
5.9402E-04 
6.0728E-04 
6. 1537E-04 
6.1829E-04 
6. 1610E-04 
6.0888 E-04 
5.9676E-04 
5.7990E-04 
5.5850E-04 
5. 3279 E-04 
5 .1 208E-02 
2.0 1756-01 
6.1996E-01 
1.4705E+00 
2. 50936+00 
3.4932E+00 
4.11856+00 
4.0875E+00 
3.7135E+00 
3.3095E+00 
2.8790E+00 
2.4258E+00 
1.9537E+00 
1.4667E+00 
1 .07456+00 
8.84886-01 
7.73746-01 
6.1786E-01 
4.1941E-01 
1.8088E-01 
-9.49816-02 
-4.051BE— 01 
-7.46486-01 
- 1. 1 155E+00 
-1.50866+00 
-1.9221E+00 

- 2. 3522E + 00 
-2.7950E+00 
-3.2271E+00 
-3.62576+00 
-3.9882E+00 
-4.31236+00 
-4.5961E+00 
-4.8378E+00 
-5.0363E+00 
-5.1905E+00 
-5.29996+00 
-5.3643E+00 
-5.3838E+00 
-5.3588E+00 
-5.2903E+00 
-5.1792E+00 
-5.0271E+00 
-4.8358E+00 
-4.60716*00 
-4.3436E+00 
-4. 04766+00 
-3.7220E+00 

- 3.36966+00 
•-2.9937E+00 
-2.5974E+00 
-2.1 841E+00 
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03 

-a 


3.40 1.0000E-01 -1.3000E+00 

3.45 1.0000E-01 -1.30C0E+00 

3.50 l.OOOOE-Ol - 1 .30QQE+00 

3.55 1.0000E-O1 -1.3000E+00 

3.60 1.0000E-01 -L.3000E+00 

3.65 l.0000£-01 -1.3000E+00 

3.70 1.0000E-01 -1 .300QE+00 

3.75 1.0000E-01 -1.3000E+00 

3.80 1.0000E-01 -1.3000E+00 

3.85 l.OOOOE-Ol -1.3000E+00 

3.90 l.OOOOE-Ol -1 .3000E+00 

3.95 l.OOOOE-Ol -1.3000E+00 

4.00 l.OOOOE-Ol -1 .3000E+00 

4.05 l.OOOOE-Ol -1 .3000E+00 

4.10 l.OOOOE-Ol -1.3000E+00 

4.15 l.OOOOE-Ol -1.3.000E+00 

4.20 l.OOOOE-Ol -1.30C0E+00 

4.25 l.OOOOE-Ol -1.3000E+00 

4.30 l.OOOOE-Ol -1.3000E+00 

4.35 l.OOOOE-Ol -1.3000E+00 

4.40 l.OOOOE-Ol -1 .3000E+00 

4.45 l.OOOOE-Ol -1.3000E+00 

4.50 l.OOOOE-Ol -1.3000E+00 

4.55 l.OOOOE-Ol -1.3000E+00 

4.60 l.OOOOE-Ol -1 . 3000E+00 

4.65 l.OOOOE-Ol -1.3000E+00 

4.70 l.OOOOE-Ol -1 .3000E+00 

4.75 l.OOOOE-Ol -1.3000E+00 

4.80 l.OOOOE-Ol -1.3000E+00 

4.85 l.OOOOE-Ol -1.3000E+00 

4.90 l.OOOOE-Ol -1.30C0E+00 

4.95 l.OOOOE-Ol -1.3000E+00 

5.00 -2.3500E+00 -1.3000E+00 

5.05 -4. 8000E+00 -1.30G0E+00 

5.10 -7.2500E+00 -1.3000E+00 

5.15 -9. 7000E+00 -1.3000E+00 

5.20 -9.5000E+00 -9.0000E-01 

5.25 -9.3000E+00 -5.0000E-01 

5.30 -7.9000E+00 -5.0000E-01 

5.35 -6.5000E+00 -5.0000E-01 

5.40 -5.1000E+00 -5.0000E-01 

5.45 -3.700 OE+OO -5.0000E-01 

5.50 -2. 3000E+00 -5.0000E-01 

5.55 -9.0000E-01 -5.0000E-01 

5.60 5. OOOOE-Ol -5.0000E-01 

5.65 1.9000E+00 -5. OOOOE-Ol 

5.70 3.3000E+00 -5.00C0E-01 

5.75 4.7000E+00 -5. OOOOE-Ol 

5.80 6. 1QOOE+OQ -5. OOOOE-OL 

5.85 7.5000E+00 -5. OOOOE-Ol 

5.90 8. 9000E+00 -5.0000E-01 

5.95 1.0300E+01 -5. OOOOE-Ol 

6.00 1.04Q0E+01 -5. OOOOE-OL 

6.05 1.0500E+01 -5. OOOOE-Ol 

6.10 1 • 060QE+0 1 -5. OOOOE-Ol 

6.15 1.0700E+01 -5. OOOOE-Ol 

6.20 9. 500QE+00 -5. OOOOE-OL 

6.25 8. JOQOE+OO -5.0000E-01 

6.30 7. 1Q00E+00 -5. OOOOE-Ol 

6.35 5. 9000E+00 -5. OOOOE-Ol 

6.40 4.7000E+00 -5. OOOOE-Ol 

6.45 3.5000E+00 -5. OOOOE-Ol 

6.50 2. 3000E+00 -5. OOOOE-Ol 


l.OOOOE+OO 5.205 1E + 00 2.31E5E + 00 -4.0017E+00 
l.OOOOE+OO 7.2009E+00 2.0961E+00 -4.0787E+00 
l.OOOOE+OO 9.2061E+00 1.8690E+00 -4.1337E+00 
l.OOOOE+OO 1.1208E+01 1.6397E+00 -4.1664E*-00 
l.OOOOEtOO 1.3194E+01 1 .4 105E+00 -4. 1768E+00 
l.OOOOE+OO 1.5153E»01 1.1839E+00 -4.1650E*-00 
1.0000E<-00 1.7071E+01 9.6227E-01 -4.1314E+00 
l.OOOOE+OO 1 . 8936E+0 1 7.4783E-01 -4.0763E+00 
l.OOOOE+OO 2. 0739E+01 5.4281E-01 -4.0004E+00 
l.OOOOE+OO 2.2468E+01 3.4931E-01 -3.9045E+00 
l.OOOOE+OO 2.4113E+01 1.6932E-01 -3.7894E+00 
l.OOOOE+OO 2.5665E+01 4.6996E-03 -3.6563E+00 
l.OOOOE+OO 2.7115E+01 -1.4282E-01 -3.5063E+00 
l.OOOOE+OO 2 . 8456E +0 1 -2.7168E-01 -3.3407E + 00 
l.OOOOE+OO 2.9680E+01 -3.8048E-01 -3.1610E+00 
l.OOOOE+OO 3.0781E+01 -4.6797E-01 -2.9686E+00 
l.OOOOE+OO 3. 1756E+01 -5.3311E-01 -2.7552E+00 
l.OOOOE+OO 3.2599E+01 -5.7504E-01 -2.5523E+00 
l.OOOOE+OO 3.3307E+01 -5.9306E-01 -2.3317E+00 
l.OOOOE+OO 3.3880E+01 -5.8670E-01 -2.1050E+00 
l.OOOOE+OO 3.4315E+01 -5.5567E-01 -1.8742E+00 
l.OOOOE+OO 3.4612E+01 -4.9986E-01 -1.6409E+00 
l.OOOOE+OO 3.4774E+01 -4.1939E-01 -1.4069E+00 
l.OOOOE+OO 3. 480 1E+01 -3.1452E-01 -1.1739E+00 
l.OOOOE+OO 3.4697E+01 -1.8574E-01 -9.4379E-01 
l.OOOOE+OO 3.4466E+01 -3.3681E-02 -7.181LE-01 
l.OOOOE+OO 3.4113E+01 1.4083E-01 -4.9853E-01 
l.OOOOE+OO 3.3642E+01 3.3680E-01 -2.S660E-01 
l.OOOOE+OO 3.3061E+01 5.5311E-01 -8.3799E-02 
l.OOOOE+OO 3.2377E+01 7.8847E-01 1.0849E-01 
l.OOOOE+OO 3.1598E+01 1.0415E+00 2.8896E-01 
l.OOOOE+OO 3.0731E+01 1.3106E+00 4.5645E-01 
l.OOOOE+OO 2.9267E+01 1.6636E+00 6.Q711E-01 
l.COQOE+OO 2.6784E+01 2.1667E+00 7.3170E-01 
l.OOOOE+OO 2.3257E+01 2.8150E+00 8.1853E-01 
l.OOOOE+OO 1.8743E+01 3.6021E+00 8.5630E-01 
l.OOOOE+OO 1. 3870E+01 4.4083E+00 8.3825E-01 
l.OOOOE+OO 9.2564E+00 5.1147E+00 7.6614E-01 
l.OOOOE + OO 5 . 130 LE +00 5.7202E+00 6.4643E-01 
l.OOOOE+OO 1.7129E+00 6.2247E+00 4.8654E-01 
l.OOOOE+OO -9. 8466E-01 6.6252E+00 2.9447E-01 
l.OOOOE+OO -2.95o8E+00 6.9199E+00 7.8374E-02 
l.OOOOE+OO -4.2024E+00 7.1078E+00 -1.5349E-01 
l.OOOOE+OO -4.7252E+00 7.1887E+00 -3.9283E-01 
l.OOOOE+OO -4.5334E+00 7.1635E+00 -6.3137E-01 
l.OOOOE+OO -3.6399E+00 7.0338E+00 -8.6095E-01 
l.OOOOE+OO -2.0619E+00 6.8021E+00 -1.0736E+00 
l.OOOOE+OO 1.79L2E-01 6.4719E+00 -1.2614E+00 
l.OOOOE+OO 3.0577E+00 6.0473E+00 -1.4169E+00 
l.OOOOE+OO 6 • 5446E + 00 5.5331E+00 -1.5330E+00 
l.OOOOE+OO 1 . 0607E + 0 1 4.9349E+00 -1.6028E+00 
l.OOOOE+OO 1.52C8E+01 4.2591E+00 -1.6L99E+00 
l.OOOOE+OO 2 . 0045E+CL 3.5492E+00 -1.5799E+00 
l.OOOOE+OO 2.4815E+01 2.8488E+00 -1.4818E+00 
l.OOOOE+OO 2.9488E+01 2.1640E+00 -1.3265E+00 
l.OOOOE+OO 3.403 IE +01 1.50C7E+00 -1.1155E+00 
l.OOOOE+OO 3.8152E+01 9 . 0 1 34E-0 1 -8 . 51 94E- 01 
l.OOOOE+OO 4. 1562E+0 L 4.0757E-01 -5.4241E-01 
l.OOOOE+OO 4.4248E+01 2.2710E-02 -1.9526E-01 
l.OOOOE+OO 4.6202E+01 -2.5081E-01 1.S099E-01 
l.OOOOE+OO 4.7418E+01 -4.1152E-01 5.7768E-01 
l.OOOOE+OO 4.7900E+01 -4.5890E-01 9.8607E-01 
l.OOOOE+OO 4.7653E+01 -3.9337E-01 1.3975E+00 


1.4342E+01 3.97 35E + 01 -4.3849E+00 -1.7574E+00 

1.4032E+01 4.0052E+01 -4.5030E+00 -1.3206E+00 

1 • 362 IE+01 4.01 15E+0 1 -4.5724E+00 -8.7744E-01 

1 . 31 1 IE +01 3.9925E + 01 -4.5931E + 00 -4.3136E-01 

- 1 . 2501 E+0 1 3.9485E + 01 -4.5654E + 00 1.4084E-02 

-1 . 1792E+0 1 3.8800E+01 -4.4901E+00 4.5541E-01 

- 1 . 0986E+01 3.7877E+01 -4.3685E + 00 8.8923E-01 

-1.0086E+01 3.6725E+01 -4.2018E+00 1.3122E+00 

-9.0936E+00 3.5353E+01 -3.9921E+00 1.7213E+00 

-8.C131E+00 3.3773E+01 -3.7414E+00 2.1133E+00 

-6 • 8482E+00 3.1999E+01 -3.4521E+00 2.4855E+00 

-5. 6033E+00 3.0046E+01 -3.1269E+00 2.8352E+00 

-4.2834E+00 2.7930E+01 -2.7688E+00 3.1600E+00 

-2. 8936E+00 2.5667E+01 -2.3809E+00 3.4576E+00 

-1.4398E+00 2.3276E+01 -1.9667E+00 3.7261E+00 

7. 2295E-02 2.0775E+01 -1.5296E+00 3.9638E+00 

1.6363E+00 1 .8 185E+01 -1.0732E+00 4.1691E+00 

3* 2457E + 00 1.5524E + 01 -6.0139E-01 4.3409E+00 

4. 8939E+00 1.2814E+01 -1.1794E-01 4.4783E+00 

6. 5741E+00 1.0074E+01 3.7327E-01 4.5807E+00 

8 . 2796E+00 7.3245E+00 8.6835E-01 4.6476E+00 

1 • 0003E+0 1 4.5866E+00 1.3634E+00 4.6790E+00 

1.1739E+01 1.8798E+00 1.8547E+00 4.6751E+00 

1 . 3478E+0 1 -7.7663E-01 2.3384E+00 4.6365E+00 

1.5216E+01 -3.3638E+00 2.8109E+00 4.5637E+00 

1 • 6946E+0 l -5.8637E+00 3.2686E+00 4.4579E+00 

1 . 8661E+01 -8.2591E+00 3.7084E+00 4.3202E+00 

2.0355E+01 -1.0534E+01 4.1269E+00 4.1521E+00 

2.2023E+01 -1.2673E+01 4.5212E+00 3.9553E+00 

2 . 3660E+01 -1.4663E+01 4.8886E+00 3.7317E+00 

2.5260E+01 -1.6491E+01 5.2265E+00 3.4834E+00 

2.6B18E+01 -1.8146E+01 5.5328E+00 3.2126E+00 

2.8319E+01 -3.9593E+01 8.5789E+00 2.8107E+00 

2.9721E+01 -6.0442E+01 1.1536E+01 2.1712E+00 

3. 0972E+01 -8.0572E+01 1.4379E+01 1.3015E+00 

3.2022E+01 -9.9863E+01 1.7086E+01 2.1072E-01 

3 . 2837E+0 l -9.4955E+01 1.5140E+01 -9.2989E-01 

3.3414E+01 - 8.95 10E + 01 1.3098E + 01 -1.9509E+00 

3 . 3773E+01 -7.5464E+01 1.H10E+01 -2.8329E+00 

3. 3944E+01 -6.1161E+01 9.0562E+00 -3.5572E+00 

3.3962E+01 -4.6694E+01 6.9557E+00 -4.1197E+00 

3.3862E+01 -3.21 60E+01 4.8260E+00 -4.5179E+00 

3.36B3E+01 - 1.7651E+01 2.6852E+00 -4.7503E+00 

3 . 3459E+0 1 -3.2614E+00 5.5144E-01 -4.8169E+00 

3.3227E+01 1.0916E+01 -1.5576E+00 -4.7189E+00 

3. 3022E+01 2.4793E+01 -3.6245E+00 -4.4589E+00 

3. 2879E+01 3.8283E+01 -5.6324E+00 -4.0406E+00 

3.2832E+01 5 .1 302E+01 -7.5652E+00 -3.4690E+00 

3. 2913E+0 1 6.3772E+01 -9.4075E+00 -2.7501E+00 

3 . 31 53 E+0 1 7.5620E + 01 -1.1145E + 01 -1.8910E+00 

3.3581E+01 8.6777E+01 -1.2763E+01 -8.9988E-01 

3 . 4227 E+0 1 9.7181E+01 -1.4250E + 01 2.1420E-01 

3. 5109E+01 9.61 78E+01 -1.4124E+01 1.3825E+00 

3.6231E+01 9.4531E+01 -1.3873E+01 2.5364E+00 

3. 7589E+01 9.2260E+01 -1.3501E+01 3.6673E+00 

3.917 8E + 01 8.9390E + 01 -1.3013E + 01 4.7667E+00 

4. 0983E+0 1 7.5351E+01 -1.0942E+01 5.7674E+00 

4.2977E+01 6.0989E+01 -8.7948E+00 6.6041E+00 

4.5123E+01 4.6403E+01 -6.589LE+00 7.2719E+00 

4.7385E+01 3. 1694E+01 -4.3444E+00 7.7675E+00 

4. 9727E+0 1 1.6960E+01 -2.0797E+00 8.0890E+00 

5.2111E+01 2.3014E+00 1.8572E-01 8.2359E+00 

5 . 4501E+0 1 -1.2184E+01 2.4332E+00 8.2088E+00 
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05 

00 


6.55 
6.60 

6.65 

6.70 

6.75 
6.80 

6.85 

6.90 

6.95 

7.00 

7.05 
7. 10 

7.15 

7.20 

7.25 

7.30 

7.35 

7.40 

7.45 

7.50 

7.55 

7.60 

7.65 

7.70 

7.75 

7.80 

7.85 

7.90 

7.95 

8.00 

8.05 
8.10 

8.15 

8.20 

8.25 

8.30 

8.35 

8.40 

8.45 

8.50 

8.55 

8.60 

8.65 

8.70 

8.75 

8.80 

8.85 

8.90 

8.95 
9.00 

9.05 
9.10 

9.15 
9.20 

9.25 

9.30 


1. 1000E+00 
1 . 00008—0 1 
1.3000E+00 
2.5000E+00 
3.7000E+00 
3.7000E+00 
3. 700 0E+00 
3.7000E+00 
3.7000E+00 
3 . 7000E + 00 
3.7000E+00 
3. 7000E+00 
3. 7000E+00 
3.7000E+00 
3.7000E+00 
3.7000E+00 
3.7000E+00 

3.7 OOOE+OO 
3.7000E+00 
3. 7000E + 00 
3. 7000E+00 
3.7000E+00 
3 . 7000E+00 
3.7000E+00 
3. 7000E+00 
3.7000E+00 
3.7000E+00 
3.7000E+00 
3. 7000E+00 
3.7000E+00 
3.7000E+00 
3.7000E+00 
3.7000E+00 
3.4429E+00 
3. 1 8576+00 
2.9286E+00 
2. 6714E+00 
2.4143E+00 
2. 1571E+Q0 
1.9000E+00 
1.6429E+00 
1 .3857E+00 
1 • 1286E + 00 
8.7143E-01 
6.1429E-01 
3.5714E-01 
1.0000E-01 
9.5455E-02 
9.0909E-02 
8.6364 E-02 
8. 1818E-02 
7.7273E-02 
7.2727E-02 
6.8182E-02 
-6.3636E-02 
-5.9091E-02 


-5.0000E-01 
-5.0000E-01 
-5.0000E-01 
-5.0000E-01 
-5.0000E-01 
9.0000E-01 
2.30Q0E+00 
3.7000E+00 
5.10C06+00 
5.0900E+00 
5.0800E+00 
5.0700E+00 
5.0600E+00 
5.05G0E+00 
5.0400E+00 
5.0300E+00 
5 .0200E+00 
5.0100E+00 
5.0000E+00 
4.9900E+00 
4.9800E+00 
4.9700E+00 
4.9600E+00 
4.9500E+00 
4.9400E+00 
4.9300E+00 
4.9200E+00 
4. 91006+00 
4.9000E+00 
3.4000E+00 
1.9000E+00 
4.0000E-01 
-1.1000E+00 
-1.0357E+00 
-9.7143E-01 
-9.0714E-01 
-8.4286E-01 
-7.7857E-01 
-7.1429E-01 
-6.50C0E-01 
-5.857 IE-01 
-5.2143E-01 
-4.5714E-01 
-3.9286E-01 
-3.2857E-01 
-2.6429E-01 
-2.0000E-01 
-1.9Q91E-01 
-1.8 182E-01 
-1.7 273E-0 1 
-1.6364E-01 
-1.5455E-01 
-1.4545E-01 
-1 . 3636E-01 
-1.2727E-01 
-1.1818E-01 


1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1 • OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1 .OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
l. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
l.COOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
1. OOOOE+OO 
l.COOOE+OO 


4. 6687E+01 
4.5020E+01 
4.2670E+C1 
3.9663E+01 
3.6027E+01 
3.2185E+01 
2.8549E+01 
2. 5 132E+01 
2. 1937E+01 
l. 8816E + 01 
1.5620E+01 
1.2350E+01 
9.0083E+00 
5-5966E+00 
2. 1196E+00 
-1.4176E+00 
-5.0088E+00 
-8.6466E+00 
- 1 • 2323E + 01 
-1 . 6029E+01 
-1.5755E+01 
-2.3491E+C1 
-2.7225E+01 
-3.0948E+01 
-3.4647E+01 
-3.8310E+01 
-4. 1927E + 01 
-4.5485E+01 
-4.8972E + 01 
-5.2531E+01 
-5.6297E+01 
-6.0242E+01 
-6.43316+01 
-6.8307E+01 
-7.1915E+01 
-7.5123E+01 
-7.7900E+01 
-8.0222E + 01 
-8.2066E+01 
—8. 34156+01 
-8.42586+01 
-8.45866+01 
-8.4394E+01 
-8. 3683E+01 
-8.2459E+01 
-8.0730E+01 
-7.8509E+01 
-7.5870E+01 
-7.2891E+01 
-6.9591E+01 
-6.5995E+01 
-6.2129E+01 
-5.3019E+01 
-5.3696E+01 
-4.9190E+01 
-4.4533E+01 


-2. 1630E-0 1 
7.0027E-02 
'♦.62476-01 
9. 57026-01 
1 . 5489E+00 
2. C676E+00 
2.3450E+00 
2. 3806E+00 
2.1754E+00 
1.8637E+00 
1.5785E+00 
1. 3199E+00 
1.0873E+00 
8.8018E-01 
6.9766E-01 
5.3868E-01 
4. 0195E-0 1 
2.8602E-01 
1. 8925E-0 1 
1.0985E-01 
4.5884E-02 
-4.6973E-03 
-4. 40546—02 
-7 .44266-02 
-9.8120E-02 
-1 . 1748E-01 
-1.3489E-01 
-1.5273E-01 
-1.7335E— 01 
-5.9847E-02 
3.2248E-01 
9.6614E-01 
1.8622E+00 
2. 8470E+00 
3.7582E+00 
4. 5891E+00 
5.3339E+00 
5.9875E+00 
6.5455E+00 
7.0044E+00 
7.3615E+00 
7.6149E+00 
7.7636E+00 
7. 8075E+00 
7.7470E+00 
7. 5837E+00 
7.3198E+00 
6.9706E+00 
6.5519E+00 
6.0674E+00 
5.5216E+00 
4.9191E+00 
4.2652E+00 
3.5652E+00 
2.8251E+00 
2.0507E+00 


1.8032E+00 
2. 1947E + 00 
2.5638E+00 
2.9024E+00 
3.2028E+00 
3.4626E+00 
3.6894E+00 
3.8961E+00 
4.Q951E+00 
4.295 IE+00 
4.4974E+00 
4.6998E+00 
4.9001E+00 
5.0960E+00 
5.2853E+00 
5.4662E+00 
5.6366E+00 
5.7947E+00 
5.5389E+00 
6.0675E+00 
6. 1793E+00 
6.2729E+00 
6.3473E+00 
6.4015E+00 
6.4350E+00 
6.4470E+00 
6.4372E+00 
6.4054E+00 
6.3517E+00 
6.2724E+00 
6. 15676+00 
5.9903E+00 
5.7599E+00 
5.4572E+00 
5.0832E+00 
4.6437E+00 
4.1450E+00 
3.59366+00 
2 .99656 + 00 
2.3608E+00 
1.6940E+00 
1.0036E+00 
2.9723E-01 
-4.1755E-01 
-1.1331E+00 
- 1 . 8420E+00 
-2.5368E+00 
-3.2109E+00 
-3. 8 586E + 00 
-4.4752E+00 
-5.0559E+00 
-5.5966E+00 
-6.0934E+00 
-6.5427E+00 
-6.9417E+00 
-7.2875E+00 


5.6860E+01 
S.9154E+01 
6. 1347E+01 
6. 3406E+01 
6.5298E+01 
6. 7004E+01 
6.8523E+01 
6.9865E+01 
7 • 1042E+0 1 
7.2061E+01 
7. 2922E+01 
7.3622E+01 
7.4156E+01 
7.4521E+01 
7. 4715E+0 1 
7.4732E+01 
7.4572E+01 
7.4231E+01 
7. 37 07E + 01 
7.2998E+01 
7.2103E+01 
7. 1022E+01 
6. 9754E+01 
6 • 8300E+0 1 
6 . 6660E+0 1 
6.4836E+01 
6.2830E+01 
6. 0644E+01 
5. 8282E+01 
5. 5744E+01 
5. 3023E+01 
5. 0109E+0 1 
4.6994E+01 
4.3677E+01 
4. 0170E+0 1 
3.6493E+01 
3.2666E+01 
2.8712E+01 
2.4653E+01 
2.0514E+01 
1.6321E+01 
1 • 2098E+0 1 
7. 8721E+00 
3 • 66 85E+00 
-4. 8669E-0 1 
-4 • 5680E+00 
-8. 5505E + 00 
-1.24 HE +01 
-1 . 61 32E+01 
-1 . 9695E+01 
—2 • 3086E+01 
-2.6290E+01 
-2. 9295E+01 
-3.2088E+01 
-3.4661E+01 
-3.7005E+01 


-2.6402E+01 
-4.0259E+01 
-5.3666E+01 
-6.6539E+01 
-7.8758E + 01 
-7.4826E+01 
-7.0551E+01 
-6.6120E+01 
-6.1678E+01 
-6 .31 64E+01 
-6.4656E+01 
-6.6125E+01 
-6.7549E+01 
-6.8902E+01 
-7.0162E+01 
-7.1306E+01 
-7.2316E+01 
-7.31 71E + 01 
-7.3855E+01 
-7.4353E+01 
-7.4652E+01 
-7.4741E + 01 
-7 .46 10E+01 
-7.4253E + 01 
-7 .3664E+01 
-7.2840E + 01 
-7.1780E+01 
-7.0487E + 01 
-6.8962E+01 
-7.3342E+01 
-7.7232E+01 
-8.0479E + 01 
-8.2942E+0L 
-7 .59 57E+01 
T-6.8264E + 01 
-5 .99 40E+01 
-5.1063E + 01 
-4.17 17E+01 
-3.1985E+01 
-2.1957E+01 
-1.1720E + 01 
-1.3598E+00 
9.0330E+00 
1.9373E+01 
2.9575E+01 
3.9555E+01 
4.9234E+01 
5.6250E+01 
6.2867E+01 
6.9037E+01 
7.4713E+01 
7.9855E+01 
8.4427E+01 
8.8397E+01 
9.17396+01 
9.4432E+01 


4.6443E+00 
6.8009E+00 
8 . 8859E+00 
1.0882E+01 
1.2775E+01 
7.9599E+00 
3* 1 238E+00 
-1.7056E+00 
-6.5013E+00 
-5.9689E+00 
-5.4370E+00 
-4.9105E+00 
-4. 3944E+00 
-3.8931E + 00 
-3.411 IE +00 
-2.952SE+00 
-2.5214E+00 
-2. 1213E+00 
-1 . 7555E + 00 
-1.4270E+00 
-1.1384E+00 
-8.92 00E-01 
-6 • B966E-01 
-5.3285E-01 
-4. 2266E-0 1 
-3. 5978E-01 
-3.4449E-01 
-3.7670E— 01 
-4.5590E-01 
4. 9858E+00 
1.0292E+01 
1 . 5434E+01 
2.0384E+01 
1.8982E+01 

1.7 44 1E+01 
1.5775E+01 
1. 3999E+01 
1.2128E+01 
1. 01 79E+01 
8. 1672E+00 
6. 1 097E+00 
4.0232E+00 
1.9244E+00 
-1.6998E-01 
-2.2437E+00 
-4.2809E+00 
-6.2660E+00 
-7.6920E+00 
-9.0460E+00 
-1.0318E+01 
-l. 1499E+01 
-1.2582E+01 
-1.3558E+01 
-1.4421E+01 
-1.5166E+01 
-1.5789E+01 


8.0100E+00 
7.6429E+00 
7.1 122E+00 
6.4239E+00 
5.5851E+00 
4.8002E+00 
4.2694E+00 
3.9910E+00 
3.9611 E+00 
4.0313E+00 
4.0554E+00 
4.0342E+00 
3.9684E+00 
3.8595E+00 
3.7091E+00 
3. 51896+00 
3.2912E+00 
3.0285E+00 
2.7334E+00 
2.4088E+00 
2.0577E+00 
1.6835E+00 
1. 28946+00 

8. 7 894E-01 
4. 55726-01 
2.3321E-02 
-4.1464E-01 
-8.5452E-01 
- 1.2928E+00 
-1.8769E+00 
-2. 7501E+00 
-3.B987E+00 
-5.3073E+00 
-6.7885E+00 
-8.15B0E+00 
-9.406 9E+00 
-1.0527E+01 

- 1.1512E+01 
-1.2355E+01 
- 1.3052E+01 
-1.3600E+01 
- 1.3996E+01 
-1.4240E+01 
- 1.4332E+01 
-1.4272E+01 
- 1.4064E+01 
-1.3711E+01 
-1.3235E+01 
-1.2659E+01 
- 1 . 1987E+01 
-1.1228E+01 
-1.0387E+01 
-9.4726E+00 

- 8.4 927E+00 
-7.4559E+00 
-6.3710E+00 
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9.35 -5.4545E-02 -1.0909E-01 
9.40 -5.0000E-02 -1.0000E-01 
9.45 -4.5455E-02 -9.0909E-02 
9.50 -4.0909E-02 -8.1818E-02 
9.55 -3.6364E-02 -7.2727E-02 
9.60 -3.1818E-02 -6.3636E-02 
9.65 -2.7273E-02 -5.4545E-02 
9.70 -2.2727E-02 -4.5455E-02 
9.75 -1.8182E-02 -3.6364E-02 
9.80 -1.3636E-02 -2.7273E-02 
9.85 -9.0909E-03 -1.8182E-02 
9.90 -4.5455E-03 -9.0909E-03 
9.95 -5.6932E-13 -1.1386E-12 
THE VALUE OF THE FIT ERROR IS 
MATRIX EE*D HAS 1 ROWS AND 
1.4898E-06 3.8844E-08 1.4690E-06 

4.7047E-08 1.8428E-07 

MATRIX RMS HAS 1 ROWS AND 7 COLUMNS 

1.2206E— 03 1.9709E-04 1.2120E-C3 4.1330E-03 

2.1690E-C4 4.2928E-04 

MATRIX GRAD HAS 1 ROWS AND 16 COLUMNS 
9.1397E-05 -9.2442E-06 1.4474E-C5 -1.9878E-05 

-9.6792E-06 4.4500E-04 9.2840E-05 -5.0273E-05 

-1.2587E-05 4.5126E-05 -1.4U1E-C3 -7.6391E-07 

2.0880E-05 

WOULD YOU BELIEVE 

8.3496613E-05 

MATRIX OC HAS 1 ROWS AND 15 COLUMNS 

-1.9149E-06 -1.8729E-06 -1.7436E-C5 9.9519E-08 -7.6805E-06 

-1.2673E-0A 4 .9877E-07 2.9333E-06 4.5277E-0° 6.6766E-07 

l . 81646-05 -3.6874E-C7 1.8826E-06 -1.1702E-05 

HAS 4 ROWS and 4 COLUMNS 
6.3028E-02 -1.14U0E+C1 -0. 

-2.0814E-01 2.1701E»C0 -0. 

-1. OOOOE+OO -3. 14946-Cl 7.8900E-02 

- 0 . - 0 - 

4 ROWS AND 3 COLUMNS 


1. OOOOE+OO -3.9758E+01 1.2483E+00 -7.5781E+00 -3.9113E + 01 9.6462E+01 -1.6285E+01 

1. OOOOE+OO -3.4898E+01 4.2431E-01 -7.8117E+00 -4.0979E+01 9.7817E+01 -1.6654E+01 

1. OOOOE+OO -2.9988E+01 -4.1489E-01 -7.5871E+Q0 -4.2602E+01 9.8493E+Q1 -1.6893E+01 

1. OOOOE+OO -2.5060E+01 -1.2628E+00 -8.1034E+00 -4.3978E+01 9.8492E+01 -1.7002E + 01 

1. OOOOE+OO -2.0I50E+01 -2.1129E+00 -8.1605E+00 -4.5I08E+01 9.7819E+01 -1.6981E+01 

l.OOOQE+OO -1.5289E+01 -2.9588E+00 -8.1584E+00 -4.5993E+0I 9.6485E+01 -1.6834E+01 

1. OOOOE+OO -1.0512E+01 -3.7942E+00 -8.0977E+00 -4.6638E+01 9.4507E+01 -1.6561E+01 

1. OOOOE+OO -5.8489E+00 -4.6129E+00 -7.9795E+00 -4.7047E+01 9.1906E+0I -1.6I68E+01 

1. OOOOE+OO -1 .331 IE +00 -5.4090E+00 -7.8053E+00 -4.7225E+0I 8.8708E+01 -1.5657E+01 

1. OOOOE+OO 3.012 5E + 00 -6.1768E+00 -7.5771E+00 -4.7I82E+01 8.4943E+01 -I.5036E + 01 

1. OOOOE+OO 7. 154 3E + 00 -6.9108E+00 -7.2971E + 00 -4.6927E+01 8.0645E+01 -1.43I0E + 01 

1. OOOOE+OO I.1069E+01 -7.6061E+00 -6.9680E+00 -4.6471E+01 7.5853E+01 -1.3485E+01 

I. OOOOE+OO I.4732E+01 -8.2578E+00 -6.5928E+00 -4.5825E+01 7.0606E+01 -I.2570E+01 

2.08800E-05 FOR THE ITERATION NUMBER 9 
7 COLUMNS 


1.7081E-05 5.6953E-07 


7.5467E-04 


9 • 3030E-06 
7.7592E-05 
3.0171E-C7 


2.6512E-06 
MATRIX A 
-2.7396E-C1 
1.3100E-03 
8.2997E-02 
1 .OOOOE+OO -0. 
MATRIX B HAS 


8.2232E+00 4.1772E+00 1.9426E-C3 

— L. 1406E + 00 -3.7614E+00 -S.0107E-04 


- 0 . 

- 0 . 

MATRIX DZl 

0. 

0. 

MATRIX CZN 

0. 

0. 

MATRIX F 
-2.7396E-01 


- 0 . 

- 0 . 

HAS 

0. 

0. 

HAS 

0. 

0. 

HAS 


3.5644E-C4 

- 0 . 

ROWS AND 

0 . 


COLUMNS 

0 . 


I ROWS ANO 

0 . 


COLUMNS 

0 . 


3 ROWS ANO 4 CCLUMNS 
6.3028E-02 -1.1400S+0I 0. 

1.3100E-03 -2.0814E-0I 2.17C1E+C0 0. 

8.29S7E-02 -l. OOOOE+OO -3.1494E-C1 7.8900E-02 

MATRIX G HAS 3 ROWS ANO 3 COLUMNS 
8.2232E+00 4.1772E+00 1.9426E-C3 

-1.1406E+00 -3.7614E+00 -S.01C7E-C4 
0. 0. 3.5644E-04 


OS 

to 


8.3496613E-05 

MATRIX S-IS HAS 15 ROWS ANO 15 COLUMNS 
1. OOOOE+OO -1.7053E-13 -S.3814E-15 8.8818E-15 1.0214E-14 

-1.2790E-13 -3.7517E-12 -1.0658E-12 3.6948E-13 -2.4025E-13 

3.0909E-13 -7.9561E-13 2.4727E-12 1.0658E-14 1.9185E-13 

-2.7569E-12 1. OOOOE+OO -2.9843E-13 8.1712E-14 3.3751E-14 


5.2471E+00 
4.0935E+00 
2.9195E+00 
1.7345E+00 
5.4806E-01 
6.3068E-01 
1.7926E+00 
2.9289E+00 
4.0310E+00 
5.0908E+00 
6. 1005E+00 
7.0531E+00 
7.9417E +00 
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-3 

O 


-6.8212E-13 
-4. 26336-14 
-7.84446-12 
-1.1227E-12 
-3.97906-13 
5.22256-13 
-9.2371E-14 
-2.5668E-13 
-6.48016-12 
-4.54756-13 
1.5632E— 12 
-3.50166-11 
1. 00006*00 
9. 32236-12 
3.00206-13 
2.93106-14 
-8.0824E-14 
1.25066-12 
2.13166-13 
-3.05536-13 
1.79066-12 
2.4869E-13 
— 3.4106E-13 
1.06586-13 
1.77646-14 
6.7502E- 14 
-5.57786-13 
-7.77166-14 
1.00006*00 
4.43386-12 
1.04456-12 
-1.88296-12 
-3.10866-14 
-1.11026-14 
3.3751E-14 
5.54226-13 
4.61856-14 
-2.4514E-13 
-5.11596-12 
-7.38966-13 
1.36426-12 
MATRIX DC* 
-1.91496-06 
-1.87296-06 
— 1.74 36E-05 
9.95196-08 
-7.68056-06 
-1.2673E-C4 
4.98776-07 
2.93336-06 
4.52776-06 
6.6766E-C7 
2.65126-06 
1.81646-05 
-3.68746-07 
1.88266-06 
-1.17026-05 


.29656-11 - 
•67 16t-l2 - 
.17956-12 
.32146-11 - 
.88856-12 
.02876-13 - 
.73596-12 - 
.50156-13 - 
.24236-13 - 
.63716-11 - 
. 29696-12 
.91176-12 
.45526-10 - 
.45566-11 
. 01986-14 - 
.00006+00 
. 3856E-13 - 
.06066-13 
.41286-1 2 
.02326-12 - 
.62906-13 
.77716-12 
.44956-12 - 
. 1316E-14 
.41066-13 
.26336-14 - 
-6956E-13 
.82656-12 - 
.76596-13 - 
.21506-12 - 
.63756-11 
.00006+00 - 
.06586-14 - 
.03746-12 - 
.92786-14 
.90806-14 - 
.79806-12 
.84746-13 
.68436-13 
•0009E-1 1 - 
.41066-12 
S 15 ROWS 


2.78536 
4.54756 
1.00006 
4.32016 
4.5020E 
8.74866 
1.6698E 
8.4270E 
8.43776 
4.09276 
1.65966 
5.79056 
3.27426 
7.82176 
4.44096 
2.27376 
1.05876 
5.1926E 
1 .00006 
4 • 774eE 
1.86306 
1.5916E 
1.26156 
3.93026 
1.27906 
9.66346 
5.43336 
8.56206 
1.70466 
1.76536 
7.84446 
6.36656 
2.60626 
7 .99366 
l.OOOOE 
1.1058E 
6.25286 
5.9686E 
2.3981E 
4.77486 
1.54616 
ANO 


-12 1 
-12 4 

+ 00 7 

-12 3 

-11 1 
-14 1 

-13 -3 
-12 1 
-14 l 
-12 1 
-11 5 

-13 3 

-11 8 
-11 2 
-16 -5 
- 13 -7 
-12 3 

-14 -4 
+ 00 -3 
-12 -6 
i-13 -8 
-12 1 
-11 -9 
-14 -2 
-13 5 

-13 -8 
-13 -2 
-13 -3 
-11 -1 
-13 -7 
-12 -1 
-12 -3 
-14 1 

-14 4 

+ 00 1 
-13 4 

-13 -8 
-13 l 
-14 7 

-12 1 
-11 -5 
1 CCL 


•5632E-13 
.67546-12 
.99366-13 
.83696-12 
. 19666-11 
.00006+00 
•8192E-13 
.37676-12 
.27906-13 
• 3927E-12 
.68436-14 
.69486-13 
. 1855E-12 
. 273 7E-13 
.55116-15 
. 1942E- 14 
.10866-14 
♦7962E-14 
. 33966- 13 
. 1817E-13 
.17126-14 
.00006+00 
.16606-13 
.57576-14 
.329 16-14 
•5265E-14 
.62686-13 
.73036-13 
.21866-12 
.99366-14 
.02326-12 
.83696-13 
.21576-14 
.66296-14 
.42UE-13 
.04126-14 
.17126-14 
. OOOOE+OO 
.81606-14 
.36426-12 
.68436-14 
UMNS 


-3. 1442E-13 
5.40016-13 
-4.54756-13 
-3.99336-12 
1.76216-12 
2.37596-14 
3.38406-13 
1.72316-13 
1.00006+00 
-1 . 1280E-12 
7.38966-13 
8.52656-13 
-2.48696-12 
1.09146-11 
-7.77166-16 
5.58446-14 
-4.61856-14 
2. 66456-15 
3.89476-13 
-2.98436-13 
2.1316E-14 
3.9702E-13 
-3.83696-13 
1.77646-14 
1 . OOOOE+OC 
-7.10546-15 
2.13386-13 
8.9195E-13 
1. 0125E-13 
-1. 86 52 E- 13 
9.85886-14 
- 1 . 6485E-12 
-6.43936-15 
-5.76216-14 
1.95406-14 
-1.19906-14 
7 .99366-15 
-8. 52656-14 
4.97386-14 
-6.14626-13 
l.OOCOE+OO 
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Parameter Group Row 

l x i 

-0 5 1 

1 1 1 

-0 5 1 

1 1 1 

-0 5 1 

1 2 1 

-0 6 1 

1 2 1 

-0 6 1 

1 2 1 

-0 6 1 

1 1 2 

-0 5 2 

1 1 2 

-0 5 2 

1 1 2 

-0 5 2 

12 2 

-0 6 2 

12 2 

-0 6 2 

1 2 2 

-0 6 2 

l 1 3 

-0 5 3 

1 1 3 

-0 5 3 

1 2 3 

-0 6 3 


Column A priori Deviation Scale 

1 -0.000099999.0000 1.0000 

1 -0.0000 -0.0000 1.0000 

2 -0.000099999.0000 1.0000 

2 -0.0000 -0.0000 1.0000 

3 -0.000099999.0000 1.0000 

3 -0.0000 -0.0000 1.0000 

1 -0.000099999.0000 1.0000 

1 -0.0000 -0.00 00 1.0000 

2 -0.000099999.0000 1.0000 

2 -0.0000 -0.0000 1.0000 

3 -0.000099999.0000 1.0000 

3 -0.0000 -0.0000 1.0000 

1 -0.000099999.0000 1.0000 

1 -0.0000 -0.0000 1.0000 

2 -0.000099999.0000 1.0000 

2 -0.0000 -0.0000 1.0000 

3 -0. 000099999.0000 1.0000 

3 -0.0000 -Q.000C 1.0000 

1 -0.000099999.0000 1.0000 

1 -0.0000 -0.0000 1.0000 

2 -0.000099999.0000 1.0000 

2 -0.0000 -0.0000 1.0000 

3 -0.000099999.0000 1.0000 

3 -0.0000 -0.0000 1.0000 

l -0. 000099999. 000C 1.0000 

1 -0.0000 -0.0000 1.0000 

3 -0.000099999.0000 1.0000 

3 -0.0000 -0.0000 1.0000 

3 -0.000099999.0000 1.0000 

3 -0.0000 -0.0000 1.0000 


03/13/72 LRC CM SCOPE 3.0 660CB-131K 0 9/01/710 

13. 01. 09. LB9 16 74* - 0313 1251 

13.01.09. LRC COMPUTER COMPLEX 

13. 01. 09. JOB, 1,300, 77777. 1000. A2329 R 

13.01. 09. KN211 1150 

13. 01. 09. USER. TAYLOR, LAWRENCE M , JR 0008 

13.01. 09. 75539N 14110 

13.01.09. RUNIS) 

13.01.16. LINECNT( 10000) 

13.01. 16. SETINDF. 

13. 01. 17. LOO. 

13. 01. 25. MEMORY 035600 CP 

13.01.25. RFL 0000085 0/S CALLS 

13. 01.25. RFL CPU 3.589856 SEC. 

13. 03. 05. STOP 

13. 03. 05. SPPRl NT (OUTPUT, 5) 

13.03.06. 0000103 0/S CALLS 

13. 03. 06. CPU 101.891553 SEC. 

13.03.07. PPU 15.654912 SEC. 

13.08.29. LB91674. 1813 LINES FRINTEO. LP25 


-0 


Answer Deviation 


-.2740 .0975 

- 0.0000 0.0000 

.0630 .0139 

-O.DOOO 0.0000 

-11.4004 .0065 

- 0.0000 0.0000 

8.2232 .0259 

-0.0000 0.0000 

4.1772 .0118 

- 0.0000 0.0000 

.0019 .0037 

- 0.0000 0.0000 

.0013 .3813 

-0.0000 0.0000 

-.2081 .0678 

-0.0000 0.0000 

2.1701 .0355 

-0.0000 0.0000 

-1.1406 .0631 

-0.0000 0.0000 

-3.7614 .0251 

-0.0000 0.0000 

-.0009 .0240 

-0.0000 0.0000 

.0830 .3457 

-0.0000 0.0000 

-.3149 .0516 

-0.0000 0.0000 

.0004 .0221 

- 0.0000 0.0000 


APPENDIX B - Concluded 
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